Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,281,042 C→A 25.2% R150S (CGT→AGT)  narJ → molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,281,0420CA25.2% 389.8 / 97.5 170R150S (CGT→AGT) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base C (78/49);  new base A (21/22);  total (99/71)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 172 total aligned reads displayed.

CACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGG  >  CP009273/1280895‑1281189
                                                                                                                                                   |                                                                                                                                                   
caCGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGc                                                                                                                                                     >  1:2852890/1‑148 (MQ=255)
 aCGGCTTGCAGTTTAACAGCCGCGAATTGCCGGACCATCTTCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCg                                                                                                                                                    <  1:2439723/148‑1 (MQ=255)
           gTTAACCAGCCGCGTATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGACGCCGTGGATGGTTTGAAAGATATCGCGCCGATTGTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGt                                                                                                                                          <  2:1208827/148‑1 (MQ=255)
           gTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGt                                                                                                                                          <  2:1208813/148‑1 (MQ=255)
             tAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTa                                                                                                                                        <  2:2988308/148‑1 (MQ=255)
             tAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTa                                                                                                                                        >  1:2210126/1‑148 (MQ=255)
             tAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTAGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTa                                                                                                                                        >  2:550900/1‑148 (MQ=255)
                aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGc                                                                                                                                     <  1:3135803/148‑1 (MQ=255)
                aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGc                                                                                                                                     <  1:907741/148‑1 (MQ=255)
                aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGc                                                                                                                                     <  1:658923/148‑1 (MQ=255)
                aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGc                                                                                                                                     <  1:2090118/148‑1 (MQ=255)
                aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGc                                                                                                                                     <  1:1562363/148‑1 (MQ=255)
                  aGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCg                                                                                                                                   >  2:1931437/1‑148 (MQ=255)
                   gCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGt                                                                                                                                  <  2:706248/148‑1 (MQ=255)
                   gCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGt                                                                                                                                  <  2:706236/148‑1 (MQ=255)
                            tGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGa                                                                                                                         >  2:527609/1‑148 (MQ=255)
                            tGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGa                                                                                                                         >  2:1841435/1‑148 (MQ=255)
                             gCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGAt                                                                                                                        <  2:1945095/148‑1 (MQ=255)
                                 gACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGACGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATctgc                                                                                                                    >  2:3054434/1‑148 (MQ=255)
                                 gACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATctgc                                                                                                                    >  2:2041851/1‑148 (MQ=255)
                                 gACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATctgc                                                                                                                    >  1:433012/1‑148 (MQ=255)
                                   ccATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATctgctg                                                                                                                  <  2:1106995/148‑1 (MQ=255)
                                   ccATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATctgctg                                                                                                                  >  2:3620994/1‑148 (MQ=255)
                                    cATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATctgctgc                                                                                                                 >  1:3750452/1‑148 (MQ=255)
                                     aTCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATctgctgct                                                                                                                <  1:2546345/148‑1 (MQ=255)
                                       cTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGa                                                                                                              >  2:3789270/1‑148 (MQ=255)
                                       cTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGa                                                                                                              >  2:2287712/1‑148 (MQ=255)
                                         gCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTTGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGaaa                                                                                                            <  1:2056164/148‑1 (MQ=255)
                                         gCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGaaa                                                                                                            <  1:1280513/148‑1 (MQ=255)
                                             cTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTgg                                                                                                        >  2:3761051/1‑148 (MQ=255)
                                              tGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGc                                                                                                       <  1:1087769/148‑1 (MQ=255)
                                               gTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCg                                                                                                      >  1:714948/1‑148 (MQ=255)
                                                tATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGa                                                                                                     <  1:3293158/148‑1 (MQ=255)
                                                 aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGac                                                                                                    >  2:3517713/1‑147 (MQ=255)
                                                 aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGaa                                                                                                    >  2:1622304/1‑148 (MQ=255)
                                                 aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGaa                                                                                                    >  1:1871567/1‑148 (MQ=255)
                                                   cTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATa                                                                                                  <  1:3088362/148‑1 (MQ=255)
                                                    tGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGACGTTGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTTTTTGATCTGCTGCTGAAAATGGCGAATAc                                                                                                 <  2:3768208/148‑1 (MQ=255)
                                                      gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCg                                                                                               <  2:1918347/148‑1 (MQ=255)
                                                       aGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGc                                                                                              >  2:1050924/1‑148 (MQ=255)
                                                          aCCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTAt                                                                                           >  2:2935487/1‑148 (MQ=255)
                                                          aCCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTAt                                                                                           <  1:1111574/148‑1 (MQ=255)
                                                           ccTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATc                                                                                          >  2:1645854/1‑148 (MQ=255)
                                                              ggCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGAc                                                                                       <  1:3022797/148‑1 (MQ=255)
                                                                   aGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGa                                                                                  <  1:296041/148‑1 (MQ=255)
                                                                   aGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGa                                                                                  <  1:1931476/148‑1 (MQ=255)
                                                                      tGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACaa                                                                               >  2:495690/1‑148 (MQ=255)
                                                                      tGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACaa                                                                               <  1:1579649/148‑1 (MQ=255)
                                                                      tGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACaa                                                                               >  2:495689/1‑148 (MQ=255)
                                                                           cAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTgg                                                                          >  2:806701/1‑148 (MQ=255)
                                                                             aaGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCg                                                                        <  2:2088549/148‑1 (MQ=255)
                                                                              aGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCgg                                                                       >  2:1675162/1‑148 (MQ=255)
                                                                                cGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaa                                                                     <  2:1158729/148‑1 (MQ=255)
                                                                                   aGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaa                                                                  <  1:1704551/148‑1 (MQ=255)
                                                                                    gacGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATAGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                 >  2:1845150/3‑148 (MQ=255)
                                                                                    gCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                 >  1:531068/1‑148 (MQ=255)
                                                                                    gCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                 >  2:1083713/1‑148 (MQ=255)
                                                                                          gAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCcg                                                           <  1:1910948/148‑1 (MQ=255)
                                                                                          gAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCcg                                                           >  1:2237683/1‑148 (MQ=255)
                                                                                                  gAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAgcg                                                   <  2:1596911/148‑1 (MQ=255)
                                                                                                   aaaGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAgcgc                                                  >  1:460895/1‑148 (MQ=255)
                                                                                                   aaaGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAgcgc                                                  >  1:460926/1‑148 (MQ=255)
                                                                                                    aaGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAcgcgcg                                                 >  2:1379261/1‑148 (MQ=255)
                                                                                                    aaGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAgcgcg                                                 >  2:3653860/1‑148 (MQ=255)
                                                                                                      gATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGcgcgcg                                               <  2:993316/148‑1 (MQ=255)
                                                                                                       atatCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCga                                              <  1:2814501/148‑1 (MQ=255)
                                                                                                         atCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatg                                            >  1:3539986/1‑148 (MQ=255)
                                                                                                         atCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatg                                            >  1:3539643/1‑148 (MQ=255)
                                                                                                         atCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatg                                            >  1:2729820/1‑148 (MQ=255)
                                                                                                           cgcgCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatgat                                          >  1:922953/1‑148 (MQ=255)
                                                                                                                 gATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgccg                                    >  1:914620/1‑148 (MQ=255)
                                                                                                                  aTTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgccgc                                   >  1:1200493/1‑148 (MQ=255)
                                                                                                                   ttCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCa                                  >  2:1353749/1‑148 (MQ=255)
                                                                                                                   ttCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCa                                  >  2:1351586/1‑148 (MQ=255)
                                                                                                                    tCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAg                                 >  1:2472358/1‑148 (MQ=255)
                                                                                                                    tCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAg                                 >  1:2472365/1‑148 (MQ=255)
                                                                                                                       ggCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGgcg                              <  1:1639042/148‑1 (MQ=255)
                                                                                                                        gCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGgcgc                             <  1:2716924/148‑1 (MQ=255)
                                                                                                                         cATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCt                            >  2:1479578/1‑148 (MQ=255)
                                                                                                                          aTTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTg                           >  1:1421625/1‑148 (MQ=255)
                                                                                                                                  gcgcgcGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGt                   <  1:2287712/148‑1 (MQ=255)
                                                                                                                                   cgcgcgTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTc                  <  2:2210126/148‑1 (MQ=255)
                                                                                                                                   cgcgcgTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTc                  <  2:3219197/148‑1 (MQ=255)
                                                                                                                                   cgcgcgTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTc                  <  2:2296597/148‑1 (MQ=255)
                                                                                                                                   cgcgcgACTGCAACAGAGTGAAAGCCTTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCTATGATACGCCGCAGGCGCTGGATGCTGTc                  <  2:2296582/148‑1 (MQ=255)
                                                                                                                                       cgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTggg              >  2:1348463/1‑148 (MQ=255)
                                                                                                                                       cgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTggg              <  2:2351172/148‑1 (MQ=255)
                                                                                                                                          cTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaag           >  2:2707546/1‑148 (MQ=255)
                                                                                                                                          cTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaag           >  1:892449/1‑148 (MQ=255)
                                                                                                                                            gCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaagaa         >  2:379648/1‑148 (MQ=255)
                                                                                                                                            gCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaagaa         >  2:379903/1‑148 (MQ=255)
                                                                                                                                            gCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaagaa         >  2:2933405/1‑148 (MQ=255)
                                                                                                                                                cAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGc     >  1:2571750/1‑148 (MQ=255)
                                                                                                                                                cAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGc     >  1:1749577/1‑148 (MQ=255)
                                                                                                                                                cAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGc     >  1:1748806/1‑148 (MQ=255)
                                                                                                                                                 aGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCa    <  2:3144025/148‑1 (MQ=255)
                                                                                                                                                 aGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCa    <  1:2832801/148‑1 (MQ=255)
                                                                                                                                                  gCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAg   <  1:1554559/148‑1 (MQ=255)
                                                                                                                                                  gCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAg   >  1:71599/1‑148 (MQ=255)
                                                                                                                                                   cGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAgg  >  2:1877908/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGG  >  CP009273/1280895‑1281189

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: