Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,489,397 C→A 87.4% 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,3970CA87.4% 801.0 / 114.5 308F676L (TTC→TTAnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (17/22);  new base A (134/134);  total (151/157)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 309 total aligned reads displayed.

ACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCG  >  CP009273/3489248‑3489542
                                                                                                                                                     |                                                                                                                                                 
   ggTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgt                                                                                                                                                >  1:2354944/1‑150 (MQ=255)
     tCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGgtgt                                                                                                                                              <  2:3857534/150‑1 (MQ=255)
      cATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtc                                                                                                                                             >  2:900772/1‑150 (MQ=255)
       aTGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtct                                                                                                                                            >  2:4604341/1‑150 (MQ=255)
       aTGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGtct                                                                                                                                            >  1:4778214/1‑150 (MQ=255)
         gCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTcc                                                                                                                                          <  2:4324488/150‑1 (MQ=255)
         gCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTcc                                                                                                                                          <  2:2781016/150‑1 (MQ=255)
         gCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTcc                                                                                                                                          >  2:2772692/1‑150 (MQ=255)
             aTGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCt                                                                                                                                      >  1:3657957/1‑150 (MQ=255)
             aTGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCt                                                                                                                                      <  1:1566368/150‑1 (MQ=255)
              tGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTg                                                                                                                                     >  2:3558258/1‑150 (MQ=255)
                 gAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTAc                                                                                                                                  >  1:4731967/1‑150 (MQ=255)
                 gAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAc                                                                                                                                  >  1:285293/1‑150 (MQ=255)
                 gAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTAc                                                                                                                                  >  2:2389488/1‑150 (MQ=255)
                  aaaGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTAcc                                                                                                                                 >  1:3562058/1‑150 (MQ=255)
                      cACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTg                                                                                                                             >  1:2335826/1‑150 (MQ=255)
                         tGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAAt                                                                                                                          <  1:637470/150‑1 (MQ=255)
                          gCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATg                                                                                                                         <  1:3814144/150‑1 (MQ=255)
                            gTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGtt                                                                                                                       <  2:3802711/150‑1 (MQ=255)
                               tGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTAGGCGTTGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAg                                                                                                                    >  1:686971/1‑150 (MQ=255)
                                 gCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGaa                                                                                                                  <  2:3529064/150‑1 (MQ=255)
                                   gAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGc                                                                                                                >  2:1443871/1‑150 (MQ=255)
                                    aaaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  2:4716029/150‑1 (MQ=255)
                                    aaaaCCTGCGTGGGGAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCt                                                                                                               <  2:2540583/150‑1 (MQ=255)
                                     aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTc                                                                                                              >  1:1106391/1‑150 (MQ=255)
                                     aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTc                                                                                                              >  1:2004265/1‑150 (MQ=255)
                                     aaaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTc                                                                                                              <  1:2664527/150‑1 (MQ=255)
                                      aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCa                                                                                                             >  1:3476446/1‑150 (MQ=255)
                                      aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCa                                                                                                             <  2:3194748/150‑1 (MQ=255)
                                      aaCCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCa                                                                                                             <  2:674159/150‑1 (MQ=255)
                                        ccTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAgg                                                                                                           >  2:3664450/1‑150 (MQ=255)
                                        ccTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAgg                                                                                                           <  2:627161/150‑1 (MQ=255)
                                          tGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGt                                                                                                         <  1:3410092/150‑1 (MQ=255)
                                          tGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGt                                                                                                         <  1:2401501/150‑1 (MQ=255)
                                           gCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTa                                                                                                        <  1:818530/150‑1 (MQ=255)
                                            cGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaa                                                                                                       <  2:1971604/150‑1 (MQ=255)
                                            cGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTaa                                                                                                       >  1:3245663/1‑150 (MQ=255)
                                                ggTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGAc                                                                                                   <  1:290462/150‑1 (MQ=255)
                                                      aCCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAATGACGTGGGt                                                                                             >  2:2966080/1‑150 (MQ=255)
                                                      aCCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGt                                                                                             >  2:771152/1‑150 (MQ=255)
                                                        cTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtat                                                                                           >  2:48612/1‑150 (MQ=255)
                                                         tGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                          >  2:3928635/1‑150 (MQ=255)
                                                         tGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatt                                                                                          >  1:1492947/1‑150 (MQ=255)
                                                          ggTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGtatta                                                                                         <  2:1195098/150‑1 (MQ=255)
                                                            tGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATc                                                                                       <  1:3143406/150‑1 (MQ=255)
                                                             gCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCg                                                                                      <  2:171785/150‑1 (MQ=255)
                                                                 cTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCAc                                                                                  <  2:1148964/150‑1 (MQ=255)
                                                                  tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACt                                                                                 >  2:1088658/1‑150 (MQ=255)
                                                                  tACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACt                                                                                 >  1:1475584/1‑150 (MQ=255)
                                                                      gCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaa                                                                             >  2:341334/1‑150 (MQ=255)
                                                                        gTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaaaa                                                                           >  1:1536356/1‑150 (MQ=255)
                                                                        gTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGaaaaa                                                                           >  2:4100106/1‑150 (MQ=255)
                                                                         ttGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAg                                                                          <  2:1799454/150‑1 (MQ=255)
                                                                            gCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCt                                                                       <  2:2106010/150‑1 (MQ=255)
                                                                                cAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGaa                                                                   <  2:1366914/150‑1 (MQ=255)
                                                                                cAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGaa                                                                   <  2:3730180/150‑1 (MQ=255)
                                                                                cAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGaa                                                                   <  2:4577224/150‑1 (MQ=255)
                                                                                 aGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAAc                                                                  >  1:3267348/1‑150 (MQ=255)
                                                                                    gTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTg                                                               >  1:2798632/1‑150 (MQ=255)
                                                                                       ggCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTAt                                                            <  2:3030786/150‑1 (MQ=255)
                                                                                         ccTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGt                                                          >  1:167783/1‑150 (MQ=255)
                                                                                           tCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGttt                                                        >  2:3969078/1‑150 (MQ=255)
                                                                                           tCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGttt                                                        <  2:2342891/150‑1 (MQ=255)
                                                                                            cGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCGGTATGTTTg                                                       <  2:249306/150‑1 (MQ=255)
                                                                                              gCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCg                                                     <  1:3545652/150‑1 (MQ=255)
                                                                                                gTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGt                                                   <  1:4505309/150‑1 (MQ=255)
                                                                                                gTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGt                                                   <  2:3426940/150‑1 (MQ=255)
                                                                                                gTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGt                                                   <  2:1025515/150‑1 (MQ=255)
                                                                                                    aaCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcg                                               <  1:755494/150‑1 (MQ=255)
                                                                                                    aaCTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAAcg                                               <  1:1393843/150‑1 (MQ=255)
                                                                                                         gAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACggcggc                                          <  2:3852691/150‑1 (MQ=255)
                                                                                                          aaaaCCGCTACATAGGCATCCGTACGCCGCACAAAATGAAGTTTGGTGTCTCCGGCTGTTCCGGTGTTTGTTCAGTAGCTCAGGGTATAGACGTGGGTTTTATCGACACTGAAAAAGGCTGGATCCTGTATGTTTGCGGTAACGGCGGCa                                         <  1:3978651/150‑1 (MQ=255)
                                                                                                           aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCAt                                        <  2:4505048/150‑1 (MQ=255)
                                                                                                           aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCAt                                        <  1:2236536/150‑1 (MQ=255)
                                                                                                           aaaCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCAt                                        <  1:3746270/150‑1 (MQ=255)
                                                                                                             aCCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGa                                      >  1:1204546/1‑150 (MQ=255)
                                                                                                                   aCAAAGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACcgc                                >  1:169839/1‑150 (MQ=255)
                                                                                                                    cAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTTCGGTAACGGCGGCATGAAACcgcg                               >  2:4217458/1‑150 (MQ=255)
                                                                                                                     aaaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGt                              >  2:3595103/1‑150 (MQ=255)
                                                                                                                     aaaGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGt                              <  2:1312605/150‑1 (MQ=255)
                                                                                                                       aGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCa                            <  2:3320128/150‑1 (MQ=255)
                                                                                                                       aGGCATCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCa                            <  2:2712327/150‑1 (MQ=255)
                                                                                                                           aTCCGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCg                        <  2:2934338/150‑1 (MQ=255)
                                                                                                                             ccGTACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGa                      >  1:4239688/1‑150 (MQ=255)
                                                                                                                                tACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATct                   <  1:949934/150‑1 (MQ=255)
                                                                                                                                tACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATct                   >  1:1855463/1‑150 (MQ=255)
                                                                                                                                tACGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATct                   >  2:2741153/1‑150 (MQ=255)
                                                                                                                                 aCGCCGCACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctg                  <  2:2246395/150‑1 (MQ=255)
                                                                                                                                  cgccgcACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgc                 >  1:586026/1‑150 (MQ=255)
                                                                                                                                   gccgcACAAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATctgct                <  1:900772/150‑1 (MQ=255)
                                                                                                                                       cacaAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcg            >  1:2559400/1‑150 (MQ=255)
                                                                                                                                        acaAAATGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTggcgg           >  1:4567817/1‑150 (MQ=255)
                                                                                                                                          aaaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGgcggcg         >  2:333259/1‑150 (MQ=255)
                                                                                                                                          aaaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGgcggcg         >  1:1017216/1‑150 (MQ=255)
                                                                                                                                           aaaTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGcggcgg        >  2:1548263/1‑150 (MQ=255)
                                                                                                                                             aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  1:2750943/1‑150 (MQ=255)
                                                                                                                                             aTGAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGat      >  1:2673604/1‑150 (MQ=255)
                                                                                                                                               gAAGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGatat    <  2:4201207/150‑1 (MQ=255)
                                                                                                                                                aaGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatc   <  2:1492947/150‑1 (MQ=255)
                                                                                                                                                 aGTTAGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATatcg  >  2:4447691/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCGGTCATGCCTATGCGAAAGCACTGCGTATGGCGAAAACCTGCGTGGGTAGCACCTGGTGCCGCTACGGCGTTGGCGACAGCGTCGGCCTCGGCGTGGAACTGGAAAACCGCTACAAAGGCATCCGTACGCCGCACAAAATGAAGTTCGGTGTCTCCGGCTGTACCCGTGAATGTTCAGAAGCTCAGGGTAAAGACGTGGGTATTATCGCCACTGAAAAAGGCTGGAACCTGTATGTTTGCGGTAACGGCGGCATGAAACCGCGTCATGCGGATCTGCTGGCGGCGGATATCG  >  CP009273/3489248‑3489542

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: