Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,792,808 T→A 6.8% I25N (ATT→AAT)  gabT → 4‑aminobutyrate aminotransferase GabT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,792,8080TA 6.8% 185.3 / 3.5 59I25N (ATT→AAT) gabT4‑aminobutyrate aminotransferase GabT
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base A (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TCGAAGTATGGCATCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGTTGAAAAAACTGT  >  NC_000913/2792665‑2792954
                                                                                                                                               |                                                                                                                                                  
tCGAAGTATGGCATCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGc                                                                                                                                              <  1:1090699/150‑1 (MQ=255)
tCGAAGTATGGCATCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGc                                                                                                                                              >  2:255649/1‑150 (MQ=255)
      tATGGCATCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACcg                                                                                                                                        <  2:531059/150‑1 (MQ=255)
          gCATCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACcgcgcg                                                                                                                                    >  1:1303837/1‑150 (MQ=255)
             tCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGaa                                                                                                                                 >  1:51085/1‑150 (MQ=255)
                  gATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTg                                                                                                                            <  2:399277/150‑1 (MQ=255)
                   aTTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGc                                                                                                                           >  2:422452/1‑150 (MQ=255)
                   aTTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGc                                                                                                                           <  2:1379696/150‑1 (MQ=255)
                       cTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCggg                                                                                                                       >  1:1266148/1‑150 (MQ=255)
                       cTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCggg                                                                                                                       >  1:896355/1‑150 (MQ=255)
                           gAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGtgtg                                                                                                                   <  1:621557/150‑1 (MQ=255)
                            aaaTCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTgg                                                                                                                  <  2:995273/150‑1 (MQ=255)
                              aTCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGGTTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGa                                                                                                                <  1:143565/150‑1 (MQ=255)
                              aTCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGa                                                                                                                >  2:869771/1‑150 (MQ=255)
                              aTCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGa                                                                                                                >  2:174804/1‑150 (MQ=255)
                                      tatGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAg                                                                                                        <  1:1292484/150‑1 (MQ=255)
                                        tgtgCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCACTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGc                                                                                                      >  2:402543/1‑150 (MQ=255)
                                              tCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAg                                                                                                >  1:576155/1‑150 (MQ=255)
                                                   cTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCt                                                                                           >  1:333536/1‑150 (MQ=255)
                                                     ttAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGTGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTg                                                                                         >  2:1181055/1‑150 (MQ=255)
                                                            gagaATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTcgc                                                                                  <  1:982467/150‑1 (MQ=255)
                                                              gaATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCgcgg                                                                                <  1:441113/150‑1 (MQ=255)
                                                                 tGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCgcgggcg                                                                             >  2:192495/1‑150 (MQ=255)
                                                                    gAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGa                                                                          >  1:8094/1‑150 (MQ=255)
                                                                        gAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGc                                                                      >  2:731354/1‑150 (MQ=255)
                                                                         aaCAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCg                                                                     >  1:185271/1‑150 (MQ=255)
                                                                          aCAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCgg                                                                    <  1:174804/150‑1 (MQ=255)
                                                                          aCAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCgg                                                                    <  1:255649/150‑1 (MQ=255)
                                                                             gCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGc                                                                 <  1:748877/150‑1 (MQ=255)
                                                                             gCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGc                                                                 >  1:1400341/1‑150 (MQ=255)
                                                                             gCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGc                                                                 <  1:280137/150‑1 (MQ=255)
                                                                              cAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCt                                                                <  2:8094/150‑1 (MQ=255)
                                                                                aTAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCa                                                              >  2:839875/1‑150 (MQ=255)
                                                                                        ttAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCggg                                                      <  2:576155/150‑1 (MQ=255)
                                                                                            tGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCAcc                                                  <  2:565991/150‑1 (MQ=255)
                                                                                             gCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCt                                                 >  1:1061849/1‑150 (MQ=255)
                                                                                                  gccgcAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGCCCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATc                                            >  1:1161069/1‑150 (MQ=255)
                                                                                                          cAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGAATTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAgg                                      <  2:752693/148‑1 (MQ=255)
                                                                                                          cAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGAATTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAgg                                      >  1:752693/1‑148 (MQ=255)
                                                                                                          cAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGAATTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAgg                                      >  1:747798/1‑148 (MQ=255)
                                                                                                          cAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGAATTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAgg                                      <  2:747798/148‑1 (MQ=255)
                                                                                                           aGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAggtgg                                   >  1:381013/1‑150 (MQ=255)
                                                                                                           aGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCACTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAggtgg                                   <  1:402543/150‑1 (MQ=255)
                                                                                                               gATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGc                               <  1:731354/150‑1 (MQ=255)
                                                                                                                    cccGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCgg                          >  2:873801/1‑150 (MQ=255)
                                                                                                                    cccGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCgg                          <  2:1303837/150‑1 (MQ=255)
                                                                                                                     ccGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGt                         >  1:574448/1‑150 (MQ=255)
                                                                                                                        tGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGa                      <  1:873801/150‑1 (MQ=255)
                                                                                                                          gCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAg                    >  1:919866/1‑150 (MQ=255)
                                                                                                                           cGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAgc                   >  1:1366458/1‑150 (MQ=255)
                                                                                                                           cGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAgc                   >  1:771543/1‑150 (MQ=255)
                                                                                                                               gggCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAg               >  2:1006125/1‑150 (MQ=255)
                                                                                                                                ggCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGt              >  1:200332/1‑150 (MQ=255)
                                                                                                                                 gCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGtt             <  2:896355/150‑1 (MQ=255)
                                                                                                                                  cAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGTTg            <  1:869771/150‑1 (MQ=255)
                                                                                                                                      ttCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGTTGaaaa        <  1:192495/150‑1 (MQ=255)
                                                                                                                                         aCCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTg                        >  2:698479/1‑131 (MQ=255)
                                                                                                                                         aCCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTg                        <  1:698479/131‑1 (MQ=255)
                                                                                                                                            cGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATGCCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGTTGAAAAAACTGt  >  1:1174551/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCGAAGTATGGCATCGAAGATTACTTAGAAATCAAATATATGTGCATCGGTCTTTAACTGGAGAATGCGAATGAACAGCAATAAAGAGTTAATGCAGCGCCGCAGTCAGGCGATTCCCCGTGGCGTTGGGCAAATTCACCCGATTTTCGCTGACCGCGCGGAAAACTGCCGGGTGTGGGACGTTGAAGGCCGTGAGTATCTTGATTTCGCGGGCGGGATTGCGGTGCTCAATACCGGGCACCTGCATCCGAAGGTGGTGGCCGCGGTGGAAGCGCAGTTGAAAAAACTGT  >  NC_000913/2792665‑2792954

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