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.G94.8% 1152.5 / 52.0 381intergenic (‑242/+111)glmU/atpCN‑acetyl glucosamine‑1‑phosphate uridyltransferase/membrane‑bound ATP synthase, F1 sector, epsilon‑subunit
Reads supporting (aligned to +/- strand):  ref base . (14/6);  new base G (174/187);  total (188/193)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 465 total aligned reads displayed.

CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGG  >  NC_000913/3912978‑3913157
                                                                                             |                                                                                        
cTTTTGTCTGCAGAAAATTAAGGCAGTACCCTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACa                                                                                           <  2:7729691/93‑1 (MQ=255)
cTTGTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACa                                                                                           <  1:6060203/93‑1 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                        <  1:632124/93‑1 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                        <  1:4385879/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gg                                                                                       <  1:1452857/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gg                                                                                       <  2:6814259/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gg                                                                                       <  1:2949983/93‑1 (MQ=255)
   ttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑ggg                                                                                      <  1:1065991/93‑1 (MQ=255)
   ttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑ggg                                                                                      <  2:4976158/93‑1 (MQ=255)
    tGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggg                                                                                     >  2:1214625/1‑93 (MQ=255)
    tGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggg                                                                                     <  1:4031550/93‑1 (MQ=255)
     ttcggCAGAAAATTAAGGCAGTACCTTAATTTAAGCGGCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑ggggg                                                                                    <  1:168304/89‑1 (MQ=255)
     gTCTGCAGAAAATTAAGGCAGTACCTTAATGTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑ggggg                                                                                    <  2:3140010/93‑1 (MQ=255)
      tCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggggg                                                                                   <  1:7862658/93‑1 (MQ=255)
      tCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggggg                                                                                   <  1:6048721/93‑1 (MQ=255)
            gAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGgtttta                                                                             >  1:3976966/1‑91 (MQ=255)
             aaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtat                                                                             >  1:4885520/1‑93 (MQ=255)
                aTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATc                                                                          <  1:4743096/93‑1 (MQ=255)
                   aaGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcgaa                                                                       >  1:6482851/1‑90 (MQ=255)
                    aGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagca                                                                      <  1:1963520/93‑1 (MQ=255)
                    aGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagca                                                                      <  2:6682756/93‑1 (MQ=255)
                    aGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagca                                                                      >  2:1714383/1‑93 (MQ=255)
                       cAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGa                                                                   <  1:3453640/93‑1 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGAt                                                                  >  1:520752/1‑93 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGTTTTATATCAGCAGGATc                                                                 <  2:3321216/93‑1 (MQ=255)
                         gTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATc                                                                 <  1:7320945/93‑1 (MQ=255)
                          tACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCt                                                                >  1:4678770/1‑93 (MQ=255)
                          tACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCt                                                                <  2:3751028/93‑1 (MQ=255)
                          tACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGTTTTATATCAGCAGGATCTa                                                               >  2:3483294/1‑93 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                              >  2:307580/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                             <  2:4099926/93‑1 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                            <  2:419741/93‑1 (MQ=255)
                               tAATTTAAGCGTCAGGCGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGATTTATATCAGCAGGATCTAtgtg                                                           >  2:2054005/1‑93 (MQ=255)
                               tAATTGAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAGTGAATAAACAGGAAAACAGGGGGGGGTTTATATCAGCAGGATCTAtgtg                                                           <  2:3917631/93‑1 (MQ=255)
                                  tttAAGCGTCAGGTGGATTTTTTTGCTCTTATTTCGATAAATGAATAAACAGAAAAACAGGGGGGGGTTTTATATCAGCAGGATCTATGTGaa                                                         <  1:6974819/93‑1 (MQ=255)
                                  tttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAAc                                                        >  1:4611835/1‑93 (MQ=255)
                                   ttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACg                                                       <  2:3590616/93‑1 (MQ=255)
                                    tAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGc                                                      <  2:2802890/93‑1 (MQ=255)
                                     aaGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCt                                                     <  2:1056446/93‑1 (MQ=255)
                                        cGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAtt                                                  <  2:5937483/93‑1 (MQ=255)
                                        cGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAtt                                                  <  2:2220472/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                <  1:6560033/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                >  2:1509905/1‑93 (MQ=255)
                                           cAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAg                                               <  1:6666626/93‑1 (MQ=255)
                                           cAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAg                                               >  2:1416165/1‑93 (MQ=255)
                                            aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                              >  2:4498983/1‑93 (MQ=255)
                                            aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                              <  2:7085226/93‑1 (MQ=255)
                                             ggTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                             >  1:211591/1‑93 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                            >  1:7543861/1‑93 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                            >  2:4655131/1‑93 (MQ=255)
                                              gTGGATGTGTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                            <  1:5860460/93‑1 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                          <  2:1895488/93‑1 (MQ=255)
                                                 gATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACggg                                         >  1:4986059/1‑93 (MQ=255)
                                                 gATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACggg                                         >  2:5429296/1‑93 (MQ=255)
                                                 gATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACggg                                         >  1:5076349/1‑93 (MQ=255)
                                                 gATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACggg                                         >  2:6023685/1‑93 (MQ=255)
                                                 gATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACggg                                         >  2:2468093/1‑93 (MQ=255)
                                                  aTGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGt                                        <  2:2903068/93‑1 (MQ=255)
                                                   tGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTc                                       <  1:2311778/93‑1 (MQ=255)
                                                    gTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTca                                      >  2:6281861/1‑93 (MQ=255)
                                                     tttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                                     >  1:6634329/1‑93 (MQ=255)
                                                     tttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                                     >  1:1330467/1‑93 (MQ=255)
                                                      ttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcaca                                    <  2:6341602/93‑1 (MQ=255)
                                                      ttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcaca                                    >  1:7294311/1‑93 (MQ=255)
                                                      ttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcaca                                    <  1:916835/93‑1 (MQ=255)
                                                       tttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCacac                                   <  1:5259722/93‑1 (MQ=255)
                                                       tttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCacac                                   <  1:1865118/93‑1 (MQ=255)
                                                        ttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCacac                                   >  1:6457913/1‑93 (MQ=255)
                                                        ttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACAcg                                  >  1:3114862/1‑93 (MQ=255)
                                                           ctctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGGCACACgcgc                               >  1:289247/1‑93 (MQ=255)
                                                              ttATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaa                            >  2:6123331/1‑93 (MQ=255)
                                                              ttATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaa                            >  2:3654395/1‑93 (MQ=255)
                                                                aTTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaa                          >  1:5791888/1‑93 (MQ=255)
                                                                 tttCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaa                         >  1:1169215/1‑93 (MQ=255)
                                                                  tacGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaa                        <  2:7251475/91‑1 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGGTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGGCACACGCGCaaaaaaaa                       >  2:6507699/1‑93 (MQ=255)
                                                                    cGATCAATGAATAAACAGGAAAACAGGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                       >  1:7679245/1‑93 (MQ=255)
                                                                    cGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag                      <  2:4862254/93‑1 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                     >  2:2605756/1‑93 (MQ=255)
                                                                      aTCAATGAATAAACAGGAAAACAGGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                     >  2:5565808/1‑93 (MQ=255)
                                                                      aTCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAACAAAagcc                    >  1:3373157/1‑93 (MQ=255)
                                                                      aTCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcc                    <  1:450110/93‑1 (MQ=255)
                                                                       tCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                   <  2:4289293/93‑1 (MQ=255)
                                                                         aaTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc                 >  2:4040/1‑93 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                >  1:5482484/1‑93 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                >  2:570949/1‑93 (MQ=255)
                                                                            gAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTg              >  2:919444/1‑93 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt             >  1:3731292/1‑93 (MQ=255)
                                                                              aTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGtt            >  1:2590839/1‑93 (MQ=255)
                                                                               tAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGttt           >  2:1568391/1‑93 (MQ=255)
                                                                                aaaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc          <  1:6918722/93‑1 (MQ=255)
                                                                                 aaCAGGAAAACAGGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc          >  2:7159346/1‑93 (MQ=255)
                                                                                 aaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTcc         >  2:5200697/1‑93 (MQ=255)
                                                                                  aCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCa        <  1:713050/93‑1 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAACGCCAGCCTGTTTcccgc      <  2:7022564/93‑4 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa      <  1:5963100/93‑1 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa      <  1:541272/93‑1 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa      <  1:3638371/93‑1 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa      >  2:755182/1‑93 (MQ=255)
                                                                                       aaaaCAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTg   <  1:6816326/93‑1 (MQ=255)
                                                                                             |                                                                                        
CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGG  >  NC_000913/3912978‑3913157

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: