Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,497,141 A→G A19A (GCA→GCG yhfL → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,497,1410AG99.0% 642.3 / ‑3.6 205A19A (GCA→GCGyhfLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (178/25);  total (180/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 211 total aligned reads displayed.

TATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG  >  NC_000913/3497052‑3497230
                                                                                         |                                                                                         
tatTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCgtgt                                                                                        >  1:2146614/1‑93 (MQ=255)
tatTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCgtgt                                                                                        <  2:6870886/93‑1 (MQ=255)
 attaattaATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTa                                                                                       <  1:2559415/93‑1 (MQ=255)
  ttaattaaTATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTAc                                                                                      >  1:900409/1‑93 (MQ=255)
    aattaatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTg                                                                                    >  1:5608571/1‑93 (MQ=255)
    aattaatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTg                                                                                    >  2:611546/1‑93 (MQ=255)
      ttaatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGt                                                                                  >  1:281207/1‑93 (MQ=255)
        aatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCa                                                                                >  2:6816357/1‑93 (MQ=255)
        aatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCa                                                                                >  1:6530645/1‑93 (MQ=255)
          tatataATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCatt                                                                              >  2:6009876/1‑92 (MQ=255)
           atataATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCatat                                                                             >  2:2298598/1‑93 (MQ=255)
            tataATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATAtt                                                                            >  1:1640344/1‑93 (MQ=255)
            tataATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGGACTGGTCATAtt                                                                            >  1:2931479/1‑93 (MQ=255)
            tataATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGGACTGGTCATAtt                                                                            >  2:7047244/1‑93 (MQ=255)
             ataataCATAGGGAATGTAAATGAACTAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTg                                                                           >  1:396351/1‑93 (MQ=255)
             ataataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTg                                                                           >  1:6626620/1‑93 (MQ=255)
              taataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGa                                                                          >  1:2463997/1‑93 (MQ=255)
              taataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGa                                                                          >  2:3508757/1‑93 (MQ=255)
               aataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCGACGTTAACTGCGTGTACTGGTCATATTGaa                                                                         >  1:3251800/1‑93 (MQ=255)
                 taCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGaaaa                                                                       <  1:5129244/93‑1 (MQ=255)
                   cATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAAcc                                                                     >  2:2524112/1‑93 (MQ=255)
                     tAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGt                                                                   >  1:2393594/1‑93 (MQ=255)
                     tAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGCACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGt                                                                   <  2:6257702/93‑1 (MQ=255)
                      aGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTg                                                                  >  2:7478323/1‑93 (MQ=255)
                       gggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGa                                                                 >  1:818015/1‑93 (MQ=255)
                       gggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGa                                                                 >  2:524434/1‑93 (MQ=255)
                        ggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGAt                                                                >  1:427762/1‑93 (MQ=255)
                        ggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGAt                                                                <  1:1525519/93‑1 (MQ=255)
                         gAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATa                                                               >  1:5448810/1‑93 (MQ=255)
                          aaTGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATaa                                                              <  1:6827051/93‑1 (MQ=255)
                          aaTGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATaa                                                              >  2:3943239/1‑93 (MQ=255)
                           aTGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAg                                                             >  2:3233302/1‑93 (MQ=255)
                            tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGa                                                            >  1:3395187/1‑93 (MQ=255)
                            tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGa                                                            >  2:3293075/1‑93 (MQ=255)
                            tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGa                                                            >  2:5726780/1‑93 (MQ=255)
                            tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGa                                                            >  2:7731585/1‑93 (MQ=255)
                             gTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGaa                                                           <  1:3872614/93‑1 (MQ=255)
                             gTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGGGCCGTACTGGGTACGTTAACTGCGTGGACGGGGCATATTGAAAACCGTGATAAGaa                                                           >  1:3091622/1‑93 (MQ=255)
                              tAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAAc                                                          >  1:7487026/1‑93 (MQ=255)
                              tAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAAc                                                          >  2:420920/1‑93 (MQ=255)
                              tAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAAc                                                          >  1:2589056/1‑93 (MQ=255)
                               aaaTGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACt                                                         >  2:3586990/1‑93 (MQ=255)
                                 aTGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGc                                                       >  1:6236638/1‑93 (MQ=255)
                                  tGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGct                                                      >  2:2158234/1‑93 (MQ=255)
                                    aaCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGctct                                                    >  1:7187742/1‑93 (MQ=255)
                                    aaCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGctct                                                    >  2:4752544/1‑93 (MQ=255)
                                     aCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCtt                                                   <  2:7870034/93‑1 (MQ=255)
                                     aCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCtt                                                   >  1:232448/1‑93 (MQ=255)
                                     aCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCtt                                                   >  1:6161922/1‑93 (MQ=255)
                                      cAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTa                                                  >  2:7359863/1‑93 (MQ=255)
                                       aaaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTAc                                                 >  2:2957421/1‑93 (MQ=255)
                                       aaaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTAc                                                 >  2:2343828/1‑93 (MQ=255)
                                        aaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACg                                                >  2:2886402/1‑93 (MQ=255)
                                        aaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACg                                                >  1:4339756/1‑93 (MQ=255)
                                        aaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACg                                                >  2:1409291/1‑93 (MQ=255)
                                          tttATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTAAGAc                                              >  1:364862/1‑93 (MQ=255)
                                           ttattaAAGTTGCACTGGTAGGGGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACt                                             >  1:3091840/1‑93 (MQ=255)
                                              ttaAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTAcc                                          >  1:6699365/1‑93 (MQ=255)
                                                 aaGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACctgc                                       >  1:3017620/1‑93 (MQ=255)
                                                  aGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACctgct                                      >  2:2363730/1‑93 (MQ=255)
                                                   gTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACctgctg                                     >  2:5784216/1‑93 (MQ=255)
                                                    ttGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCtgctgc                                    >  2:4879536/1‑93 (MQ=255)
                                                      gCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCAc                                  >  2:187497/1‑93 (MQ=255)
                                                      gCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCAc                                  <  1:7348082/93‑1 (MQ=255)
                                                      gCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCAc                                  >  1:7634516/1‑93 (MQ=255)
                                                       cACTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCAcc                                 <  2:3773126/93‑1 (MQ=255)
                                                        aCTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCAccc                                >  2:4871657/1‑93 (MQ=255)
                                                        aCTGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCAccc                                >  2:2437336/1‑93 (MQ=255)
                                                          tGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCgg                              >  2:4004339/1‑93 (MQ=255)
                                                          tGGTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCgg                              >  2:5047246/1‑93 (MQ=255)
                                                           ggTAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGc                             >  2:2811787/1‑93 (MQ=255)
                                                             tAGGTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCaa                           <  1:4875584/93‑1 (MQ=255)
                                                               ggTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatt                         >  1:3820801/1‑93 (MQ=255)
                                                               ggTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatt                         >  2:5651878/1‑93 (MQ=255)
                                                               ggTGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatt                         >  1:1157009/1‑93 (MQ=255)
                                                                 tGCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttc                       >  2:4582642/1‑93 (MQ=255)
                                                                  gCAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttct                      >  1:3970129/1‑93 (MQ=255)
                                                                   cAGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttcta                     >  1:1469820/1‑93 (MQ=255)
                                                                    aGTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctat                    >  1:2202565/1‑93 (MQ=255)
                                                                     gTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatt                   >  1:1747474/1‑93 (MQ=255)
                                                                     gTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatt                   >  1:1025997/1‑93 (MQ=255)
                                                                     gTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatt                   >  1:3676512/1‑93 (MQ=255)
                                                                     gTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatt                   >  1:7589313/1‑93 (MQ=255)
                                                                     gTACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatt                   >  2:1451182/1‑93 (MQ=255)
                                                                      tACTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctattt                  >  1:4107147/1‑93 (MQ=255)
                                                                       aCTGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatttc                 >  2:1193225/1‑93 (MQ=255)
                                                                        cTGGCTACGTTAACTGCATGTACTGGTCATATTGAACACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatttct                >  2:7546667/1‑93 (MQ=255)
                                                                         tGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAAtttctatttcta               >  1:3561623/1‑93 (MQ=255)
                                                                         tGGCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAAtttctatttcta               >  2:7075545/1‑93 (MQ=255)
                                                                          ggCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTaa              >  2:5158166/1‑93 (MQ=255)
                                                                          ggCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTaa              >  1:7496278/1‑93 (MQ=255)
                                                                          ggCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTaa              >  2:6430133/1‑93 (MQ=255)
                                                                          ggCTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTaa              >  1:4470659/1‑93 (MQ=255)
                                                                            cTACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTaaaa            >  2:6545964/1‑93 (MQ=255)
                                                                             tACGTTAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAAt           >  1:7406943/1‑93 (MQ=255)
                                                                                 ttAACTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCAtt       >  2:1732761/1‑93 (MQ=255)
                                                                                    aCTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGc    >  1:5017289/1‑93 (MQ=255)
                                                                                    aCTGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGc    <  2:5999913/93‑1 (MQ=255)
                                                                                      tGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCgg  >  2:2415096/1‑93 (MQ=255)
                                                                                      tGCGTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCgg  >  1:5155471/1‑93 (MQ=255)
                                                                                         |                                                                                         
TATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG  >  NC_000913/3497052‑3497230

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