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.2% 844.4 / 42.2 275intergenic (‑242/+111)glmU/atpCN‑acetyl glucosamine‑1‑phosphate uridyltransferase/membrane‑bound ATP synthase, F1 sector, epsilon‑subunit
Reads supporting (aligned to +/- strand):  ref base . (11/5);  new base G (135/124);  total (146/129)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 309 total aligned reads displayed.

CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGG  >  NC_000913/3912978‑3913157
                                                                                             |                                                                                        
cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACa                                                                                           >  2:2817312/1‑93 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAg                                                                                         >  1:5760796/1‑93 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAg                                                                                         <  1:2518869/93‑1 (MQ=255)
 ttttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAg                                                                                         <  1:4871894/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgg                                                                                       <  2:204910/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgg                                                                                       <  1:977271/93‑1 (MQ=255)
  tttGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgg                                                                                       <  1:5325778/93‑1 (MQ=255)
        tGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggggggt                                                                                 <  2:3207507/93‑3 (MQ=255)
        tGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggggggt                                                                                 >  1:5530236/1‑91 (MQ=255)
           aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGgtttt                                                                              >  1:3320823/1‑92 (MQ=255)
            gAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGgtttta                                                                             >  1:1107390/1‑91 (MQ=255)
            gAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGgtttta                                                                             >  2:10259/1‑91 (MQ=255)
              aaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTtata                                                                            <  1:385981/93‑1 (MQ=255)
                aTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATc                                                                          <  2:5136525/93‑1 (MQ=255)
                 ttAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATca                                                                         >  2:3971913/1‑93 (MQ=255)
                   aaGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagc                                                                       >  2:3525428/1‑93 (MQ=255)
                     ggCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATcagcag                                                                     <  2:443601/93‑1 (MQ=255)
                      gCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAgg                                                                    <  1:2195426/93‑1 (MQ=255)
                       caatACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGCGGGGGTTTTATATCAGCAGGa                                                                   <  2:4370284/90‑1 (MQ=255)
                       cAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGa                                                                   >  1:3524759/1‑93 (MQ=255)
                       cAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGa                                                                   <  1:4228421/93‑1 (MQ=255)
                       cAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGa                                                                   >  1:3358975/1‑93 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGAt                                                                  >  1:6145816/1‑93 (MQ=255)
                        aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGAt                                                                  >  1:794334/1‑93 (MQ=255)
                         gTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATc                                                                 >  2:6178969/1‑93 (MQ=255)
                           aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTa                                                               >  1:2389082/1‑93 (MQ=255)
                           aCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTa                                                               >  1:755821/1‑93 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                              <  1:3281707/93‑1 (MQ=255)
                            ccTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAt                                                              >  1:4334273/1‑93 (MQ=255)
                             cTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATAAATAAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtg                                                             <  2:1163025/93‑1 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                            >  2:3450813/1‑93 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                            >  2:754319/1‑93 (MQ=255)
                              ttAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTAtgt                                                            >  1:6060727/1‑93 (MQ=255)
                                aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGa                                                          <  2:4749718/93‑1 (MQ=255)
                                aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGa                                                          <  2:4550046/93‑1 (MQ=255)
                                aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGa                                                          <  2:2675328/93‑1 (MQ=255)
                                aaTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGa                                                          >  2:2853817/1‑93 (MQ=255)
                                 aTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGaa                                                         <  2:1639060/93‑1 (MQ=255)
                                 aTTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGaa                                                         <  1:5105563/93‑1 (MQ=255)
                                   ttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACg                                                       <  1:4092755/93‑1 (MQ=255)
                                   ttAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACg                                                       <  1:4623761/93‑1 (MQ=255)
                                    tAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGc                                                      <  1:2330134/93‑1 (MQ=255)
                                    tAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGc                                                      <  2:6090718/93‑1 (MQ=255)
                                      aGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTa                                                    <  1:5381403/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                <  2:3651632/93‑1 (MQ=255)
                                          tCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCa                                                <  1:4715415/93‑1 (MQ=255)
                                            aGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                              >  1:1450740/1‑93 (MQ=255)
                                             ggTGGATGTTTTTGCTCTTAGTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGa                                             <  2:2915298/93‑1 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                            <  2:3936103/93‑1 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                            >  1:3636152/1‑93 (MQ=255)
                                              gTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGAc                                            >  1:5230145/1‑93 (MQ=255)
                                               tGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACg                                           >  2:327607/1‑93 (MQ=255)
                                                ggATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACgg                                          >  2:5506110/1‑93 (MQ=255)
                                                     tttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                                     >  1:681460/1‑93 (MQ=255)
                                                      ttttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcaca                                    <  1:5688856/93‑1 (MQ=255)
                                                       tttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCacac                                   >  1:5593825/1‑93 (MQ=255)
                                                       tttGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCacac                                   <  2:2159059/93‑1 (MQ=255)
                                                         tGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACAcgc                                 >  2:375828/1‑93 (MQ=255)
                                                          gCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACAcgcg                                >  1:6206726/1‑93 (MQ=255)
                                                          gCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACAcgcg                                >  1:2628818/1‑93 (MQ=255)
                                                           ctctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACgcgc                               >  1:1776898/1‑93 (MQ=255)
                                                            tctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCa                              >  1:6091028/1‑93 (MQ=255)
                                                            tctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCa                              >  1:4955462/1‑93 (MQ=255)
                                                             ctTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaa                             <  2:4832693/93‑1 (MQ=255)
                                                              ttATTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaa                            >  2:5311122/1‑93 (MQ=255)
                                                                aTTTCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaa                          <  2:1329045/93‑1 (MQ=255)
                                                                 tttCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaa                         <  1:4008790/93‑1 (MQ=255)
                                                                  ttCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaa                        <  1:1065849/93‑1 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa                       >  1:4746765/1‑93 (MQ=255)
                                                                   tCGATCAATGAATAAACAGGAAAACAGGGGGGGGTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaagaaaaa                       >  1:1275073/1‑93 (MQ=255)
                                                                    cGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag                      >  2:5420862/1‑93 (MQ=255)
                                                                    cGATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag                      >  1:3592696/1‑93 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                     >  1:2829085/1‑93 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                     <  1:3109465/93‑1 (MQ=255)
                                                                     gATCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc                     <  2:1479541/93‑1 (MQ=255)
                                                                      aTCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcc                    >  1:2106557/1‑93 (MQ=255)
                                                                      aTCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcc                    <  1:400083/93‑1 (MQ=255)
                                                                       tCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                   <  1:594768/93‑1 (MQ=255)
                                                                       tCAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                   <  1:1352303/93‑1 (MQ=255)
                                                                       tCAATGAATAAACAGGAAAACAGGGGGGGCTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                   >  1:200919/1‑93 (MQ=255)
                                                                       tCAATGAATAAACAGAAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca                   <  2:43280/93‑1 (MQ=255)
                                                                        cAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccag                  >  2:5969276/1‑93 (MQ=255)
                                                                        cAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccag                  >  2:2497326/1‑93 (MQ=255)
                                                                        cAATGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTCAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccag                  >  2:911287/1‑93 (MQ=255)
                                                                         aaTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc                 >  2:1785369/1‑93 (MQ=255)
                                                                         aaTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagc                 >  1:2290229/1‑93 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                <  1:4458306/93‑1 (MQ=255)
                                                                          aTGAATAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagccagcc                <  1:319511/93‑1 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt             >  2:5013156/1‑93 (MQ=255)
                                                                             aaTAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGt             >  2:5726108/1‑93 (MQ=255)
                                                                               tAAACAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGttt           <  2:4092724/93‑1 (MQ=255)
                                                                                aaaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc          <  1:5831728/93‑1 (MQ=255)
                                                                                aaaCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTc          >  1:370475/1‑93 (MQ=255)
                                                                                  aCAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCa        >  2:4883106/1‑93 (MQ=255)
                                                                                   cAGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAg       <  2:1090056/93‑1 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa      >  2:5873519/1‑93 (MQ=255)
                                                                                    aGGAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGa      >  2:1209657/1‑93 (MQ=255)
                                                                                      gAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACt    >  1:3930083/1‑93 (MQ=255)
                                                                                      gAAAACAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACt    <  1:3429692/93‑1 (MQ=255)
                                                                                       aaaaCAGGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAACCAGCCTGTTTTCAGACTg   >  1:2721812/1‑93 (MQ=255)
                                                                                             |                                                                                        
CTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCCAGACTGG  >  NC_000913/3912978‑3913157

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