Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,913,076 (G)6→7 intergenic (‑248/+105) glmU ← / ← atpC N‑acetyl glucosamine‑1‑phosphate uridyltransferase/membrane‑bound ATP synthase, F1 sector, epsilon‑subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,913,0701.G96.4% 975.2 / 25.5 306intergenic (‑242/+111)glmU/atpCN‑acetyl glucosamine‑1‑phosphate uridyltransferase/membrane‑bound ATP synthase, F1 sector, epsilon‑subunit
Reads supporting (aligned to +/- strand):  ref base . (7/4);  new base G (167/128);  total (174/132)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 349 total aligned reads displayed.

CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTG  >  NC_000913/3912978‑3913156
                                                                                             |                                                                                      
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAg                                                                                       >  2:1719692/1‑93 (MQ=255)
   ttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggg                                                                                     <  1:4385548/93‑1 (MQ=255)
    tGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggg                                                                                    <  2:2537690/93‑1 (MQ=255)
    tGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggg                                                                                    <  2:4103500/93‑1 (MQ=255)
     gTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggggg                                                                                   >  1:5080185/1‑93 (MQ=255)
     gTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggggg                                                                                   >  1:2240877/1‑93 (MQ=255)
      tCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggggg                                                                                  <  1:2511334/93‑1 (MQ=255)
      tCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggggg                                                                                  <  1:2456190/93‑1 (MQ=255)
       cTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggggggg                                                                                 <  1:1860159/93‑2 (MQ=255)
        tGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATAAATAAACAGGAAAACAgggggggt                                                                                <  1:3497074/93‑3 (MQ=255)
          cAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGgttt                                                                              <  2:4907271/93‑1 (MQ=255)
           aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGgtttt                                                                             >  1:1631610/1‑92 (MQ=255)
           aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGgtttt                                                                             >  1:4757919/1‑92 (MQ=255)
             aaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtat                                                                            >  1:5176312/1‑93 (MQ=255)
             aaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtat                                                                            >  1:4694770/1‑93 (MQ=255)
              aaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtata                                                                           >  1:722119/1‑93 (MQ=255)
              aaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtata                                                                           <  2:2074452/93‑1 (MQ=255)
               aaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTatat                                                                          <  1:6453616/93‑1 (MQ=255)
                 ttAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGCATAAACAGGACAACAGGGGGGGTTTTATATca                                                                        >  1:6639996/1‑93 (MQ=255)
                   aaGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagc                                                                      >  1:6284693/1‑93 (MQ=255)
                    aGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagca                                                                     >  1:1026256/1‑93 (MQ=255)
                      gCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAgg                                                                   >  1:5981589/1‑93 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGAt                                                                 >  2:4146741/1‑93 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGAt                                                                 >  1:4526755/1‑93 (MQ=255)
                         gTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATc                                                                <  2:5061788/93‑1 (MQ=255)
                          tACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCt                                                               >  2:654108/1‑93 (MQ=255)
                           aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTa                                                              >  1:939022/1‑93 (MQ=255)
                           aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTa                                                              >  2:2119688/1‑93 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                             <  2:4858567/93‑1 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                             >  2:4671451/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                            <  2:6186352/93‑1 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                            >  1:2366690/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                            <  1:382834/93‑1 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                            >  1:2116017/1‑93 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATTTCAGCAGGATCTAtgt                                                           >  1:6173322/1‑93 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                           >  1:3259111/1‑93 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                           <  2:2872446/93‑1 (MQ=255)
                                 aTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGaa                                                        >  2:752807/1‑93 (MQ=255)
                                  tttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAAc                                                       >  1:2623405/1‑93 (MQ=255)
                                   ttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACg                                                      >  1:2668374/1‑93 (MQ=255)
                                    tAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGc                                                     >  1:4008966/1‑93 (MQ=255)
                                      aGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTa                                                   >  2:3705539/1‑93 (MQ=255)
                                            aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                             >  2:2341805/1‑93 (MQ=255)
                                            aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                             >  2:6468088/1‑93 (MQ=255)
                                             ggTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                            <  2:5718742/93‑1 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                           >  1:4234819/1‑93 (MQ=255)
                                               tGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACg                                          <  2:3846881/93‑1 (MQ=255)
                                               tGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACg                                          <  2:6047242/93‑1 (MQ=255)
                                               tGGATGTTTTTGCTCTTATCTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACg                                          <  2:4303260/93‑1 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCCGGATCTATGTGAACGCTATTCAGGACgg                                         >  2:2944209/1‑93 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                         >  2:6337556/1‑93 (MQ=255)
                                                 gATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACggg                                        >  1:3189832/1‑93 (MQ=255)
                                                  aTGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGt                                       <  1:1705593/93‑1 (MQ=255)
                                                   tGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTc                                      >  2:3316716/1‑93 (MQ=255)
                                                   tGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTc                                      >  2:1181192/1‑93 (MQ=255)
                                                     tttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                                    >  1:1113068/1‑93 (MQ=255)
                                                     tttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                                    >  2:1584467/1‑93 (MQ=255)
                                                      ttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcaca                                   <  1:1775604/93‑1 (MQ=255)
                                                           ctctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACgcgc                              >  2:2438055/1‑93 (MQ=255)
                                                           ctctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTCATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACgcgc                              >  1:4150152/1‑93 (MQ=255)
                                                            tctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCa                             >  2:3920666/1‑93 (MQ=255)
                                                             ctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaa                            >  1:3299878/1‑93 (MQ=255)
                                                              ttATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaa                           >  2:1090154/1‑93 (MQ=255)
                                                               tATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaa                          <  2:314888/93‑1 (MQ=255)
                                                                  ttCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaa                       <  1:6537900/93‑1 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                      <  1:5675074/93‑1 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                      <  1:5416490/93‑1 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                      <  1:494978/93‑1 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                      >  1:482387/1‑93 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                      >  2:275764/1‑93 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                      >  2:5263649/1‑93 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                    >  1:128768/1‑93 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGACACACGCGCAAAAAAAagc                    >  2:435289/1‑93 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAAGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                    <  1:2002868/93‑1 (MQ=255)
                                                                       tCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                  <  2:3057091/93‑1 (MQ=255)
                                                                       tCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                  <  2:5840315/93‑1 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc               <  1:975522/93‑1 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc               >  2:1512035/1‑93 (MQ=255)
                                                                           tGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCt              >  2:675177/1‑93 (MQ=255)
                                                                           tGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCt              <  2:2787417/93‑1 (MQ=255)
                                                                           tGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCt              >  2:2213158/1‑93 (MQ=255)
                                                                           tGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCt              >  2:35420/1‑93 (MQ=255)
                                                                           ggAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCt              <  2:3488561/92‑1 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt            <  1:1527444/93‑1 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATACCAGCAGGATCTATGTGAACCCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt            <  2:3378078/93‑1 (MQ=255)
                                                                              aTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGtt           >  2:4472169/1‑93 (MQ=255)
                                                                               tAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGttt          >  1:50347/1‑93 (MQ=255)
                                                                               tAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGttt          >  2:1846828/1‑93 (MQ=255)
                                                                                aaaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc         <  1:5603095/93‑1 (MQ=255)
                                                                                aaaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCCATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc         <  1:3408878/93‑1 (MQ=255)
                                                                                 aaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTcc        >  2:4749142/1‑93 (MQ=255)
                                                                                  aCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCa       <  2:841374/93‑1 (MQ=255)
                                                                                  aCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCa       >  1:2801335/1‑93 (MQ=255)
                                                                                   cAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAg      >  2:4526340/1‑93 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa     >  2:6006476/1‑93 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa     >  2:3790424/1‑93 (MQ=255)
                                                                                     ggAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTGCAGAc    >  1:2087415/1‑93 (MQ=255)
                                                                                     ggAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGAc    >  2:1610583/1‑93 (MQ=255)
                                                                                     agAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGAc    <  2:5506279/92‑1 (MQ=255)
                                                                                      gAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACt   <  1:5900947/93‑1 (MQ=255)
                                                                                             |                                                                                      
CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTG  >  NC_000913/3912978‑3913156

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: