Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,624,717 A→G intergenic (‑220/‑56) ydeI ← / → ydeJ BOF family protein YdeI/conserved protein YdeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,624,7170AG88.1% 1250.0 / 134.8 388intergenic (‑220/‑56)ydeI/ydeJBOF family protein YdeI/conserved protein YdeJ
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base G (184/158);  total (204/184)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 388 total aligned reads displayed.

TATACAGAGAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCG  >  NC_000913/1624618‑1624815
                                                                                                   |                                                                                                  
tataCAGAGAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGAc                                                                                                   >  2:2948293/1‑101 (MQ=255)
 ataCAGAGAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCg                                                                                                  <  2:2947571/101‑1 (MQ=255)
    cAGAGAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGc                                                                                               <  1:2947737/101‑1 (MQ=255)
    cAGAGAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGc                                                                                               <  1:2947771/101‑1 (MQ=255)
     agagAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCa                                                                                              >  1:2948300/1‑101 (MQ=255)
     agagAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCa                                                                                              <  2:2947622/101‑1 (MQ=255)
       agaATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATa                                                                                            <  1:2947754/101‑1 (MQ=255)
        gaATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATAc                                                                                           >  1:2948305/1‑101 (MQ=255)
        gaATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATAc                                                                                           >  2:2948304/1‑101 (MQ=255)
        gaATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATAc                                                                                           <  1:2947596/101‑1 (MQ=255)
            ttAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGaa                                                                                       <  1:2947800/101‑1 (MQ=255)
            ttAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGaa                                                                                       >  2:2948311/1‑101 (MQ=255)
             tAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGaaa                                                                                      >  1:2948313/1‑101 (MQ=255)
               aCTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGaaaaa                                                                                    >  1:2948318/1‑101 (MQ=255)
                  cACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCa                                                                                 <  2:2947819/101‑1 (MQ=255)
                    cTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCaaa                                                                               <  1:2947230/101‑1 (MQ=255)
                    cTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCaaa                                                                               <  1:2947477/101‑1 (MQ=255)
                     tAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAAt                                                                              >  2:2948323/1‑101 (MQ=255)
                      aaaGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATc                                                                             >  1:2948327/1‑101 (MQ=255)
                       aaGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCa                                                                            >  2:2948330/1‑101 (MQ=255)
                          ttAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCt                                                                         >  2:2948334/1‑101 (MQ=255)
                           tAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCtt                                                                        >  1:2948338/1‑101 (MQ=255)
                              gaagaTTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTAc                                                                     <  1:2948062/101‑1 (MQ=255)
                                agaTTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACat                                                                   >  1:2948346/1‑101 (MQ=255)
                                agaTTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACat                                                                   <  1:2947825/101‑1 (MQ=255)
                                agaTTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACat                                                                   >  2:2948344/1‑101 (MQ=255)
                                 gaTTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACata                                                                  >  2:2948350/1‑101 (MQ=255)
                                 gaTTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACata                                                                  <  1:2947670/101‑1 (MQ=255)
                                   ttGAAAATTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACAtata                                                                <  2:2947358/101‑1 (MQ=255)
                                     gAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATaaa                                                              <  1:2947682/101‑1 (MQ=255)
                                      aaaaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATaaaa                                                             <  2:2947846/101‑1 (MQ=255)
                                      aaaaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATaaaa                                                             <  1:2947726/101‑1 (MQ=255)
                                      aaaaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATaaaa                                                             >  1:2948355/1‑101 (MQ=255)
                                       aaaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAg                                                            >  2:2948361/1‑101 (MQ=255)
                                       aaaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAg                                                            <  2:2947764/101‑1 (MQ=255)
                                       aaaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAg                                                            >  1:2948359/1‑101 (MQ=255)
                                       aaaGTCTTAAACATATTTTCAGAATAATCGGATTTATACGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAg                                                            <  1:2947832/101‑1 (MQ=255)
                                        aaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGc                                                           >  1:2948365/1‑101 (MQ=255)
                                        aaGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGc                                                           >  2:2948369/1‑101 (MQ=255)
                                         aGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCg                                                          <  1:2947815/101‑1 (MQ=255)
                                            cTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTgg                                                       >  1:2948376/1‑101 (MQ=255)
                                            cTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTgg                                                       >  2:2948378/1‑101 (MQ=255)
                                            cTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTgg                                                       <  2:2947916/101‑1 (MQ=255)
                                             ttAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGa                                                      >  1:2948381/1‑101 (MQ=255)
                                             ttAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGa                                                      <  1:2947873/101‑1 (MQ=255)
                                              tAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGAc                                                     >  1:2948388/1‑101 (MQ=255)
                                               aaaCATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACa                                                    >  2:2948390/1‑101 (MQ=255)
                                               aaaCATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACa                                                    <  2:2948030/101‑1 (MQ=255)
                                               aaaCATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACa                                                    <  2:2947940/101‑1 (MQ=255)
                                                aaCATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACaa                                                   >  2:2948394/1‑101 (MQ=255)
                                                 aCATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACaaa                                                  <  2:2948154/101‑1 (MQ=255)
                                                   atatTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACaaaaa                                                >  2:2948399/1‑101 (MQ=255)
                                                     atTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAAc                                              >  1:2948403/1‑101 (MQ=255)
                                                      ttttCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACa                                             >  2:2948406/1‑101 (MQ=255)
                                                      ttgtCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATaaa                                                              <  2:2947691/81‑1 (MQ=255)
                                                        ttCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGt                                           <  2:2947921/101‑1 (MQ=255)
                                                        ttCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGt                                           <  2:2947654/101‑1 (MQ=255)
                                                        ttCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGt                                           >  1:2948409/1‑101 (MQ=255)
                                                         tCAGATTAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTg                                          <  1:2947907/101‑1 (MQ=255)
                                                         tCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTg                                          <  2:2947780/101‑1 (MQ=255)
                                                          cAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGa                                         <  2:2947736/101‑1 (MQ=255)
                                                          cAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGa                                         <  1:2947714/101‑1 (MQ=255)
                                                          cAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGa                                         >  2:2948414/1‑101 (MQ=255)
                                                          cAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGa                                         <  2:2947803/101‑1 (MQ=255)
                                                            gAATAATCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAAc                                       >  1:2948417/1‑101 (MQ=255)
                                                              ataatCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACAt                                     >  2:2948420/1‑101 (MQ=255)
                                                               taatCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACAtt                                    >  1:2948425/1‑101 (MQ=255)
                                                                  tCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAAt                                 >  2:2948427/1‑101 (MQ=255)
                                                                  tCGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAAt                                 <  1:2947890/101‑1 (MQ=255)
                                                                   cGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATa                                <  1:2947936/101‑1 (MQ=255)
                                                                   cGGATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATa                                <  1:2947635/101‑1 (MQ=255)
                                                                    ggATTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACAtt                                    >  2:2948433/1‑96 (MQ=255)
                                                                      aTTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATagag                             >  2:2948438/1‑101 (MQ=255)
                                                                      aTTTATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATagag                             <  1:2947531/101‑1 (MQ=255)
                                                                        ttATATGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGAt                           >  2:2948441/1‑101 (MQ=255)
                                                                        ttATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGAt                           >  1:2948442/1‑101 (MQ=255)
                                                                         tataTGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATa                          <  2:2948003/101‑1 (MQ=255)
                                                                          atatGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATaa                         <  1:2947786/101‑1 (MQ=255)
                                                                            atGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATaaaa                       >  1:2948451/1‑101 (MQ=255)
                                                                            atGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATaaaa                       >  2:2948448/1‑101 (MQ=255)
                                                                             tGTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAAt                      <  1:2947387/101‑1 (MQ=255)
                                                                              gTTTGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAAtt                     >  2:2948455/1‑101 (MQ=255)
                                                                              gTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAAtt                     <  1:2948016/101‑1 (MQ=255)
                                                                               tttGAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTg                    >  1:2948458/1‑101 (MQ=255)
                                                                                  gAAAATTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTAc                 <  1:2948273/101‑1 (MQ=255)
                                                                                   aaaaTTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAAtt                     <  2:2947892/96‑1 (MQ=255)
                                                                                   aaaaTTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACa                >  2:2948463/1‑101 (MQ=255)
                                                                                     aaTTATTATATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTAc                 >  2:2948468/1‑98 (MQ=255)
                                                                                        tattatATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTg           >  1:2948474/1‑101 (MQ=255)
                                                                                        tattatATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTg           <  2:2947788/101‑1 (MQ=255)
                                                                                        tattatATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTg           >  2:2948473/1‑101 (MQ=255)
                                                                                          ttatATTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTa         <  2:2947905/101‑1 (MQ=255)
                                                                                           tataTTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAg        <  2:2947864/101‑1 (MQ=255)
                                                                                           tataTTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAg        >  2:2948480/1‑101 (MQ=255)
                                                                                            atatTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGa       >  1:2948483/1‑101 (MQ=255)
                                                                                             tatTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGAt      >  2:2948488/1‑101 (MQ=255)
                                                                                              atTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATa     >  1:2948492/1‑101 (MQ=255)
                                                                                              atTGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATa     <  2:2947486/101‑1 (MQ=255)
                                                                                               ttGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATAc    <  1:2947962/101‑1 (MQ=255)
                                                                                                tGGGCGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATAcc   >  1:2948497/1‑101 (MQ=255)
                                                                                                   |                                                                                                  
TATACAGAGAATTTAACTCACTAAAGTTAAGAAGATTGAAAAGTCTTAAACATATTTTCAGAATAATCGGATTTATATGTTTGAAAATTATTATATTGGACGAGCATACAGAAAAAGCAAATCACCTTTACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCG  >  NC_000913/1624618‑1624815

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