Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,274,979 G→T 14.9% intergenic (+175/‑341) narK → / → narG nitrate/nitrite transporter/nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,274,9790GT14.9% 668.7 / 67.7 221intergenic (+175/‑341)narK/narGnitrate/nitrite transporter/nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (105/83);  new base T (18/15);  total (123/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 222 total aligned reads displayed.

TGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAGGTATAAAAC  >  CP009273/1274833‑1275125
                                                                                                                                                  |                                                                                                                                                  
tgtgCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATga                                                                                                                                                   >  2:657992/1‑147 (MQ=255)
  tgctgcGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGt                                                                                                                                                 <  1:3143218/148‑1 (MQ=255)
    ctgcggctgcgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCa                                                                                                                                               <  2:1150980/148‑1 (MQ=255)
       cggctgcgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATATAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTAGTCAttt                                                                                                                                            >  1:1945919/1‑148 (MQ=255)
             cgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTac                                                                                                                                      >  2:1870906/1‑148 (MQ=255)
               ttcattcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTacaa                                                                                                                                    <  2:2035123/148‑1 (MQ=255)
                tcattcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTacaac                                                                                                                                   >  1:2465672/1‑148 (MQ=255)
                tcattcaCCCCAGTCACTTACTTTAGTAAGATCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTacaac                                                                                                                                   >  2:1592340/1‑148 (MQ=255)
                 cattcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTacaaca                                                                                                                                  <  1:2947890/148‑1 (MQ=255)
                 cattcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTacaaca                                                                                                                                  <  2:919406/148‑1 (MQ=255)
                 cagtcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTacaaca                                                                                                                                  <  2:1934656/145‑1 (MQ=255)
                  attcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACAt                                                                                                                                 >  2:137276/1‑148 (MQ=255)
                   ttcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATa                                                                                                                                <  1:66145/148‑1 (MQ=255)
                      aCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTa                                                                                                                             <  1:1108123/148‑1 (MQ=255)
                      aCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTa                                                                                                                             <  2:384890/148‑1 (MQ=255)
                        cccAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAAACGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCATTTTTTATGCGTCATTTAGTTACAACATACTAAt                                                                                                                           >  2:1937587/1‑148 (MQ=255)
                           aGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGtt                                                                                                                        <  2:2722511/148‑1 (MQ=255)
                                  tacttTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGt                                                                                                                 <  2:1086315/148‑1 (MQ=255)
                                   acttTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGtt                                                                                                                <  1:574829/148‑1 (MQ=255)
                                    cttTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGttt                                                                                                               <  2:2423255/148‑1 (MQ=255)
                                      ttAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTAt                                                                                                             <  2:1950461/148‑1 (MQ=255)
                                       tAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTAtt                                                                                                            >  1:416569/1‑148 (MQ=255)
                                        aGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTAttt                                                                                                           >  2:1741697/1‑148 (MQ=255)
                                         gTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTc                                                                                                          >  2:1766016/1‑148 (MQ=255)
                                          tAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCg                                                                                                         <  1:90154/148‑1 (MQ=255)
                                           aaGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGc                                                                                                        >  1:3257711/1‑148 (MQ=255)
                                             gCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCg                                                                                                      >  1:3040895/1‑148 (MQ=255)
                                             gCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCg                                                                                                      >  1:515258/1‑148 (MQ=255)
                                               tCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGa                                                                                                    <  2:702081/148‑1 (MQ=255)
                                                ccTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGAt                                                                                                   <  1:2814006/148‑1 (MQ=255)
                                                  tGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGAttt                                                                                                 >  1:1146255/1‑148 (MQ=255)
                                                   gggATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTc                                                                                                <  1:7212/148‑1 (MQ=255)
                                                    ggATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCa                                                                                               >  2:1116534/1‑148 (MQ=255)
                                                    ggATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCa                                                                                               >  1:1969759/1‑148 (MQ=255)
                                                     gATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCAt                                                                                              >  1:3009213/1‑148 (MQ=255)
                                                      attcattcACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCAtt                                                                                             >  1:2285989/1‑148 (MQ=255)
                                                           ttcaCTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAagag                                                                                        >  2:2379118/1‑148 (MQ=255)
                                                           gtcaCTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAagag                                                                                        <  2:1449382/147‑1 (MQ=255)
                                                                   ccgccTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAat                                                                                <  1:1517265/148‑1 (MQ=255)
                                                                   ccgccTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAat                                                                                <  2:2116481/148‑1 (MQ=255)
                                                                   ccgccTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAat                                                                                <  1:3031817/148‑1 (MQ=255)
                                                                      ccTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATg                                                                             <  2:1360022/148‑1 (MQ=255)
                                                                        ttCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGtt                                                                           >  2:979754/1‑148 (MQ=255)
                                                                        ttCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATCTGtt                                                                           <  2:2532055/148‑1 (MQ=255)
                                                                         tCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTa                                                                          <  1:1245346/148‑1 (MQ=255)
                                                                            tGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTAccc                                                                       >  2:545802/1‑148 (MQ=255)
                                                                             gTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCa                                                                      >  1:2714977/1‑148 (MQ=255)
                                                                                aaCCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATgg                                                                   >  2:1082844/1‑148 (MQ=255)
                                                                                 aCCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATggg                                                                  >  1:1994596/1‑148 (MQ=255)
                                                                                 aCCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATggg                                                                  >  2:513628/1‑148 (MQ=255)
                                                                                  ccGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATgggg                                                                 <  1:541168/148‑1 (MQ=255)
                                                                                   cGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGa                                                                <  1:2865232/148‑1 (MQ=255)
                                                                                      aTTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATa                                                             >  1:1326060/1‑148 (MQ=255)
                                                                                      aTTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATa                                                             >  2:2062575/1‑148 (MQ=255)
                                                                                       ttATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATAc                                                            >  2:1673213/1‑148 (MQ=255)
                                                                                         atataGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTc                                                          <  2:3295053/148‑1 (MQ=255)
                                                                                          tataGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACACATGGGGAATACTca                                                         >  1:2235747/1‑147 (MQ=255)
                                                                                           ataGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCt                                                        >  1:2170411/1‑148 (MQ=255)
                                                                                           ataGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCt                                                        <  1:101689/148‑1 (MQ=255)
                                                                                           ataGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCt                                                        <  1:1054101/148‑1 (MQ=255)
                                                                                            taGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCtt                                                       >  2:2759183/1‑148 (MQ=255)
                                                                                             agagTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTa                                                      >  1:930277/1‑148 (MQ=255)
                                                                                              gagTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTaa                                                     >  1:2693939/1‑148 (MQ=255)
                                                                                               agTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAAt                                                    >  1:716347/1‑148 (MQ=255)
                                                                                                gTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATa                                                   >  1:1733877/1‑148 (MQ=255)
                                                                                                 tAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATAc                                                  >  1:2612493/1‑148 (MQ=255)
                                                                                                 tAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATAc                                                  <  1:3191845/148‑1 (MQ=255)
                                                                                                       atTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCt                                            >  1:3062210/1‑148 (MQ=255)
                                                                                                        tttGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTg                                           >  1:2879437/1‑148 (MQ=255)
                                                                                                         ttGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGc                                          <  2:2465672/148‑1 (MQ=255)
                                                                                                          tGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCa                                         <  2:229132/148‑1 (MQ=255)
                                                                                                               aTCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATaaaa                                    >  2:1529347/1‑148 (MQ=255)
                                                                                                                tCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATaaaaa                                   <  1:1150068/148‑1 (MQ=255)
                                                                                                                tCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATaaaaa                                   >  2:2966459/1‑148 (MQ=255)
                                                                                                                tCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATaaaaa                                   <  1:2800691/148‑1 (MQ=255)
                                                                                                                 ccTTTGCGCGGCATGATGTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAAt                                  >  2:3083113/1‑148 (MQ=255)
                                                                                                                  cTTTGCGCGGCTTGATGTCGCGATTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGTTTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATc                                 <  2:213202/148‑1 (MQ=255)
                                                                                                                  cTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATc                                 <  2:3252175/148‑1 (MQ=255)
                                                                                                                     tGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTa                              <  1:1735406/148‑1 (MQ=255)
                                                                                                                      gcgcGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTaa                             >  1:2774253/1‑148 (MQ=255)
                                                                                                                       cgcgGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAAt                            >  1:3296578/1‑148 (MQ=255)
                                                                                                                       cgcgGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAAt                            <  1:1445178/148‑1 (MQ=255)
                                                                                                                         cgGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAg                          >  2:2907818/1‑148 (MQ=255)
                                                                                                                          ggCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGt                         <  2:1969832/148‑1 (MQ=255)
                                                                                                                           gCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGtt                        >  2:2271036/1‑148 (MQ=255)
                                                                                                                              tgatgTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTaa                     <  1:351838/148‑1 (MQ=255)
                                                                                                                                atgTCGCGCTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAAt                   >  2:35342/1‑148 (MQ=255)
                                                                                                                                atgTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAAt                   >  2:1261223/1‑148 (MQ=255)
                                                                                                                                atgTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAAt                   >  2:3135057/1‑148 (MQ=255)
                                                                                                                                   tCGCGCTTTTTTTATGCGTCATTTAGTTACAACATATTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAAc                <  1:2733844/148‑1 (MQ=255)
                                                                                                                                     gcgcTTTTTTTATTCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTa              <  1:979754/148‑1 (MQ=255)
                                                                                                                                     gcgcTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTa              <  2:797040/148‑1 (MQ=255)
                                                                                                                                      cgcTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTAc             <  1:621562/148‑1 (MQ=255)
                                                                                                                                        cTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAg           <  1:2595881/148‑1 (MQ=255)
                                                                                                                                         tttttttATTCTTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAgg          <  2:1945919/148‑1 (MQ=255)
                                                                                                                                         tttttttATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAgg          <  2:716414/148‑1 (MQ=255)
                                                                                                                                          ttttttATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAGGt         >  1:1794298/1‑148 (MQ=255)
                                                                                                                                          ttttttATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAGGt         >  1:1415606/1‑148 (MQ=255)
                                                                                                                                             tttATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAGGtata      >  2:3059018/1‑148 (MQ=255)
                                                                                                                                                 tGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAGGTATAAAAc  >  2:521519/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTATTTCGCCGGATTTCATTAAGAGCCATTAATATGTTACCCATGGGGAATACTCCTTAATACCCATCTGCATAAAAATCTTAATAGTTTAAATAACTACAGGTATAAAAC  >  CP009273/1274833‑1275125

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