Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,693 C→T 6.3% intergenic (‑307/+20) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6930CT 6.3% 203.3 / 3.7 63intergenic (‑307/+20)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (41/18);  new base T (2/2);  total (43/20)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

TTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACA  >  NC_000913/1764548‑1764837
                                                                                                                                                  |                                                                                                                                                
ttAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAGAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAGGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGac                                                                                                                                               >  1:663017/1‑149 (MQ=255)
    gTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAg                                                                                                                                           <  2:509312/149‑1 (MQ=255)
     tATCTCATAGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAgg                                                                                                                                          <  2:663017/149‑1 (MQ=255)
        ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCtt                                                                                                                                       >  4:649561/1‑149 (MQ=255)
         tcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTa                                                                                                                                      >  1:586230/1‑149 (MQ=255)
          cATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTaa                                                                                                                                     <  4:504424/149‑1 (MQ=255)
           aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                                    >  1:129964/1‑149 (MQ=255)
                        ggTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCa                                                                                                                       >  2:739771/1‑149 (MQ=255)
                               atCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                                                                                >  2:117149/1‑149 (MQ=255)
                                   aaGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTa                                                                                                            >  3:645513/1‑149 (MQ=255)
                                      gAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                     >  4:730755/1‑121 (MQ=39)
                                      gAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                     <  3:730755/121‑1 (MQ=39)
                                         cATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCg                                                                                                      >  1:68474/1‑149 (MQ=255)
                                          aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                                                                     >  3:589555/1‑149 (MQ=255)
                                          aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                                                                     >  1:74229/1‑149 (MQ=255)
                                          aTTAAGGCTGTGATGAAAAAGAAAAAGTTATATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                                                                     >  4:614956/1‑149 (MQ=255)
                                              aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAcg                                                                                                 >  3:248110/1‑148 (MQ=255)
                                                gCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCActct                                                                                               <  4:248110/149‑1 (MQ=255)
                                                        gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGc                                                                                                                                                  <  2:394865/90‑1 (MQ=255)
                                                        gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGc                                                                                                                                                  >  1:394865/1‑90 (MQ=255)
                                                        gaaaaagaaaaaGTTATATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                       <  3:614956/149‑1 (MQ=255)
                                                           aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                                                                     >  2:418334/1‑148 (MQ=255)
                                                           aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                                                                     <  1:418334/148‑1 (MQ=255)
                                                           aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAg                                                                                    >  4:123008/1‑149 (MQ=255)
                                                           aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAg                                                                                    >  2:544683/1‑149 (MQ=255)
                                                              gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                   >  2:310452/1‑147 (MQ=255)
                                                              gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                   <  1:310452/147‑1 (MQ=255)
                                                                aaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCtt                                                                               >  4:311137/1‑149 (MQ=255)
                                                                    gTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                           <  3:311137/149‑1 (MQ=255)
                                                                     ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                          >  1:801156/1‑149 (MQ=255)
                                                                        tATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                                                       <  3:123008/149‑1 (MQ=255)
                                                                         aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCg                                                                      >  2:241914/1‑149 (MQ=255)
                                                                          ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATc                                                                                                                     >  2:311570/1‑101 (MQ=255)
                                                                          ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATc                                                                                                                     <  1:311570/101‑1 (MQ=255)
                                                                          ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                     >  4:323100/1‑149 (MQ=255)
                                                                             aaaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa                                                                  <  2:68474/149‑1 (MQ=255)
                                                                             aaaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCGGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa                                                                  <  3:649561/149‑1 (MQ=255)
                                                                                aTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa                                                               >  3:351550/1‑149 (MQ=255)
                                                                                        ttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa                                                       <  4:645513/149‑1 (MQ=255)
                                                                                            gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaaa                                                   >  4:62901/1‑149 (MQ=255)
                                                                                               aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                >  3:269748/1‑149 (MQ=255)
                                                                                                 aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTaa                                              >  3:660377/1‑149 (MQ=255)
                                                                                                      aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgc                                         >  3:761955/1‑149 (MQ=255)
                                                                                                              cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                 <  4:351550/149‑1 (MQ=255)
                                                                                                                          aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGa                     >  3:23044/1‑149 (MQ=255)
                                                                                                                           tAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAt                    <  1:739771/149‑1 (MQ=255)
                                                                                                                              gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgggg                                                                        >  2:593632/1‑94 (MQ=255)
                                                                                                                              gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgggg                                                                        <  1:593632/94‑1 (MQ=255)
                                                                                                                              gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                 >  1:35109/1‑149 (MQ=255)
                                                                                                                              gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                 >  2:68551/1‑149 (MQ=255)
                                                                                                                               aaaaCCCTTAAGTCTGTGTGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                     <  3:281747/96‑1 (MQ=255)
                                                                                                                               aaaaCCCTTAAGTCTGTGTGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                     >  4:281747/1‑96 (MQ=255)
                                                                                                                                    ccTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat           >  1:355602/1‑149 (MQ=255)
                                                                                                                                    ccTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat           >  2:736748/1‑149 (MQ=255)
                                                                                                                                    ccTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat           >  3:755114/1‑149 (MQ=255)
                                                                                                                                    ccTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat           >  2:616216/1‑149 (MQ=255)
                                                                                                                                     cTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTattatt          >  3:742265/1‑149 (MQ=255)
                                                                                                                                     cTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTattatt          >  4:277271/1‑149 (MQ=255)
                                                                                                                                      ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt         <  4:589555/149‑1 (MQ=255)
                                                                                                                                       tAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt        >  2:793514/1‑149 (MQ=255)
                                                                                                                                         aGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg      >  3:300686/1‑149 (MQ=255)
                                                                                                                                           tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa    <  1:117149/149‑1 (MQ=255)
                                                                                                                                           tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa    >  1:152652/1‑149 (MQ=255)
                                                                                                                                           tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa    >  4:803540/1‑149 (MQ=255)
                                                                                                                                             tgtgCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa  >  1:254202/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACA  >  NC_000913/1764548‑1764837

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