Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,455 +G intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G98.3% 215.9 / ‑4.2 59intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (25/33);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTT  >  NC_000913/3560305‑3560597
                                                                                                                                                        |                                                                                                                                              
gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGc                                                                                                                                                 >  2:12882/1‑151 (MQ=255)
     aatgtcaaaaaaCAGCGCCGACCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCGCGGGGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGA‑TGGcgcagg                                                                                                                                            <  1:394497/142‑6 (MQ=255)
       tatCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGcgcag                                                                                                                                             >  1:181056/1‑144 (MQ=255)
        atCGAAAAACGGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGAt                                                                                                                                          <  1:336107/150‑1 (MQ=255)
          cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATc                                                                                                                                         >  1:855634/1‑149 (MQ=255)
              aaccggCGTCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCTCCGCCATGA‑TGGCGCAGGATCAgg                                                                                                                                      <  2:473244/143‑1 (MQ=255)
                     gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                            >  1:680702/1‑151 (MQ=255)
                      tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                           >  1:211335/1‑151 (MQ=255)
                        gggCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGc                                                                                                                                <  2:180670/142‑1 (MQ=255)
                                     tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                            >  2:514925/1‑151 (MQ=255)
                                       tggtCCGCCACTTTGCGGGCGACTCGCTCTTTTTCTTCGGTCTGGGTCGCCTTGCGATCGTGCCACGGCGAGTTAACGTAACTGGAAGGCAGCGCCGCACCGCCATGAGGGGG‑CAGGACCATGTTTTGCTTCTCCCGCCCATCGAGGTtgc                                                                                                          >  1:129106/4‑151 (MQ=255)
                                        gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                         >  2:272596/1‑151 (MQ=255)
                                         cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                        >  1:102720/1‑151 (MQ=255)
                                           ccGCCCCTTTCGGGGCGATGCACTCTTTTTCTTCGGTTTGGGTGCCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCGCCCCCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                      <  2:1116636/151‑1 (MQ=255)
                                           ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                      >  2:1203520/1‑151 (MQ=255)
                                           ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                        <  2:1213022/149‑1 (MQ=255)
                                             gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                     >  2:656367/1‑150 (MQ=255)
                                                   ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                               <  2:1024372/150‑1 (MQ=255)
                                                       gggCGTTGCGCTCTTTTTCTTCGGTCGGGGTGCCCTTGCGGTGGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                          <  1:81579/151‑1 (MQ=255)
                                                       gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                                                                           <  1:55286/150‑1 (MQ=255)
                                                         gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATGGTGCCACGGCGTGTTAACAGAGCTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                        <  1:317382/151‑1 (MQ=255)
                                                            aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                      >  2:1009550/1‑150 (MQ=255)
                                                             tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                    >  1:678541/1‑151 (MQ=255)
                                                              gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                   >  1:1150510/1‑151 (MQ=255)
                                                               cgcTCTTTTTCTCGGATCTGGTTGGCTTTGGCTTCGCCCCACGCCGTGTTAACCAAACGGGAGGGCAGCACCGCACCGCCATGA‑TGGCGCAGTATCAGGTTTTGCTACGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa                                                                                     <  1:752218/148‑1 (MQ=255)
                                                                    tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                              >  1:83113/1‑150 (MQ=255)
                                                                     ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                            >  2:851907/1‑151 (MQ=255)
                                                                      tttCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                            <  2:742723/150‑1 (MQ=255)
                                                                       ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                                                                          >  2:854087/1‑151 (MQ=255)
                                                                       ttctgggttCTGGGTGGCCTGGCAATCATGCCGCGGCGTGTTAACCGAACTGGAGGGCAGGGCGTCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGGGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                            <  2:88976/141‑1 (MQ=255)
                                                                        tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                         >  1:1016051/1‑151 (MQ=255)
                                                                         cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                        >  1:1111183/1‑151 (MQ=255)
                                                                          ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctct                                                                       >  1:275316/1‑151 (MQ=255)
                                                                               tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                    >  2:169974/1‑149 (MQ=255)
                                                                                  gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                                <  1:1068064/150‑1 (MQ=255)
                                                                                  gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                 <  2:179369/149‑1 (MQ=255)
                                                                                   ggTGGCCTTTCGACCGTGCCACGGCGTGTTAACCGAATTGGAAGGCAGCGCCGCCCCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                               <  2:187610/150‑1 (MQ=255)
                                                                                      ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                           >  2:422358/1‑151 (MQ=255)
                                                                                        ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                         <  2:1104676/151‑1 (MQ=255)
                                                                                         cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                         >  1:824388/1‑150 (MQ=255)
                                                                                                  gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                <  1:827884/150‑1 (MQ=255)
                                                                                                      cccGGCGTTTTAACCAACCTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                           <  1:1176813/149‑1 (MQ=255)
                                                                                                       aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                          >  1:195528/1‑151 (MQ=255)
                                                                                                        cGGCGAGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                         >  1:344010/1‑151 (MQ=255)
                                                                                                            gtgtAAACCGAACTGGAAGGATGCCCCACACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                       <  1:444560/149‑1 (MQ=255)
                                                                                                                          ggAAGGCAGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                        <  2:268144/150‑1 (MQ=255)
                                                                                                                          ggAAGCCGGCGCCGCTCCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTAGata                       <  1:268882/151‑1 (MQ=255)
                                                                                                                            aaGGCAGCGCTGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                     <  2:181062/151‑1 (MQ=255)
                                                                                                                                 aGCGCCGCACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                >  2:880314/1‑151 (MQ=255)
                                                                                                                                   cgccgcACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt               <  1:135048/150‑1 (MQ=255)
                                                                                                                                   cgccgcACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                <  2:62637/149‑1 (MQ=255)
                                                                                                                                     ccgcACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg            >  1:50046/1‑151 (MQ=255)
                                                                                                                                     ccgcACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt             <  1:669199/150‑1 (MQ=255)
                                                                                                                                        cACCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg         >  2:918819/1‑151 (MQ=255)
                                                                                                                                         aCCGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa           <  2:491240/148‑1 (MQ=255)
                                                                                                                                          ccGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt       <  1:710930/151‑1 (MQ=255)
                                                                                                                                          ccGCCATGA‑TGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGCCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAGCCCTGCTGTTTAACCAGTTCGATAATGCCGTTGTGACGt        <  2:582303/150‑1 (MQ=255)
                                                                                                                                           cGCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg      <  2:173741/151‑1 (MQ=255)
                                                                                                                                            gCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt     <  2:617672/151‑1 (MQ=255)
                                                                                                                                            gCCATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt     <  2:1075069/151‑1 (MQ=255)
                                                                                                                                              cATGA‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt   <  1:272589/151‑1 (MQ=255)
                                                                                                                                               atga‑tgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt  <  1:684319/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                              
GTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑TGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTT  >  NC_000913/3560305‑3560597

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