Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,902,059 C→T 6.6% V323I (GTC→ATC)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,0590CT 6.6% 996.8 / 40.1 287V323I (GTC→ATC) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (156/112);  new base T (7/12);  total (163/124)
Fisher's exact test for biased strand distribution p-value = 9.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

GTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGAC  >  CP009273/2901913‑2902207
                                                                                                                                                  |                                                                                                                                                    
gTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCt                                                                                                                                                   >  1:1759453/1‑150 (MQ=255)
    tcACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGaca                                                                                                                                               >  2:565620/1‑150 (MQ=255)
      aCGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACg                                                                                                                                             >  1:3530755/1‑150 (MQ=255)
          gcAAAACGCGCGGTCGTAATCATTCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttc                                                                                                                                         <  1:3856867/150‑1 (MQ=255)
          gcAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttc                                                                                                                                         >  2:1718865/1‑150 (MQ=255)
           cAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttct                                                                                                                                        >  2:2063564/1‑150 (MQ=255)
             aaaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttcttc                                                                                                                                      <  1:3435970/150‑1 (MQ=255)
             aaaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttcttc                                                                                                                                      >  1:1480688/1‑150 (MQ=255)
              aaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCa                                                                                                                                     >  2:2339068/1‑150 (MQ=255)
               aCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAg                                                                                                                                    <  1:155414/150‑1 (MQ=255)
               aCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAg                                                                                                                                    <  2:3809483/150‑1 (MQ=255)
                   gcgGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCa                                                                                                                                >  1:2184725/1‑150 (MQ=255)
                    cgGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCAc                                                                                                                               >  2:716917/1‑150 (MQ=255)
                     ggTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTTATGTTGACGCTGACACGATTCTTCAGCCCAcc                                                                                                                              >  1:4167424/1‑150 (MQ=255)
                     ggTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTTATGTTGACGCTGACACGATTCTTCAGCCCAcc                                                                                                                              >  1:4097949/1‑150 (MQ=255)
                     ggTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCAcc                                                                                                                              >  2:448816/1‑150 (MQ=255)
                     ggTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCAcc                                                                                                                              >  2:446716/1‑150 (MQ=255)
                        cGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCgtg                                                                                                                           >  2:4045775/1‑150 (MQ=255)
                        cGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCgtg                                                                                                                           >  1:2379270/1‑150 (MQ=255)
                        cGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCgtg                                                                                                                           >  1:2379161/1‑150 (MQ=255)
                          tAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGtt                                                                                                                         >  2:1101816/1‑150 (MQ=255)
                          tAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGtt                                                                                                                         <  2:1803490/150‑1 (MQ=255)
                          tAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGtt                                                                                                                         <  2:1802223/150‑1 (MQ=255)
                          tAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGtt                                                                                                                         >  2:2009019/1‑150 (MQ=255)
                          tAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGtt                                                                                                                         >  2:2008829/1‑150 (MQ=255)
                            aTCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTc                                                                                                                       >  1:770583/1‑150 (MQ=255)
                             tCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCa                                                                                                                      >  1:2102245/1‑150 (MQ=255)
                               aTGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGt                                                                                                                    <  2:2601974/150‑1 (MQ=255)
                               aTGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGt                                                                                                                    >  2:3637513/1‑150 (MQ=255)
                                tGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTGAAACAACTTGTGAATCGATCAGTATGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGATTCAGTg                                                                                                                   >  2:1373004/1‑150 (MQ=255)
                                  ccTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCt                                                                                                                 >  2:1111067/1‑150 (MQ=255)
                                    ttCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTc                                                                                                               <  2:1046572/150‑1 (MQ=255)
                                     tCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCg                                                                                                              <  2:29222/150‑1 (MQ=255)
                                     tCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCg                                                                                                              <  2:3455574/150‑1 (MQ=255)
                                        aCGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATc                                                                                                           >  2:2103134/1‑150 (MQ=255)
                                          gCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCAc                                                                                                         >  1:3658955/1‑150 (MQ=255)
                                           ccACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACt                                                                                                        <  1:1151185/150‑1 (MQ=255)
                                            cACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTg                                                                                                       <  1:3923246/150‑1 (MQ=255)
                                            cACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTg                                                                                                       <  1:1666131/150‑1 (MQ=255)
                                                  aGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTtat                                                                                                 >  1:3905143/1‑150 (MQ=255)
                                                  aGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTtat                                                                                                 >  1:3905087/1‑150 (MQ=255)
                                                     cGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAg                                                                                              >  2:2958335/1‑150 (MQ=255)
                                                     cGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAg                                                                                              >  2:3796966/1‑150 (MQ=255)
                                                     cGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAg                                                                                              >  2:1536029/1‑150 (MQ=255)
                                                          ccGCCAGGTTCGATGATTGCGTCCAGAACATTAATGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATTGAACAGTTTGATGTTGTCGCTGACCCGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGTTCACTGTTTTATAAGCATcc                                                                                         <  2:586971/150‑1 (MQ=255)
                                                          ccGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATcc                                                                                         <  1:4280230/150‑1 (MQ=255)
                                                          ccGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATcc                                                                                         <  2:4508074/150‑1 (MQ=255)
                                                           cGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACACCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCg                                                                                        <  1:4383619/150‑1 (MQ=255)
                                                             ccAGGTACGAGGAGTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGc                                                                                      <  2:2505349/150‑1 (MQ=255)
                                                              cAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCa                                                                                     >  1:132715/1‑150 (MQ=255)
                                                              cAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCa                                                                                     >  1:132281/1‑150 (MQ=255)
                                                                 gTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGTATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGtt                                                                                  <  2:855198/150‑1 (MQ=255)
                                                                   aCGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATACGGCAGTTCa                                                                                >  2:2299712/1‑150 (MQ=255)
                                                                    cGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCaa                                                                               >  1:4321901/1‑150 (MQ=255)
                                                                     gAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAAt                                                                              >  2:1713190/1‑150 (MQ=255)
                                                                           tGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACtt                                                                        <  1:3803400/150‑1 (MQ=255)
                                                                               tCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCg                                                                    >  2:1289363/1‑150 (MQ=255)
                                                                                   gACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCa                                                                <  1:624719/150‑1 (MQ=255)
                                                                                       tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                            >  2:1533544/1‑150 (MQ=255)
                                                                                       tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                            >  1:4470104/1‑150 (MQ=255)
                                                                                       tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                            >  1:4470242/1‑150 (MQ=255)
                                                                                       tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                            >  2:1983026/1‑150 (MQ=255)
                                                                                         tAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGACGACCATACCGa                                                          >  2:4400033/1‑150 (MQ=255)
                                                                                           aGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATg                                                        <  1:2741418/150‑1 (MQ=255)
                                                                                              aTTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTg                                                     >  2:3681308/1‑150 (MQ=255)
                                                                                                 tCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACt                                                  >  1:3277062/1‑150 (MQ=255)
                                                                                                  cAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACtt                                                 <  2:3360175/150‑1 (MQ=255)
                                                                                                   aaCGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTc                                                <  2:2184155/150‑1 (MQ=255)
                                                                                                   aaCGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTc                                                >  2:3308982/1‑150 (MQ=255)
                                                                                                   aaCGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTc                                                >  2:247681/1‑150 (MQ=255)
                                                                                                     cgccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCAc                                              >  2:3432232/1‑150 (MQ=255)
                                                                                                      gccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACt                                             <  1:2819993/150‑1 (MQ=255)
                                                                                                           gcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCg                                        <  1:929402/150‑1 (MQ=255)
                                                                                                             gTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCggg                                      <  1:4304462/150‑1 (MQ=255)
                                                                                                             gTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCggg                                      >  2:524050/1‑150 (MQ=255)
                                                                                                               ttCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGtt                                    >  2:2844317/1‑150 (MQ=255)
                                                                                                                tCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                   <  2:367161/150‑1 (MQ=255)
                                                                                                                 cAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCg                                  <  2:3887458/150‑1 (MQ=255)
                                                                                                                  aaCATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGc                                 <  1:1210554/150‑1 (MQ=255)
                                                                                                                    cATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGctt                               >  1:1559477/1‑150 (MQ=255)
                                                                                                                     aTCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttc                              <  1:459350/150‑1 (MQ=255)
                                                                                                                       cTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                            >  2:1231797/1‑150 (MQ=255)
                                                                                                                       cTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                            <  2:2108545/150‑1 (MQ=255)
                                                                                                                       cTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                            <  2:1898114/150‑1 (MQ=255)
                                                                                                                         tgtgAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCg                          <  1:2182264/150‑1 (MQ=255)
                                                                                                                         tgtgAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCg                          <  2:2121623/150‑1 (MQ=255)
                                                                                                                          gtgAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGa                         >  1:3615110/1‑150 (MQ=255)
                                                                                                                           tgAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGaa                        >  1:3132433/1‑150 (MQ=255)
                                                                                                                           tgAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGaa                        >  2:2779919/1‑150 (MQ=255)
                                                                                                                            gAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAg                       >  2:1126825/1‑150 (MQ=255)
                                                                                                                             aaTCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGa                      >  2:3088276/1‑150 (MQ=255)
                                                                                                                             aaTCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGa                      >  2:3011205/1‑150 (MQ=255)
                                                                                                                              aTCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGAt                     >  2:1419552/1‑150 (MQ=255)
                                                                                                                               tCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATa                    >  1:1002551/1‑150 (MQ=255)
                                                                                                                                  aTCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAAcc                 >  1:2633904/1‑150 (MQ=255)
                                                                                                                                  aTCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAAcc                 <  1:2070053/150‑1 (MQ=255)
                                                                                                                                      gTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGtt             >  1:3504360/1‑150 (MQ=255)
                                                                                                                                         tGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTccc          <  1:1841684/150‑1 (MQ=255)
                                                                                                                                               tGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCgg    >  2:4186061/1‑150 (MQ=255)
                                                                                                                                                 aCGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGAc  <  1:568303/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGAC  >  CP009273/2901913‑2902207

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: