Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 703,652 A→G intergenic (‑41/‑292) nagB ← / → nagE glucosamine‑6‑phosphate deaminase/N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913703,6520AG100.0% 239.8 / NA 63intergenic (‑41/‑292)nagB/nagEglucosamine‑6‑phosphate deaminase/N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/26);  total (37/26)

GACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTAGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGACAAAACTCACAAAAGACACGCGT  >  NC_000913/703512‑703798
                                                                                                                                            |                                                                                                                                                  
gACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCttt                                                                                                                                           <  1:362702/150‑1 (MQ=255)
    aTCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGc                                                                                                                                       <  2:173206/150‑1 (MQ=255)
    aTCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCg                                                                                                                                      <  1:360497/151‑1 (MQ=255)
    aTCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGAGTGGGCATCTTTAAGc                                                                                                                                       <  1:300993/150‑1 (MQ=255)
            tCGGTTTGAACGCATTGATACGATTGGCGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCt                                                                                                                                <  2:936880/149‑1 (MQ=255)
            tCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGGCCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAAcc                                                                                                                                 <  1:103353/148‑1 (MQ=255)
                tttGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGAtt                                                                                                                           <  2:102670/150‑1 (MQ=255)
                tttGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGAtt                                                                                                                           <  1:395162/150‑1 (MQ=255)
                      cGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAg                                                                                                                    >  2:441638/1‑151 (MQ=255)
                        caatGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTCCAGCGGTAGTCAGGGGGATCAGCCTCATTATTCCCTTCAATAAGTAACATGTAAGTCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGTGCAGaa                                                                                                                  >  1:92793/4‑150 (MQ=255)
                        cATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGa                                                                                                                   >  2:499377/1‑150 (MQ=255)
                        cATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGa                                                                                                                   >  1:307158/1‑150 (MQ=255)
                        cATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGa                                                                                                                   >  1:834511/1‑150 (MQ=255)
                          ttGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAg                                                                                                                >  1:571230/1‑151 (MQ=255)
                            gATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAgg                                                                                                               >  1:149049/1‑150 (MQ=255)
                                cGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCa                                                                                                          >  2:1009125/1‑151 (MQ=255)
                                        gATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGAc                                                                                                   >  2:40435/1‑150 (MQ=255)
                                           atGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGa                                                                                               >  2:236532/1‑151 (MQ=255)
                                                     gcCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAAt                                                                                     >  1:790623/1‑151 (MQ=255)
                                                        cATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGAAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGAt                                                                                  >  2:99585/1‑151 (MQ=255)
                                                          tttGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATaa                                                                                >  2:227085/1‑151 (MQ=255)
                                                              ccGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATa                                                                            <  1:2728/151‑1 (MQ=255)
                                                                gACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAg                                                                          >  1:826034/1‑151 (MQ=255)
                                                                 aCCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAg                                                                          <  1:815785/150‑1 (MQ=255)
                                                                  ccTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGtt                                                                        >  2:961705/1‑151 (MQ=255)
                                                                     gTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGtttt                                                                      >  1:228704/1‑150 (MQ=255)
                                                                     gTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGtttt                                                                      >  2:1028904/1‑150 (MQ=255)
                                                                     gTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGtttt                                                                      >  1:454410/1‑150 (MQ=255)
                                                                     gTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGtttt                                                                      >  2:444211/1‑150 (MQ=255)
                                                                        cAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTgg                                                                  >  1:519088/1‑151 (MQ=255)
                                                                             gTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAg                                                             >  2:604035/1‑151 (MQ=255)
                                                                               aGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGc                                                            >  2:1030336/1‑150 (MQ=255)
                                                                                 tCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGcc                                                           <  2:933231/149‑1 (MQ=255)
                                                                                       ggATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAg                                                   >  2:259123/1‑151 (MQ=255)
                                                                                        gATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAgg                                                  >  2:918974/1‑151 (MQ=255)
                                                                                        gATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAgg                                                  >  2:51747/1‑151 (MQ=255)
                                                                                        gATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAg                                                   <  2:834791/150‑1 (MQ=255)
                                                                                        gATCAGTCTCATTATTCACCTCAATAAGTAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAggg                                                 >  2:295212/1‑151 (MQ=255)
                                                                                         aTCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAAGCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAggg                                                 <  1:675700/151‑1 (MQ=255)
                                                                                         aTCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAggg                                                 <  2:74904/151‑1 (MQ=255)
                                                                                         aTCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAAGCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAggg                                                 <  1:1050984/151‑1 (MQ=255)
                                                                                             gTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGGATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTg                                             >  2:606896/1‑151 (MQ=255)
                                                                                             gTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTg                                             >  2:1021479/1‑151 (MQ=255)
                                                                                                tcATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGata                                          <  1:183065/151‑1 (MQ=255)
                                                                                                tcATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGata                                          >  1:518941/1‑151 (MQ=255)
                                                                                                tcATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGata                                          <  1:458137/151‑1 (MQ=255)
                                                                                                 cATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGataa                                         >  2:934510/1‑151 (MQ=255)
                                                                                                 cATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGataa                                         <  1:961362/151‑1 (MQ=255)
                                                                                                  attattCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGataat                                        <  2:210285/151‑1 (MQ=255)
                                                                                                   ttattCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGataata                                       >  1:254763/1‑151 (MQ=255)
                                                                                                     attCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGAtaataa                                      >  2:965459/1‑150 (MQ=255)
                                                                                                     attCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGAtaataa                                      >  2:614008/1‑150 (MQ=255)
                                                                                                       tCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGa                                   <  2:665538/151‑1 (MQ=255)
                                                                                                          ccTCAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGatat                                <  1:218630/151‑1 (MQ=255)
                                                                                                             cAATAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTg                             >  1:92282/1‑151 (MQ=255)
                                                                                                               aTAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGt                           <  2:657696/151‑1 (MQ=255)
                                                                                                               aTAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGt                           <  1:805735/151‑1 (MQ=255)
                                                                                                                tAAGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTg                          <  2:171475/151‑1 (MQ=255)
                                                                                                                 aaGTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGa                         <  1:806814/151‑1 (MQ=255)
                                                                                                                   gTAAAATGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGACa                       >  1:409000/1‑151 (MQ=255)
                                                                                                                      aaaTGTAAGCCGTTGGCGGATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTTAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGACaaaa                    >  1:712013/1‑151 (MQ=255)
                                                                                                                                   tGGCGGATTGGGTATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGACAAAACTCACAAAAGaca       <  2:61387/151‑1 (MQ=255)
                                                                                                                                        gATTGGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGACAAAACTCACAAAAGACACGCGt  >  1:831902/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GACGATCGGCAGTCGGTTTGAACGCATTGATACGATTGACGATATGGCGAGCAGCCCATTTGCCGACCTGTTCAGCGGTAGTCAGGGGGATCAGTCTCATTATTCACCTCAATAAGTAAAATGTAAGCCGTTGGCGGATTAGGCATCTTTAAGCGTAACCTGGATTTGCGCAGACAGGCGTCAATCCGACCTGATTTTTTGAATGATAAAATAAGTTTTCTGGTTTAGCCAGTAAAAGGGAGTGATAATAACGATATTTGGTGACAAAACTCACAAAAGACACGCGT  >  NC_000913/703512‑703798

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