Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,644 A→G W15R (TGG→CGG)  lplA ← lipoate‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,6440AG99.7% 1225.8 / ‑5.2 379W15R (TGG→CGG) lplAlipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (170/208);  total (171/208)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 399 total aligned reads displayed.

CTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTT  >  NC_000913/4621554‑4621733
                                                                                          |                                                                                         
 tacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCgg                                                                                        >  1:4581264/1‑93 (MQ=255)
 tacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCgg                                                                                        >  2:7699350/1‑93 (MQ=255)
  acCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCggg                                                                                       >  2:4845142/1‑93 (MQ=255)
  acCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCggg                                                                                       <  1:2365210/93‑1 (MQ=255)
   ccGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGt                                                                                      <  1:7209083/93‑1 (MQ=255)
    cGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTc                                                                                     <  1:5797184/93‑1 (MQ=255)
    cGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTc                                                                                     <  2:6749148/93‑1 (MQ=255)
       gtCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTa                                                                                  >  1:1906507/1‑93 (MQ=255)
        tCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGATAAACCGCGGGTCGTaa                                                                                 >  2:5059818/1‑93 (MQ=255)
         cGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAag                                                                                <  1:9134466/93‑1 (MQ=255)
         cGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAag                                                                                <  2:632653/93‑1 (MQ=255)
              ttGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCa                                                                           >  2:2503285/1‑93 (MQ=255)
              ttGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCa                                                                           >  2:7355500/1‑93 (MQ=255)
               tGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCag                                                                          >  2:1059239/1‑93 (MQ=255)
                gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCaga                                                                         <  2:4257356/93‑1 (MQ=255)
                 cgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCagag                                                                        <  1:1773183/93‑1 (MQ=255)
                  gcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTTGCGAAAAATAAACTCTTCCATCGCGAGGTTAAACCGCGGGTCGTAAGCGTCAgaga                                                                       >  1:7019607/1‑93 (MQ=255)
                  gcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAgaga                                                                       <  2:9131324/93‑1 (MQ=255)
                    cAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATg                                                                     <  2:7986975/93‑1 (MQ=255)
                      gagaAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAg                                                                   >  1:5198383/1‑93 (MQ=255)
                        aaaaCATATCGCGCTTCGTGGCGTGCATATGGCGAAAAATTCACTCTTCCACGGCCAGGTTAAACTGCGTGTCGTAAGAGTCAGAGTTGAGCa                                                                 <  2:6952050/92‑1 (MQ=255)
                         aaaCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAg                                                                <  2:5548429/93‑1 (MQ=255)
                          aaCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAgg                                                               >  1:6460397/1‑93 (MQ=255)
                          aaCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAgg                                                               <  2:3357604/93‑1 (MQ=255)
                           aCAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGc                                                              <  1:2203948/93‑1 (MQ=255)
                               aaCGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGGGTaa                                                          >  2:2643049/1‑93 (MQ=255)
                               aaCGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTaa                                                          >  1:8290791/1‑93 (MQ=255)
                               aaCGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTaa                                                          >  1:6325789/1‑93 (MQ=255)
                                 cgcgCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtg                                                        <  2:860510/93‑1 (MQ=255)
                                 cgcgCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtg                                                        <  1:7635330/93‑1 (MQ=255)
                                   cgcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgtg                                                      >  2:8768598/1‑93 (MQ=255)
                                    gcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTgg                                                     >  2:5406574/1‑93 (MQ=255)
                                    gcTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTgg                                                     >  1:4333169/1‑93 (MQ=255)
                                     cTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                    >  1:1973513/1‑93 (MQ=255)
                                     cTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                    <  1:5025939/93‑1 (MQ=255)
                                      tGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGAc                                                   <  1:2465122/93‑1 (MQ=255)
                                       gCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGAct                                                  >  2:3176765/1‑92 (MQ=255)
                                        cGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACAt                                                 <  1:3001772/93‑1 (MQ=255)
                                        cGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACAt                                                 <  1:1112151/93‑1 (MQ=255)
                                         gTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATa                                                <  2:3637633/93‑1 (MQ=255)
                                         gTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATa                                                >  2:980836/1‑93 (MQ=255)
                                           ggCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAAc                                              <  1:4257937/93‑1 (MQ=255)
                                            gCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGCGTCGGAGATGAGCAGGCGTAATGTGGACATCACg                                             >  2:2061983/1‑93 (MQ=255)
                                               ggCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTCATGTGGACATAACGAtt                                          >  2:7855715/1‑93 (MQ=255)
                                                  aTTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTcc                                       >  2:3864351/1‑93 (MQ=255)
                                                  aTTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTcc                                       >  2:616134/1‑93 (MQ=255)
                                                  aTTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTcc                                       <  2:8170767/93‑1 (MQ=255)
                                                  aTTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGGCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTcc                                       >  2:5109399/1‑93 (MQ=255)
                                                   tttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                      <  2:4354982/93‑1 (MQ=255)
                                                   tttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                      <  1:6683340/93‑1 (MQ=255)
                                                   tttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                      >  1:9620056/1‑93 (MQ=255)
                                                   tttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                      <  1:5106310/93‑1 (MQ=255)
                                                   tttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                      >  1:10018402/1‑93 (MQ=255)
                                                    ttGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCtt                                     >  1:1035379/1‑93 (MQ=255)
                                                     tGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCttt                                    >  1:4840512/1‑93 (MQ=255)
                                                       gCGAAAAATACACTCTTCCACCTCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTTGACATAACGATTTCCTTTCt                                  <  1:3427698/93‑1 (MQ=255)
                                                        cGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                 >  2:7226069/1‑93 (MQ=255)
                                                        cGAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                 <  1:7960920/93‑1 (MQ=255)
                                                         gAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTg                                <  2:8902193/93‑1 (MQ=255)
                                                         gAAAAATACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTg                                >  1:7371236/1‑93 (MQ=255)
                                                          aaaaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                               <  2:9537987/93‑1 (MQ=255)
                                                          aaaaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                               <  2:714286/93‑1 (MQ=255)
                                                          aaaaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                               >  1:5901721/1‑93 (MQ=255)
                                                          aaaaaTACACTCTACCACCGCCAGGTTACACCGCGGGTCGTAAGAGTCAGAGATGAGCAGACGTAATGTGGACATAACGATTTCCTTTCTTgt                               <  1:8976247/93‑1 (MQ=255)
                                                           aaaaTACACTCTGCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                              <  1:1402153/93‑1 (MQ=255)
                                                            aaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtag                             <  1:5506522/93‑1 (MQ=255)
                                                            aaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtag                             <  1:3529554/93‑1 (MQ=255)
                                                            aaaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtag                             >  2:1205781/1‑93 (MQ=255)
                                                             aaTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                            <  2:5758261/93‑1 (MQ=255)
                                                              aTACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagta                           <  2:4771177/93‑1 (MQ=255)
                                                               tACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTaa                          <  1:695790/93‑1 (MQ=255)
                                                               tACACTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTaa                          >  1:6007414/1‑93 (MQ=255)
                                                                acacTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAAt                         <  2:6246051/93‑1 (MQ=255)
                                                                 cacTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATg                        >  1:2847093/1‑93 (MQ=255)
                                                                  acTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAAttg                       >  2:2196724/1‑91 (MQ=255)
                                                                  acTCTTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATgg                       <  2:813501/93‑1 (MQ=255)
                                                                   ctctTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                      >  2:2763582/1‑93 (MQ=255)
                                                                   ctctTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                      >  2:4918197/1‑93 (MQ=255)
                                                                    tctTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATgggc                     >  1:8390566/1‑92 (MQ=255)
                                                                    tctTCCACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGt                     >  2:8434843/1‑93 (MQ=255)
                                                                     ctcccACCGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTc                    <  2:6826519/90‑1 (MQ=255)
                                                                           ccGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtt              <  2:2011635/93‑1 (MQ=255)
                                                                           ccGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtt              >  2:2585494/1‑93 (MQ=255)
                                                                           ccGCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtt              >  1:2320415/1‑93 (MQ=255)
                                                                             gCCAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt            >  1:9753910/1‑93 (MQ=255)
                                                                              ccAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt           >  2:3049886/1‑93 (MQ=255)
                                                                              ccAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt           >  1:3144035/1‑93 (MQ=255)
                                                                               cAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGTTTTCCTTTCTTTTAGTAATGGGTCACTCtttttt          >  1:3338275/1‑93 (MQ=255)
                                                                               cAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttttt          <  1:9919672/93‑1 (MQ=255)
                                                                               cAGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttttt          <  2:9043377/93‑1 (MQ=255)
                                                                                aGGTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtc         >  2:7523206/1‑93 (MQ=255)
                                                                                  gTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc       >  1:947527/1‑93 (MQ=255)
                                                                                  gTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc       <  1:8704229/93‑1 (MQ=255)
                                                                                  gTTAAACCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc       >  1:938165/1‑93 (MQ=255)
                                                                                     aaaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCt    <  1:5573376/93‑1 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   >  1:4004081/1‑93 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   >  2:5998603/1‑93 (MQ=255)
                                                                                      aaCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt   >  1:8853228/1‑93 (MQ=255)
                                                                                       aCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCttt  >  2:165839/1‑93 (MQ=255)
                                                                                       aCCGCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTGCTTGTAGTAATGGGTCACTCTTTTTTCTCGCttt  <  1:4088417/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTT  >  NC_000913/4621554‑4621733

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