Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 777,282 A→T 6.8% intergenic (+174/‑259) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273777,2820AT 6.8% 249.6 / 5.5 98intergenic (+174/‑259)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (39/50);  new base T (3/4);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGA  >  CP009273/777134‑777422
                                                                                                                                                    |                                                                                                                                            
ttttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTa                                                                                                                                              >  1:943764/1‑149 (MQ=32)
 tttttGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                             <  1:634999/149‑1 (MQ=32)
        tATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGGTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCT                                                                                                                                      >  2:692105‑M1/1‑149 (MQ=255)
        tATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGTGAGTAATCCTGCTCGCACCACCAATGTAAAAAAGCAACCTAAAAGCGCT                                                                                                                                      >  1:165751‑M1/1‑149 (MQ=255)
        tATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCCCGACCCACCAATGTAAACAAGCCCCCTAAAGGCGCT                                                                                                                                      >  1:45069‑M1/1‑149 (MQ=255)
        tATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCT                                                                                                                                      >  1:1114198‑M1/1‑149 (MQ=255)
         aTCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTT                                                                                                                                     >  2:207219‑M1/1‑149 (MQ=255)
          tCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTT                                                                                                                                    <  1:1217546‑M1/149‑1 (MQ=255)
          tCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCTTT                                                                                                                                    >  2:1063447‑M1/1‑149 (MQ=255)
            tGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                  >  1:938576‑M1/1‑149 (MQ=255)
            tGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                  >  2:1005570‑M1/1‑149 (MQ=255)
            tGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                  >  2:194382‑M1/1‑149 (MQ=255)
            tGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTT                                                                                                                                  >  1:162782‑M1/1‑149 (MQ=255)
             gCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTT                                                                                                                                 <  2:32258‑M1/149‑1 (MQ=255)
             gCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCCTATTT                                                                                                                                 >  2:935768‑M1/1‑149 (MQ=255)
                    tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGACATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGCCCCTTTTTTACTATTc                                                                                                                          >  2:383892/1‑149 (MQ=9)
                     cAAAGATTCGCACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGCCGCCGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                                                                                                                         <  1:289496/149‑1 (MQ=12)
                     cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                                                                                                                         <  2:1023133/149‑1 (MQ=32)
                     cAAACATTCGCAGCTGCAGCAGGTTTTAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCa                                                                                                                         <  1:1090575/149‑1 (MQ=255)
                      aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAg                                                                                                                        >  2:385073/1‑149 (MQ=32)
                           tttgtaCCCGCAGCAGGTTTGGGTTGAGCGCAGGGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATc                                                                                                                   <  2:314910/144‑1 (MQ=255)
                           ttCGAACCTGACGGAGGTTTGAGTTGAGCGCAGCCAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACTCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATc                                                                                                                   <  2:815303/149‑1 (MQ=11)
                           tcccaaCCTACAGCAATTTTCAGTCGAGCGCACCGAAACAACGGAGCCCCACGCGGCCACGGCCCAAAGAGCGATCGAAGCTAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCATGCATc                                                                                                                   <  2:521586/145‑1 (MQ=11)
                              gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTc                                                                                                                >  1:1209755/1‑149 (MQ=39)
                              gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCAGCCTc                                                                                                                >  1:628386/1‑149 (MQ=37)
                                   tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAAttt                                                                                                           <  2:922917/149‑1 (MQ=255)
                                      cgcagGTTTGGGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAc                                                                                                        <  1:249039/148‑1 (MQ=255)
                                      agcagTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCCAAAGGCGCTTTTTTACTATTCAGGCATCCTCAATTTCAc                                                                                                        >  2:868637/1‑149 (MQ=255)
                                        cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAGTTTCACtt                                                                                                      >  1:194346/1‑149 (MQ=255)
                                        cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACtt                                                                                                      >  1:1117339/1‑149 (MQ=255)
                                          ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCTTTTTTacattggtgggtcgtgcaggatgactcgc                                                                                                    <  1:1063447‑M1/149‑30 (MQ=255)
                                             ttgagttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTaa                                                                                                 >  1:319575/1‑149 (MQ=255)
                                               gggttgggCGCGGGGAAACAACGGAGCCGGTCAGGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCTCCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                               <  2:697290/147‑1 (MQ=17)
                                               gggttgagCGCCGCGCAACAAGGGCGCCGCTCGCGGCGACGGCTCGAACGGCGAGCGACGCGTCTATGCCTGTATTACCCACCAATGTAAAAAAGCGCCCTAACCGCGCTTTCTTGCTATTCAGGCATCCTTAATTTCACtttttccac                                                                                               <  2:571627/147‑7 (MQ=9)
                                               gagttgcgCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                               <  2:735217/149‑1 (MQ=255)
                                               gagttgagCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGAGCGAGCGAAGCGACTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                               <  2:886736/149‑1 (MQ=255)
                                               gagttgagCGCAGCGAAACAACGGCGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAAc                                                                                               <  2:894774/149‑1 (MQ=255)
                                                agttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCTTTTTTacattggtgggtcgtgcaggatgactcgcttcgct                                                                                              <  2:1103703‑M1/149‑36 (MQ=255)
                                                 gttgagCGCAGCGAAACAACGGAGCCGCTAGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCt                                                                                             >  1:556480/1‑149 (MQ=255)
                                                 ggtgagcgGAGCGAAGCAAGGGGGCCGGTCGCGGCGACGGCCCGATGGGCGGGGGAGGCGAGTCATCCTGCAGGACCCACGTATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTacattggtgggtcgtgcaggatgactcgcttcgctc                                                                                             <  2:1084301‑M1/147‑37 (MQ=255)
                                                  ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGTGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtg                                                                                            >  1:550191/1‑149 (MQ=255)
                                                  ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTAACTTTGTAAACCtg                                                                                            >  1:582197/1‑149 (MQ=255)
                                                  ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCCCTTTGTCAACCtc                                                                                            >  1:618843/1‑148 (MQ=255)
                                                  ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtg                                                                                            >  1:489856/1‑149 (MQ=255)
                                                  ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtg                                                                                            >  1:458886/1‑149 (MQ=255)
                                                   tgagCTCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGTGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAACGCGCCTTTAGGGCGCTTTTTTacattggtgggtcgtgcaggatgactcgcttcgctcgc                                                                                           <  1:1121106‑M1/149‑39 (MQ=255)
                                                    gagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgat                                                                                          <  2:331148/149‑1 (MQ=255)
                                                      gcgcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCtgctga                                                                                        >  1:103728/1‑149 (MQ=255)
                                                       cgcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGAc                                                                                       >  1:1128919/1‑149 (MQ=255)
                                                        gcAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACa                                                                                      <  2:497805/149‑1 (MQ=255)
                                                              caacaccGGAGCCGCTCGTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGCCa                                                                                <  2:1038066/148‑1 (MQ=255)
                                                                acaacGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCaga                                                                              <  1:385073/149‑1 (MQ=255)
                                                                 caacGGTGCCGCTCGCGGCGGCGGACCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCagag                                                                             <  1:553121/149‑1 (MQ=255)
                                                                 caacGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCagag                                                                             <  2:336256/149‑1 (MQ=255)
                                                                     ggAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACCTTGTAAACCTGATGACATCGTCAGAGCTTa                                                                         <  1:650342/149‑1 (MQ=255)
                                                                     ggAGCCGCCCGCGGCGGCGCCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACCGCCCTTTGTTTCTATTATGGGATCCTCAATTTCACTTTGTAAACCTTATGACATCGTCAGAGCTTa                                                                         <  1:601591/149‑1 (MQ=18)
                                                                      gtgCCGTTCGCGGCGCCGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTAc                                                                        <  2:173923/147‑1 (MQ=255)
                                                                      gAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTAc                                                                        >  1:983307/1‑149 (MQ=255)
                                                                       aGCCGCTCGCGGCGCCGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGAGGACATCGTCAGTGCTTACt                                                                       <  2:843255/149‑1 (MQ=255)
                                                                             tCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaa                                                                 <  2:642446/149‑1 (MQ=255)
                                                                              cgcgGCGACGGCCCGAAGGGCGAGCGAAGCTAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag                                                                <  1:500528/149‑1 (MQ=255)
                                                                                cggcgACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagca                                                              >  1:1156430/1‑149 (MQ=255)
                                                                                   cgACGGCCCGAAGGGCGAGCTAAACTACTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAAct                                                           <  2:1178495/149‑1 (MQ=255)
                                                                                          ccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATTCCGGAATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAAGCCTTACTGTGCAAGCAACTCTATGTc                                                    >  1:289881/1‑149 (MQ=25)
                                                                                           cGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCg                                                   <  2:372798/149‑1 (MQ=255)
                                                                                            gAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCgg                                                  <  1:655541/149‑1 (MQ=255)
                                                                                             aaGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGt                                                 <  2:279349/149‑1 (MQ=255)
                                                                                              aGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTg                                                <  2:559777/149‑1 (MQ=255)
                                                                                                 gcgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGaa                                             <  1:328722/149‑1 (MQ=255)
                                                                                                 gcgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGaa                                             >  2:1154635/1‑149 (MQ=255)
                                                                                                 gcgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGaa                                             <  2:556480/149‑1 (MQ=255)
                                                                                                 gcgagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAAGCGCTTTTTTGATATTCAGGCATCCTCAATTTACCTTTGTAAACCTGATCAAATCGTCAAAGCTTCCTGTGACGGCAAATCTATGTCGGTGGaa                                             >  1:723852/1‑149 (MQ=21)
                                                                                                 gcgagcgaAGCGAGTCATCCGGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTCTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGaa                                             <  2:556511/149‑1 (MQ=255)
                                                                                                  cgagggaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAAt                                            <  2:1209811/149‑1 (MQ=255)
                                                                                                   gagcgaAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACACTCTGTCGCTGGAAtt                                           >  2:730613/1‑149 (MQ=255)
                                                                                                    agcgaagcgaGTCATCCTGCACGACCCACCCCTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCTTCCCCAAGTTCACTTTGTAAACCTCATGAACTGGTCAGTGCTTACTGTGCAAGCAACTCTATGTCGATGGCAtga                                          >  2:619741/1‑147 (MQ=21)
                                                                                                      cgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAGCGCTAAGTCCGGGGAATTAgg                                        >  2:967992/1‑149 (MQ=37)
                                                                                                      cgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCCTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAgg                                        <  2:619935/149‑1 (MQ=255)
                                                                                                               tCATCCTGCACGACCCACCAATGTACCAAGGAAACTTTTAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCAATTTGTAAACCTGATGAAATCGTCAGAGCGTACTGTGCAACCAGCTCTCCGTCGGTGTGAGTAGGCGTAAAATg                               >  1:299396/1‑149 (MQ=255)
                                                                                                               tCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTATACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTGGGCGTAAAATg                               <  2:147059/149‑1 (MQ=255)
                                                                                                                cATCCTGCACGACCCGCCAATGTAATAAAGCGCCCTTAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGa                              <  1:748721/149‑1 (MQ=255)
                                                                                                                cATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGa                              <  1:202515/149‑1 (MQ=255)
                                                                                                                 aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc                             <  2:1209755/149‑1 (MQ=255)
                                                                                                                  tCCTGCACTACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACg                            <  1:1004232/149‑1 (MQ=255)
                                                                                                                   cgtGCACGACCCACCAGTAAAAAATGTCGCCCTAATGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGc                           <  2:1178117/147‑1 (MQ=255)
                                                                                                                    cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa                          >  2:209575/1‑149 (MQ=255)
                                                                                                                        ccGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACACAGTCAGAGCTTACTGTGCAAGCAACCCTATGTCCGTGGAATTAGGCGTAAAATGCCggacccc                      >  2:40366/2‑143 (MQ=38)
                                                                                                                        aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCt                      <  1:1085687/149‑1 (MQ=255)
                                                                                                                        aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGATTACCTTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCt                      <  1:659528/149‑1 (MQ=255)
                                                                                                                        aCGACCCACCAAGGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCt                      <  1:5844/149‑1 (MQ=255)
                                                                                                                          gACCCACCAATGTAAAAAAACGACCTAAAGGCGCTGTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCCGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGc                    >  2:912603/1‑149 (MQ=255)
                                                                                                                               accaATGTAAAAAAGCGCCCTAAAGGCGCGTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACAtt               <  2:491072/149‑1 (MQ=255)
                                                                                                                               accaATGGAAAAAAGCGCCCTAAAGGCGCTTTTTTTCTAATCAGGCATCCTCAATTTCCCTTTGTAAACCTGAAGACATCGTCAGAGCTTACTGTTCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACAtt               <  1:21717/149‑1 (MQ=255)
                                                                                                                                   aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGc           <  1:271950/149‑1 (MQ=255)
                                                                                                                                      tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCGTCCTGCACATTAGGCGTa        >  1:93544/1‑149 (MQ=255)
                                                                                                                                      tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTa        >  1:93388/1‑149 (MQ=255)
                                                                                                                                        aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAAt      >  2:814983/1‑149 (MQ=255)
                                                                                                                                          aatccacccTGAAGGCGCTTTTTTGCTATTCAGGCATCCTCAACTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACAGTAGGCGTAATTc    <  2:237486/143‑1 (MQ=255)
                                                                                                                                          aaaGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTc    <  1:880263/149‑1 (MQ=255)
                                                                                                                                            aGCGCCCTAAAGGCGCTTTTTTacattggtgggtcgtgtaggatgactcgcttcgctcgcccttcgggccgtcgccgcgagcggctccgttgtttcgctgcgctcaactcaaacctgctgcaggttcgaatcttttagtatcgcagata  >  2:267227‑M1/1‑22 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGA  >  CP009273/777134‑777422

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: