Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,255 T→C 6.4% intergenic (+374/‑204) 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,2550TC 6.4% 262.1 / 6.3 95intergenic (+374/‑204)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base T (44/44);  new base C (4/2);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATT  >  NC_000913/375108‑375401
                                                                                                                                                   |                                                                                                                                                  
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:1724873/1‑150 (MQ=39)
                                           ccTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                       >  2:220310/1‑150 (MQ=39)
                                           ccTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                       >  1:579276/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)
                                                                                                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:1168214/1‑150 (MQ=255)
                                                                                                          cAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCt                                        >  2:1190278/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:408519/1‑150 (MQ=255)
                                                                                                                          aCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                        <  1:787149/150‑1 (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:1037450/150‑1 (MQ=255)
                                                                                                                                           aaaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTc       <  1:1276131/150‑1 (MQ=255)
                                                                                                                                            aaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAGAGTTCt      >  1:194986/1‑150 (MQ=255)
                                                                                                                                             aCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTg     >  1:1615008/1‑150 (MQ=255)
                                                                                                                                             aCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTg     <  1:237774/150‑1 (MQ=255)
                                                                                                                                             aCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGAGAAACGGCTAGTCCGGCCTACATTTCGCCCGTAGGCAGGCATTAAATAGTTCTg     >  2:1137246/1‑150 (MQ=255)
                                                                                                                                              cGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                                                          >  1:239899/1‑80 (MQ=34)
                                                                                                                                              cGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                                                          <  2:239899/80‑1 (MQ=34)
                                                                                                                                               gCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCCGTCATTAAATAGTTCTGAt   <  1:1465154/150‑1 (MQ=255)
                                                                                                                                               gCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGAt   <  1:1401234/150‑1 (MQ=255)
                                                                                                                                               gCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCGGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGAt   <  1:1316141/150‑1 (MQ=255)
                                                                                                                                                ccTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGAta  >  1:1040550/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATT  >  NC_000913/375108‑375401

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: