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.G92.5% 699.6 / 47.5 228intergenic (‑242/+111)glmU/atpCN‑acetyl glucosamine‑1‑phosphate uridyltransferase/membrane‑bound ATP synthase, F1 sector, epsilon‑subunit
Reads supporting (aligned to +/- strand):  ref base . (9/8);  new base G (72/139);  total (81/147)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 267 total aligned reads displayed.

CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTT  >  NC_000913/3912978‑3913160
                                                                                             |                                                                                           
cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACa                                                                                              >  1:3684437/1‑93 (MQ=255)
cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACa                                                                                              >  1:1695529/1‑93 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                           >  2:4404410/1‑93 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                           <  1:3647204/93‑1 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                           >  1:4655809/1‑93 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                           >  1:1481447/1‑93 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑g                                                                                           >  2:3485950/1‑93 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gg                                                                                          <  1:4721707/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gg                                                                                          <  2:182860/93‑1 (MQ=255)
    tGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggg                                                                                        <  2:1991497/93‑1 (MQ=255)
     gTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑ggggg                                                                                       <  1:3425549/93‑1 (MQ=255)
      tCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggggg                                                                                      >  1:4749041/1‑93 (MQ=255)
       cTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑ggggggg                                                                                     <  2:3161935/93‑2 (MQ=255)
        tGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑gggggggt                                                                                    >  2:2091896/1‑91 (MQ=255)
          cAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGgttt                                                                                  >  1:3962486/1‑93 (MQ=255)
          cAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGgttt                                                                                  <  1:3150133/93‑1 (MQ=255)
             aaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtat                                                                                >  1:2151261/1‑93 (MQ=255)
             aaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGAGCAATGAATAAACAGGAAAACAGGGGGGGTTTtat                                                                                <  2:2974744/93‑1 (MQ=255)
                aTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATc                                                                             <  1:2324392/93‑1 (MQ=255)
                aTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATc                                                                             <  1:2650456/93‑1 (MQ=255)
                   aaGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGGTTTTATATcag                                                                           <  1:1578312/93‑1 (MQ=255)
                   aaGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagc                                                                          <  2:1325621/93‑1 (MQ=255)
                   aaGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagc                                                                          >  1:2283318/1‑93 (MQ=255)
                   aaGGCAGTACATTAATTTAATGGTCAGGGGAATGTTTTTGCTCTTATTTCGATAAATGAATAATCAGGAAAACAGGGGGGGTTTTATATcagc                                                                          <  2:2569173/93‑1 (MQ=255)
                     ggCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagcag                                                                        <  2:1180793/93‑1 (MQ=255)
                     ggCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagcag                                                                        >  2:2348469/1‑93 (MQ=255)
                       cAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGa                                                                      <  1:1549599/93‑1 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGAt                                                                     <  1:40512/93‑1 (MQ=255)
                          tACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCt                                                                   <  2:2744375/93‑1 (MQ=255)
                           aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTa                                                                  <  2:2381965/93‑1 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                                 <  2:3470549/93‑1 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                                 >  2:4046943/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                                >  1:1430340/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                                >  2:3195890/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACACGGGGGGTTTTATATCAGCACGATCTAtg                                                                >  2:2460498/1‑93 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                               <  1:823115/93‑1 (MQ=255)
                                aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGa                                                             >  1:1000512/1‑93 (MQ=255)
                                aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTATTATATCAGCAGGATCTATGTGa                                                             >  2:2593981/1‑93 (MQ=255)
                                 aTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGaa                                                            >  1:3900269/1‑93 (MQ=255)
                                 aTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGaa                                                            <  1:2105097/93‑1 (MQ=255)
                                  tttAAGCGTCGGGTGTATGTTTTTGCTCTTATTGCGATCAATGATTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAAc                                                           <  1:869977/93‑1 (MQ=255)
                                  tttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAAc                                                           <  2:2893005/93‑1 (MQ=255)
                                    tAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGc                                                         <  2:360810/93‑1 (MQ=255)
                                      aGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTa                                                       >  2:4072405/1‑93 (MQ=255)
                                       gCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAt                                                      <  1:3758012/93‑1 (MQ=255)
                                       gCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAt                                                      <  1:2935370/93‑1 (MQ=255)
                                         gTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTc                                                    <  2:1562184/93‑1 (MQ=255)
                                         gTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTc                                                    <  2:520615/93‑1 (MQ=255)
                                         gTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTc                                                    <  2:1556474/93‑1 (MQ=255)
                                         gTCAGGTGGATGTTGTTACACTTATTTCGATCAATGAATAAACAGGAAAGCAGGGGGGGTTTTATTTCAGCAGGATCTATATGAACGCTATTc                                                    <  2:3467446/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATGAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                   <  1:1088203/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                   <  2:1664512/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                   <  1:2462485/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                   <  1:2772605/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                   <  1:2881195/93‑1 (MQ=255)
                                             ggTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                                >  2:72806/1‑93 (MQ=255)
                                             ggTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                                <  2:3654550/93‑1 (MQ=255)
                                             ggTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                                <  1:456129/93‑1 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                               <  2:1055732/93‑1 (MQ=255)
                                               tGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACg                                              <  1:4136799/93‑1 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                             >  2:1388422/1‑93 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                             <  2:2616530/93‑1 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                             <  2:4346404/93‑1 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                             >  1:4421449/1‑93 (MQ=255)
                                                  aTGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGt                                           <  1:6541/93‑1 (MQ=255)
                                                     tttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                                        <  2:262506/93‑1 (MQ=255)
                                                       tttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCacac                                      <  2:2047869/93‑1 (MQ=255)
                                                          gCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACAcgcg                                   >  1:1892672/1‑93 (MQ=255)
                                                               tATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaa                              >  2:2022669/1‑93 (MQ=255)
                                                                aTTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaa                             <  2:1216791/93‑1 (MQ=255)
                                                                    cGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag                         <  1:2979976/93‑1 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                        >  1:4458148/1‑93 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                        <  2:611245/93‑1 (MQ=255)
                                                                      aTCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcc                       >  1:2589624/1‑93 (MQ=255)
                                                                        cAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccag                     <  2:2183948/93‑1 (MQ=255)
                                                                        cAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccag                     >  2:3901506/1‑93 (MQ=255)
                                                                        cAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAagccagc                    <  1:1408171/93‑1 (MQ=255)
                                                                         aaTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc                    <  2:3762058/93‑1 (MQ=255)
                                                                         aaTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc                    >  2:3037163/1‑93 (MQ=255)
                                                                         aaTGAATAAACAGGAAAACAGGGGGGGGTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc                    >  2:805134/1‑93 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                   >  2:4607590/1‑93 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                   >  2:1725501/1‑93 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                   <  2:3858733/93‑1 (MQ=255)
                                                                           tGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCt                  <  2:1869549/93‑1 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt                >  1:3357729/1‑93 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt                <  1:3549581/93‑1 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCGGt                >  2:2515241/1‑93 (MQ=255)
                                                                              aTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACTGGTCACACGCGCAAAAAAAAGCACGCCTggt               >  2:3953159/1‑91 (MQ=255)
                                                                              aTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGtt               <  1:1240482/93‑1 (MQ=255)
                                                                              aTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGtt               <  1:3079904/93‑1 (MQ=255)
                                                                               tAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGttt              >  2:3208888/1‑93 (MQ=255)
                                                                                aaaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc             >  2:2264058/1‑93 (MQ=255)
                                                                                 aaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTcc            <  1:3495032/93‑1 (MQ=255)
                                                                                 aaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTcc            <  1:1949691/93‑1 (MQ=255)
                                                                                   cAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAg          >  1:2393841/1‑93 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa         >  1:43423/1‑93 (MQ=255)
                                                                                     ggAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGAc        <  2:1242644/93‑1 (MQ=255)
                                                                                     ggAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGAc        <  1:1753319/93‑1 (MQ=255)
                                                                                        aaaCAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTgg     >  2:4019255/1‑93 (MQ=255)
                                                                                        aaaCAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTgg     <  2:146429/93‑1 (MQ=255)
                                                                                             |                                                                                           
CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACA‑‑GGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGGCTT  >  NC_000913/3912978‑3913160

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