Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,489,397 C→A 30.2% F676L (TTC→TTA nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,489,3970CA30.2% 461.7 / 215.3 247F676L (TTC→TTAnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (79/93);  new base A (36/39);  total (115/132)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

ACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCGATCG  >  CP009273/3489248‑3489546
                                                                                                                                                     |                                                                                                                                                     
aCCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTc                                                                                                                                                       >  2:2465243/1‑150 (MQ=255)
   ggTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGgt                                                                                                                                                    >  2:1268802/1‑150 (MQ=255)
     tCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgtgt                                                                                                                                                  <  1:1450498/150‑1 (MQ=255)
          ccTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCg                                                                                                                                             <  1:4557224/150‑1 (MQ=255)
                 gAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTAc                                                                                                                                      <  1:1047729/150‑1 (MQ=255)
                       aCTGCGTATCGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGa                                                                                                                                <  2:108528/150‑1 (MQ=255)
                        cTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGaa                                                                                                                               <  1:1997110/150‑1 (MQ=255)
                           cGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGt                                                                                                                            >  1:2073255/1‑150 (MQ=255)
                             tATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTc                                                                                                                          >  2:1930328/1‑150 (MQ=255)
                               tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAg                                                                                                                        <  2:3501334/150‑1 (MQ=255)
                               tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAg                                                                                                                        >  1:2366592/1‑150 (MQ=255)
                                    aaatgcTGCGTGGGTAGCTCCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                                   <  2:4438408/145‑1 (MQ=255)
                                    aaaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                                   >  1:1702682/1‑150 (MQ=255)
                                      taCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCa                                                                                                                 <  2:4797265/149‑1 (MQ=255)
                                      aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCa                                                                                                                 <  1:133096/150‑1 (MQ=255)
                                       aCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAg                                                                                                                <  2:4180275/150‑1 (MQ=255)
                                         cTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAggg                                                                                                              <  1:150422/150‑1 (MQ=255)
                                             gTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                          >  1:4314618/1‑150 (MQ=255)
                                             gTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                          <  1:4252596/150‑1 (MQ=255)
                                             gTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                          >  1:1774378/1‑150 (MQ=255)
                                             gTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaaa                                                                                                          <  2:2627447/150‑1 (MQ=255)
                                               gggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGa                                                                                                        <  2:1245415/150‑1 (MQ=255)
                                               gggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGa                                                                                                        >  1:4843636/1‑150 (MQ=255)
                                               gggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCGACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGa                                                                                                        <  2:1469629/150‑1 (MQ=255)
                                                ggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGAc                                                                                                       <  2:753986/150‑1 (MQ=255)
                                                ggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGAc                                                                                                       <  1:3792959/150‑1 (MQ=255)
                                                 gTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACg                                                                                                      <  2:2186237/150‑1 (MQ=255)
                                                  tAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGt                                                                                                     >  1:1246354/1‑150 (MQ=255)
                                                    gCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTgg                                                                                                   >  2:3789523/1‑150 (MQ=255)
                                                     cACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTggg                                                                                                  <  1:2579103/150‑1 (MQ=255)
                                                      aCCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGt                                                                                                 <  2:3622617/150‑1 (MQ=255)
                                                       ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGta                                                                                                >  2:3908525/1‑150 (MQ=255)
                                                       ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGta                                                                                                >  1:1870486/1‑150 (MQ=255)
                                                       ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGta                                                                                                <  2:4641078/150‑1 (MQ=255)
                                                       ccTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGta                                                                                                <  2:770137/150‑1 (MQ=255)
                                                        cTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtat                                                                                               <  1:3706248/150‑1 (MQ=255)
                                                        cTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtat                                                                                               <  2:2421953/150‑1 (MQ=255)
                                                         tGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                              <  2:2636950/150‑1 (MQ=255)
                                                         tGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                              >  1:1280795/1‑150 (MQ=255)
                                                          ggTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatta                                                                                             <  1:3855829/150‑1 (MQ=255)
                                                           gTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGACGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtattat                                                                                            >  2:1658199/1‑150 (MQ=255)
                                                               cGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGcc                                                                                        >  2:1790472/1‑150 (MQ=255)
                                                                  tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACt                                                                                     >  1:607443/1‑150 (MQ=255)
                                                                  tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACt                                                                                     >  1:3983907/1‑150 (MQ=255)
                                                                  tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACt                                                                                     >  1:3984094/1‑150 (MQ=255)
                                                                    cGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGa                                                                                   <  2:4867256/150‑1 (MQ=255)
                                                                      gCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCACCACTGaaa                                                                                 >  1:2635875/1‑150 (MQ=255)
                                                                        gTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaaaa                                                                               <  1:4224612/150‑1 (MQ=255)
                                                                         ttGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAg                                                                              <  2:1870486/150‑1 (MQ=255)
                                                                           ggCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGc                                                                            <  1:4944249/150‑1 (MQ=255)
                                                                            gCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCt                                                                           <  2:3708082/150‑1 (MQ=255)
                                                                               aCAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGa                                                                        >  2:4794233/1‑150 (MQ=255)
                                                                                 aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                      >  2:5185286/1‑150 (MQ=255)
                                                                                 aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                      >  1:2572580/1‑150 (MQ=255)
                                                                                 aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                      <  2:1993609/150‑1 (MQ=255)
                                                                                 aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                      <  2:2865151/150‑1 (MQ=255)
                                                                                 aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                      >  2:4065251/1‑150 (MQ=255)
                                                                                  gCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAcc                                                                     <  2:197639/150‑1 (MQ=255)
                                                                                   cGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCt                                                                    >  2:3017327/1‑150 (MQ=255)
                                                                                     tCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGt                                                                  <  1:4057556/150‑1 (MQ=255)
                                                                                      cGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTa                                                                 <  1:711190/150‑1 (MQ=255)
                                                                                      cGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTa                                                                 >  1:3573116/1‑150 (MQ=255)
                                                                                       ggCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTAt                                                                >  1:4950067/1‑150 (MQ=255)
                                                                                        gCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATg                                                               >  1:70393/1‑150 (MQ=255)
                                                                                          cTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGtt                                                             >  2:4458582/1‑150 (MQ=255)
                                                                                              gCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCg                                                         <  2:295800/150‑1 (MQ=255)
                                                                                              gCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCg                                                         <  2:4960451/150‑1 (MQ=255)
                                                                                               cGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCgg                                                        >  1:377558/1‑150 (MQ=255)
                                                                                                gTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGt                                                       >  2:2761668/1‑150 (MQ=255)
                                                                                                     aCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcgg                                                  <  1:4432040/150‑1 (MQ=255)
                                                                                                      ttGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcggc                                                 >  2:601851/2‑150 (MQ=255)
                                                                                                       tGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTa                                                      <  1:1588883/144‑1 (MQ=255)
                                                                                                       tGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTa                                                      >  2:1588869/1‑144 (MQ=255)
                                                                                                       tGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcggcg                                                >  2:812629/1‑150 (MQ=255)
                                                                                                        ggAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcggcgg                                               >  1:2352602/1‑150 (MQ=255)
                                                                                                        ggAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcggcgg                                               >  1:2885473/1‑150 (MQ=255)
                                                                                                           aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCAt                                            <  1:333751/150‑1 (MQ=255)
                                                                                                              ccGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaa                                         >  1:4682552/1‑150 (MQ=255)
                                                                                                               cGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGaaa                                        <  2:1774378/150‑1 (MQ=255)
                                                                                                                  tACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcg                                     >  2:230975/1‑150 (MQ=255)
                                                                                                                   aCAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgc                                    >  2:4270600/1‑150 (MQ=255)
                                                                                                                      aaGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTc                                 <  2:4576564/150‑1 (MQ=255)
                                                                                                                      aaGGCATCCGGACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTc                                 <  2:3560064/150‑1 (MQ=255)
                                                                                                                             ccGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGa                          >  1:2734485/1‑150 (MQ=255)
                                                                                                                               gTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATc                        >  1:3149753/1‑150 (MQ=255)
                                                                                                                                    ccgcACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgctg                   <  2:2561664/150‑1 (MQ=255)
                                                                                                                                    ccgcACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgctg                   <  2:2561834/150‑1 (MQ=255)
                                                                                                                                     cgcACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgctcg                  >  2:4625232/1‑148 (MQ=255)
                                                                                                                                      gcACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggc                 >  2:3939678/1‑150 (MQ=255)
                                                                                                                                      gcACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggc                 >  2:2479559/1‑150 (MQ=255)
                                                                                                                                      gcACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggc                 <  2:2507398/150‑1 (MQ=255)
                                                                                                                                        acaAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcgg               <  2:3137166/150‑1 (MQ=255)
                                                                                                                                         caAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcggc              <  2:1280795/150‑1 (MQ=255)
                                                                                                                                          aaaaTGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGgcggcg             <  2:23308/150‑1 (MQ=255)
                                                                                                                                            aaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGa           >  2:3115506/1‑150 (MQ=255)
                                                                                                                                              tGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGata         <  2:2243692/150‑1 (MQ=255)
                                                                                                                                                  gTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcga     <  1:3030827/150‑1 (MQ=255)
                                                                                                                                                  gTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcga     >  1:3278834/1‑150 (MQ=255)
                                                                                                                                                     cGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcgatcg  <  1:4042405/150‑1 (MQ=255)
                                                                                                                                                     cGGTGTCTCCGGCTGGACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcgatcg  <  2:3517871/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCGATCG  >  CP009273/3489248‑3489546

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: