Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,248 A→T 6.7% intergenic (+367/‑211) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,2480AT 6.7% 208.4 / 5.0 84intergenic (+367/‑211)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base T (4/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCT  >  NC_000913/375104‑375397
                                                                                                                                                |                                                                                                                                                     
gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGACTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAGCTGCGCACGACTGCCGGATGCGGCGTAAACGcc                                                                                                                                                  <  1:282978/150‑1 (MQ=255)
    cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAACCGCCTTAt                                                                                                                                              >  1:167139/1‑150 (MQ=37)
    cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATATGCGCACGACTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                              >  2:749784/1‑150 (MQ=37)
      cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCACATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                            <  2:858726/150‑1 (MQ=37)
      aaTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCACATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                            <  2:1338106/149‑1 (MQ=37)
       aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGACTACGGGTGGCGCGAGAATTTGTAGGCCAGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTTCCGGATGCGGCGTACACCCCTTATCCg                                                                                                                                           >  2:260942/1‑150 (MQ=11)
        tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                          <  1:1482600/150‑1 (MQ=17)
         cAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                         <  2:1387010/150‑1 (MQ=37)
         cAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCATTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTAGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                         >  2:222692/1‑150 (MQ=255)
          cggCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                        <  1:179260/149‑1 (MQ=37)
          aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATAAGGCGTTTACGCCGCATCCGGca                                                                                                                                        <  2:204731/150‑2 (MQ=11)
          aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGAACGACTGCCGGATGCGGAGTAAACGCCTTTTGCGGcc                                                                                                                                        >  1:833355/1‑150 (MQ=255)
           ggCATCTGCGCACGACTGCGGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                       <  1:532661/150‑1 (MQ=37)
                 tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGTTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGGATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                                 <  2:1100851/150‑1 (MQ=25)
                 tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                                 >  1:594897/1‑150 (MQ=38)
                     cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                             >  2:1005054/1‑150 (MQ=39)
                        gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCgg                                                                                                                          >  2:1030473/1‑150 (MQ=39)
                          cggCCGGATGCGGCGTGAACGCCTTCTCCGGCCTACGGGTGGCGCTAGAATTTGTAGGCCTGATACGACGCGTTAGCGTCGCATCAGGCAGCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                        <  1:813111/148‑1 (MQ=255)
                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                        >  2:750156/3‑150 (MQ=14)
                               ggataaggcgTTCACGCCGCATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                                                   >  2:1465154/7‑150 (MQ=14)
                                gATGCGGCGTGAACGCCTTTTCCTTCCTTCGTTTGGCGCTAGAAGTTGTTGGGCTGATAAGACGCCGTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                                  <  2:265275/150‑1 (MQ=11)
                                     ggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                             <  1:1419830/150‑1 (MQ=39)
                                     ggcgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                             >  1:1315325/1‑150 (MQ=39)
                                        gTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCCTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGa                                                                                                          >  1:443065/1‑150 (MQ=255)
                                         tcaaCTCCTAATCCGGCCTACCGGCGGCGCGAGAATTTGTACGCCGGATAACACGCCGTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTAACCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                         <  1:1170710/148‑1 (MQ=14)
                                          taaCGCCTTATCCGGCCTACGGTTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                        <  2:1002179/149‑1 (MQ=37)
                                            acgggtgaTCGGGCCTACGGGTGGCGCGAGTCTTAGTTGCCCTGTTAAGACGCGTTTGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAg                                                                                                      <  1:300675/143‑1 (MQ=14)
                                             cGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                                                                     >  2:1276131/1‑150 (MQ=39)
                                               ccTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                   >  1:579276/1‑150 (MQ=39)
                                               ccTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                   >  1:1724873/1‑150 (MQ=39)
                                               ccTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                   >  2:220310/1‑150 (MQ=39)
                                               ccTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAAGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                   >  1:364369/1‑150 (MQ=255)
                                                cTTATCGGGCCTACGGGGGGCGAGAGAATTTGTAGGGCTGATACGACGCGGTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                  <  2:1091013/150‑1 (MQ=11)
                                                   aTCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGtt                                                                                               >  2:1643199/1‑150 (MQ=39)
                                                    tCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                                                              >  1:1439004/1‑150 (MQ=39)
                                                            aCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                                                      >  2:579978/1‑150 (MQ=39)
                                                                ttGGCGCGAGAATTTGTAGGCCTGATCAGACGCGTTAGCGTCGCATCTGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                                  <  2:1335509/149‑1 (MQ=255)
                                                                gTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTAATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                                  <  1:908548/150‑1 (MQ=38)
                                                                gTGGCGCGAGAATTGGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTCATCCGGCCGACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCTTCAg                                                                                  <  1:908546/150‑1 (MQ=255)
                                                                     gcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                                             >  2:412730/1‑150 (MQ=35)
                                                                     gcgAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                                             >  2:1477153/1‑150 (MQ=35)
                                                                       gagaATTTGTAGGCCTGATACGACGCGTAAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTg                                                                           <  1:1116623/150‑1 (MQ=12)
                                                                               gTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGAc                                                                   >  2:1037450/1‑150 (MQ=37)
                                                                                  ggCCTGATAAGACGCGTTAGCGTCGCATCCGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGc                                                                <  1:370215/150‑1 (MQ=37)
                                                                                  ggCCTGATAAGACGCGTTAGCGTCGCATCCGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAAATTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCAGCACACGACGGc                                                                >  1:597747/1‑150 (MQ=255)
                                                                                   gCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCgcacaccctgc                                                               >  2:1592/1‑143 (MQ=9)
                                                                                   gCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCCCACGCCTGcc                                                               >  1:534802/1‑150 (MQ=18)
                                                                                     cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATATGCACACGACTGCCgg                                                             >  2:1556270/1‑150 (MQ=21)
                                                                                       gttatgACGCGTTTTCGTCGAATCATGAATATGCTTACGATTGCCGTATGGGGTGTTATCGCCTTATCCGGCCTACGGATGGCGGGTGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                           <  2:406132/145‑1 (MQ=11)
                                                                                            gACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCt                                                                                                            <  2:24881/96‑1 (MQ=33)
                                                                                             aCGCGTTAGCGTCGCATCAGGCTTCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                     >  1:750057/1‑150 (MQ=37)
                                                                                              cgcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                    >  1:1125042/1‑150 (MQ=38)
                                                                                                cgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGcc                                                                                                                                                  <  2:1516324/54‑1 (MQ=35)
                                                                                                    aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                                                               <  1:589064/117‑1 (MQ=32)
                                                                                                       gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACTCCGTATCTGGCCTACGGATTGCGCGGGAATTTGTAGTCCTGATAAGACGCGTTAGCGTCGCATCCGGCATCTGCACACGACTGCAGGATGGGATTAACGTCTTGt                                           >  2:1096315/1‑150 (MQ=14)
                                                                                                            aTCAGGCATCTGCGCACTACTGCTGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGTGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGc                                      <  2:614991/150‑1 (MQ=255)
                                                                                                             tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                     <  1:222692/150‑1 (MQ=255)
                                                                                                              cAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCt                                    >  2:1190278/1‑150 (MQ=255)
                                                                                                              cAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCt                                    >  2:1168214/1‑150 (MQ=255)
                                                                                                              cAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCAAACGACTGCCGGATGAGATAAACGTCTTGTCCGGCCt                                    >  2:539995/1‑150 (MQ=255)
                                                                                                                     cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTc                             <  1:220310/150‑1 (MQ=255)
                                                                                                                       gcgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                   <  1:1081442/78‑1 (MQ=34)
                                                                                                                          cACGACTGCCGGATGCGGCGTTAACGCATTATCCGGCCTACGGATGGCGCGGGATTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCg                        <  2:800947/150‑1 (MQ=255)
                                                                                                                          cACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCg                        >  2:237774/1‑150 (MQ=255)
                                                                                                                           ccgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGt                       >  2:1351024/7‑150 (MQ=17)
                                                                                                                           aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGt                       <  2:1724873/150‑1 (MQ=255)
                                                                                                                              aCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                    <  1:787149/150‑1 (MQ=255)
                                                                                                                              aCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                    >  1:408519/1‑150 (MQ=255)
                                                                                                                               cTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGAGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                   >  2:572282/1‑150 (MQ=255)
                                                                                                                                 gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                 <  2:579276/150‑1 (MQ=255)
                                                                                                                                 gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                 <  1:1005054/150‑1 (MQ=255)
                                                                                                                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                 >  1:1648087/1‑150 (MQ=18)
                                                                                                                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTCCGCCCGTAGGCAg                 >  2:484395/1‑150 (MQ=18)
                                                                                                                                   cGGATGGGGCGTATACTCAGTTTCCGGCTTTTGTTTGGCGCGTGAATTTGTAGGACTTATAAGACGCGTTAGCTTCGGATAAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTc               <  2:1604702/150‑1 (MQ=255)
                                                                                                                                    ggataaggcgTTTACGCCGCATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa              >  1:832893/7‑150 (MQ=18)
                                                                                                                                      aTGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGTGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCAtt            <  2:608809/150‑1 (MQ=255)
                                                                                                                                       tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTa           <  2:1315325/150‑1 (MQ=255)
                                                                                                                                        gcggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaa          <  2:1439004/150‑1 (MQ=255)
                                                                                                                                         cggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa         >  1:1390756/1‑150 (MQ=255)
                                                                                                                                          ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt        >  1:799935/1‑150 (MQ=255)
                                                                                                                                          ggcgTAAACGCCTTAGCCGGCCTACGGAGGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt        <  1:749439/150‑1 (MQ=255)
                                                                                                                                          ggcgTAAACGACTTATCCGGCCTACGGATGGCGCGGTAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt        <  2:244552/150‑1 (MQ=255)
                                                                                                                                            cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAg      <  1:1018053/150‑1 (MQ=255)
                                                                                                                                             ttAAACGCCTTATCCGGCCTTCGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCTGTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGt     <  2:776806/149‑1 (MQ=255)
                                                                                                                                               aaaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTc   >  2:687962/1‑150 (MQ=255)
                                                                                                                                               aaaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTc   <  1:1276131/150‑1 (MQ=255)
                                                                                                                                               aaaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTc   <  1:1037450/150‑1 (MQ=255)
                                                                                                                                                aaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAGAGTTCt  >  1:194986/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCT  >  NC_000913/375104‑375397

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