Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,909,857 A→T 6.3% intergenic (‑191/+37) pyrG ← / ← mazG CTP synthetase/nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,909,8570AT 6.3% 323.1 / 7.0 109intergenic (‑191/+37)pyrG/mazGCTP synthetase/nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (54/48);  new base T (3/4);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 110 total aligned reads displayed.

GACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTCGC  >  NC_000913/2909711‑2910003
                                                                                                                                                  |                                                                                                                                                  
  cGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAAt                                                                                                                                               >  1:1719754/1‑150 (MQ=255)
  cGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAAt                                                                                                                                               >  2:348085/1‑150 (MQ=255)
     gAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGac                                                                                                                                            <  1:348085/150‑1 (MQ=255)
     gAAAGCAGTATACCCGAACATGACCTGTGCCACAAAATTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGACATTTAAAAAATGaa                                                                                                                                            >  2:446592/1‑149 (MQ=255)
      aaaGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTTTTAAATGaca                                                                                                                                           <  1:1299489/150‑1 (MQ=255)
      aaaGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGaca                                                                                                                                           >  2:1216925/1‑150 (MQ=255)
       aaGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCACTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGacac                                                                                                                                          >  1:511297/1‑150 (MQ=255)
        atcAGTATAACCGAACATGACCTTTGCCACAAACTTTCTTTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTTTTAAATGacaag                                                                                                                                         <  1:845595/148‑3 (MQ=255)
        aGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTTTTAAATGacaag                                                                                                                                         <  1:606245/150‑3 (MQ=255)
           aGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAAt                                                                                                                                      <  1:728723/150‑1 (MQ=255)
            gTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTTTTAAATGacaagcgct                                                                                                                                     <  2:1644829/150‑7 (MQ=255)
              ataCCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGgc                                                                                                                                   <  1:1562256/150‑1 (MQ=255)
               taCCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGgcg                                                                                                                                  <  1:1480847/150‑1 (MQ=255)
                 cccGAACATTACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCt                                                                                                                                >  1:917672/1‑150 (MQ=255)
                 cccGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCt                                                                                                                                <  1:169706/150‑1 (MQ=255)
                   cGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCATGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTg                                                                                                                              >  1:269794/1‑150 (MQ=255)
                      aCATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGAcc                                                                                                                           >  1:997684/1‑150 (MQ=255)
                       cATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACcg                                                                                                                          <  1:695827/150‑1 (MQ=255)
                       cATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACcg                                                                                                                          >  2:1471035/1‑150 (MQ=255)
                         tGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATg                                                                                                                                              >  1:954354/1‑128 (MQ=255)
                         tGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATg                                                                                                                                              <  2:954354/128‑1 (MQ=255)
                         tGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACcgcg                                                                                                                        >  2:817965/1‑150 (MQ=255)
                         tGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACcgcg                                                                                                                        <  1:1391369/150‑1 (MQ=255)
                          gACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGt                                                                                                                       >  2:135428/1‑150 (MQ=255)
                           aCCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTa                                                                                                                      >  1:1282253/1‑150 (MQ=255)
                            ccTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTaa                                                                                                                     <  2:803240/150‑1 (MQ=255)
                                tgCCACAAACTTTCATTATCCCTCCTCTTCGCCATCGCACTATTGAAATCAGGTGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTcc                                                                                                                 <  1:1125239/150‑1 (MQ=255)
                                 tccACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTa                                                                                                                      >  1:1368087/2‑144 (MQ=255)
                                 gCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTa                                                                                                                      <  2:1368087/144‑1 (MQ=255)
                                 gCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTccc                                                                                                                >  1:1439443/1‑150 (MQ=255)
                                 gCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTccc                                                                                                                >  1:1053725/1‑150 (MQ=255)
                                  ccACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCt                                                                                                               >  2:1374162/1‑150 (MQ=255)
                                        aCTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAgaga                                                                                                         >  1:1237848/1‑150 (MQ=255)
                                        aCTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAgaga                                                                                                         >  1:418825/1‑150 (MQ=255)
                                         cTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATACCCTTAGAGAt                                                                                                        >  2:1077810/1‑150 (MQ=255)
                                             cATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAAt                                                                                                    <  2:958570/150‑1 (MQ=255)
                                              aTTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAAttt                                                                                                  >  2:432427/1‑150 (MQ=255)
                                                tATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTc                                                                                                 <  2:269794/150‑1 (MQ=255)
                                                   cccTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTg                                                                                              >  2:1711346/1‑150 (MQ=255)
                                                           ttCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAACTGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTa                                                                                      <  1:226731/150‑1 (MQ=255)
                                                           ttCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTa                                                                                      <  1:1597475/150‑1 (MQ=255)
                                                             cGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACATATGGCGCTTGACCGCGTAATTCCGTTAGAGATCATTTTCCTGCCGTTTTAcc                                                                                    <  1:109576/150‑1 (MQ=255)
                                                             cGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTAcc                                                                                    <  1:1676532/150‑1 (MQ=255)
                                                               ccAGCGCACTATTTAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTg                                                                                  <  2:1486650/150‑1 (MQ=255)
                                                               ccAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTg                                                                                  >  1:1109276/1‑150 (MQ=255)
                                                                cAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGt                                                                                 <  1:787745/150‑1 (MQ=255)
                                                                  gcgcACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTg                                                                               <  1:1446508/150‑1 (MQ=255)
                                                                       cTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGa                                                                          >  1:732135/1‑150 (MQ=255)
                                                                       cTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGa                                                                          >  2:411240/1‑150 (MQ=255)
                                                                        tATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATTACACAAATGGCGCTTGACCACGAAAATCCCTAAGACATCAATATCCAGCCGTTATACCTATTGCCAGAc                                                                         >  2:310888/1‑150 (MQ=255)
                                                                            gAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGActtct                                                                     >  2:1647217/1‑150 (MQ=255)
                                                                             aaaTCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGActtctt                                                                    <  1:441021/150‑1 (MQ=255)
                                                                              aaTCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACGAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACttcttc                                                                   >  2:36242/1‑150 (MQ=255)
                                                                               aTCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTcc                                                                  >  2:841189/1‑150 (MQ=255)
                                                                               aTCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTcc                                                                  <  1:10751/150‑1 (MQ=255)
                                                                                 cAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACATATGGCGCTTGACCGCGTAATTCGCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCAt                                                                <  1:100849/150‑1 (MQ=255)
                                                                                   ggAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTg                                                              >  2:1451692/1‑150 (MQ=255)
                                                                                    gAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGt                                                             >  1:188543/1‑150 (MQ=255)
                                                                                     aGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGtt                                                            <  1:1451692/150‑1 (MQ=255)
                                                                                      gTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGttt                                                           >  2:169493/1‑150 (MQ=255)
                                                                                       tGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTc                                                          <  1:1216925/150‑1 (MQ=255)
                                                                                          taaaaTAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAg                                                       <  1:1170008/149‑1 (MQ=255)
                                                                                             aaTAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTACCAGACTTCTTCCATTGTTTCGAGGTc                                                    <  2:1109276/150‑1 (MQ=255)
                                                                                               tAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCaa                                                  >  2:265889/1‑150 (MQ=255)
                                                                                                   cccTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACAcc                                              >  2:1140678/1‑150 (MQ=255)
                                                                                                     cTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTg                                            <  1:1711346/150‑1 (MQ=255)
                                                                                                      tttGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGt                                           >  2:1122079/1‑150 (MQ=255)
                                                                                                              aaaTGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCa                                   >  2:978602/1‑150 (MQ=255)
                                                                                                               aaTGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAg                                  >  1:359097/1‑150 (MQ=255)
                                                                                                                     tttGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCAcg                            <  2:66665/150‑1 (MQ=255)
                                                                                                                        gACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCAcgcgc                         <  1:1374162/150‑1 (MQ=255)
                                                                                                                         aCGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCAcgcgcg                        >  2:667021/1‑150 (MQ=255)
                                                                                                                           gCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCTTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGc                      <  2:582080/150‑1 (MQ=255)
                                                                                                                            cAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGca                     <  2:1237848/150‑1 (MQ=255)
                                                                                                                            cAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGca                     <  1:323628/150‑1 (MQ=255)
                                                                                                                            cAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGca                     <  2:418825/150‑1 (MQ=255)
                                                                                                                              aaTCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGcaac                   >  2:1033062/1‑150 (MQ=255)
                                                                                                                              aaTCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTAAGAGATAAATTTCCTGAATTTTAACATTTTGAAAGACTAATTCCATTGTTTCGATGTCAACAAATGTAATTTCCAGTCCAAGCGCGGcaac                   >  2:109129/1‑150 (MQ=255)
                                                                                                                               aTCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGcaaca                  >  1:1030748/1‑150 (MQ=255)
                                                                                                                                 cAAGCGCTTGTCATTTAAAAAATGACACAATTGGCTCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGTAACAAt                <  2:511297/150‑1 (MQ=255)
                                                                                                                                 cAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAAt                <  1:667021/150‑1 (MQ=255)
                                                                                                                                 cAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAAt                >  2:1675834/1‑150 (MQ=255)
                                                                                                                                 cAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACAGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGGCATTTCCAGACCCCGCGCGGCAACAAt                >  2:273375/1‑150 (MQ=255)
                                                                                                                                  aaGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATa               >  2:891287/1‑150 (MQ=255)
                                                                                                                                   aGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAAGAATAc              >  1:1410498/1‑150 (MQ=255)
                                                                                                                                    gctcTTGACATTTAAAAAATGCCACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGCTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACg             <  1:827360/147‑1 (MQ=255)
                                                                                                                                     tttgtGTCATTTTTTAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGc            >  2:446249/5‑150 (MQ=255)
                                                                                                                                     cgcTTGTCATTTAAAAAATGACACGAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGc            <  1:36242/150‑1 (MQ=255)
                                                                                                                                       tgtGTCATTTTTTAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTc          >  2:26427/3‑150 (MQ=255)
                                                                                                                                       tgtGTCATTTTTTAAATGACACAAATGGCGCCTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTTCCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTc          >  1:302731/3‑150 (MQ=255)
                                                                                                                                        ttGTCATTTAAAAAATGACAGATTTGGCGCTTGACCGCGTATTTCCCTTAGAGATAAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACTCCTGTCTTTTCCAGTCCACGCGCTGCAACAATACGCTcc         <  2:128003/150‑1 (MQ=255)
                                                                                                                                        ttGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTcc         <  1:446592/150‑1 (MQ=255)
                                                                                                                                          gTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCAc       <  1:135428/150‑1 (MQ=255)
                                                                                                                                           tCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGCCTGCTTCGATTGTTTCGAGGTCAACACCTGTCTTTTCCAGTCCACGCGCGGCAACAATACGATCCACt      >  2:140303/1‑150 (MQ=255)
                                                                                                                                            cATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACtt     >  1:506780/1‑150 (MQ=255)
                                                                                                                                             aTTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTc    >  1:968011/1‑150 (MQ=255)
                                                                                                                                             aTTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTc    >  1:402243/1‑150 (MQ=255)
                                                                                                                                               ttAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTcgc  <  2:78556/150‑1 (MQ=255)
                                                                                                                                               ttAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTcgc  <  1:169493/150‑1 (MQ=255)
                                                                                                                                               ttAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTcgc  <  2:447743/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GACGGGAAAGCAGTATACCCGAACATGACCTGTGCCACAAACTTTCATTATCCCTCCTCTTCGCCAGCGCACTATTGAAATCAGGAGTGAGAAAATAGCCCCTTTGGGGTAAATGTTTTTGACGCAAATCAAGCGCTTGTCATTTAAAAAATGACACAAATGGCGCTTGACCGCGTAATTCCCTTAGAGATCAATTTCCTGCCGTTTTACCTGTTGCCAGACTTCTTCCATTGTTTCGAGGTCAACACCTGTCATTTCCAGTCCACGCGCGGCAACAATACGCTCCACTTCGC  >  NC_000913/2909711‑2910003

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