Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6860GT42.6% 98.4 / 31.5 70intergenic (+312/+1061)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (8/31);  new base T (29/0);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 4.34e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAT  >  NC_000913/2903538‑2903832
                                                                                                                                                    |                                                                                                                                                     
cccTGTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTa                                                                                                                                                     <  2:432815/151‑1 (MQ=255)
    gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTATaaaa                                                                                                                                                  >  2:157747/1‑150 (MQ=255)
    gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTaaaa                                                                                                                                                  <  1:46261/150‑1 (MQ=255)
     tttttACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTaaaa                                                                                                                                                  <  1:902995/149‑1 (MQ=255)
           cAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACTc                                                                                                                                            >  1:551400/1‑149 (MQ=255)
            aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTTTTAAAAACaaatatt                                                                                                                                         >  2:350291/1‑145 (MQ=255)
             aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTa                                                                                                                                        <  1:288539/151‑1 (MQ=255)
              ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAg                                                                                                                                      <  2:176396/151‑1 (MQ=255)
               gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAacacccctt                                                                                                                                         >  1:950038/1‑143 (MQ=255)
                  tatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTaaaaaaaccccttatcc                                                                                                                                     >  1:192785/1‑136 (MQ=255)
                   atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTATcccc                                                                                                                                  >  1:659360/1‑148 (MQ=255)
                     tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                                                                                                                                <  1:564263/150‑1 (MQ=255)
                     tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaa                                                                                                                                <  1:740673/150‑1 (MQ=255)
                         tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTccatatcccccaaaa                                                                                                                              >  2:754706/1‑135 (MQ=255)
                         tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACACATTATCCTCaaaaaa                                                                                                                            >  2:78833/1‑148 (MQ=255)
                               cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGAc                                                                                                                       <  1:662511/149‑1 (MQ=255)
                                 aTGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTTACCTCCAaacaaaaacaaa                                                                                                                   >  1:299018/1‑150 (MQ=255)
                                   gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACACATTAGACTCAAAAAAAAAACaaa                                                                                                                   >  1:790558/1‑148 (MQ=255)
                                   gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAAAACTAATTTACCTTaaaaaaaaaaaaaaa                                                                                                                  >  2:489044/1‑148 (MQ=255)
                                      tttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAAACCTTAAACTCAAAAAAAAAAAAAACTTcc                                                                                                              >  2:757865/1‑150 (MQ=255)
                                         aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTAAGCTCAAAAAAAAAACaaaaattc                                                                                                              >  1:12247/1‑142 (MQ=255)
                                             tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCATTATTCTTCAAACAAAAAAAAAcattcaaaatt                                                                                                        >  2:297924/1‑139 (MQ=255)
                                              aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACCCCTCATTTACCTCAAAACAaaaaaaaaaatccccatttttt                                                                                                     >  1:248400/1‑137 (MQ=255)
                                                  gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACtt                                                                                                <  1:613415/151‑1 (MQ=255)
                                                    cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTaa                                                                                              <  1:592690/151‑1 (MQ=255)
                                                      aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAAcc                                                                                            <  1:465817/151‑1 (MQ=255)
                                                       tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAACACATCATTATTCCCAAAAAAAAAAAAAAACTTCCGATTTTTTACttta                                                                                              >  2:42566/1‑146 (MQ=255)
                                                       tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAAAACTCACTATCCTCCAAAAAAAAACAAAAATCCAGAATTATTATTTaaac                                                                                            >  1:665242/1‑148 (MQ=255)
                                                         cGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAAt                                                                                         <  2:187620/151‑1 (MQ=255)
                                                          gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACACATTAGCCTCAAAAAAAAAACAAAAATCCCAATTTATTTCTTAAc                                                                                             >  1:609542/1‑146 (MQ=255)
                                                                cACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTAACCTCAAAAAAAAAAAAAAACTTCAGATTTTTTTCTTTAACaaaaaaaaa                                                                                   >  1:287336/1‑143 (MQ=255)
                                                                 aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACAATCATTTTACTCAAAAAAAAAACAAACTTTCAGATTTATTTTTTTAACAAAGAATaa                                                                                   >  2:52550/1‑149 (MQ=255)
                                                                 aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACa                                                                                 <  2:432861/151‑1 (MQ=255)
                                                                  cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACACAATTCCCTCACAAAAAAAAAAAAATTCCAGATTTTTTATTTTAACaaataaaaaaa                                                                                 >  2:441074/1‑141 (MQ=255)
                                                                  cTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACCCATTTTCCTCAAAAAAAAAAAAAAAATCCCGATTTTTTTCTTTAACaaaaaaaaaaa                                                                                 >  1:637960/1‑141 (MQ=255)
                                                                     tATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGc                                                                              <  2:665075/150‑1 (MQ=255)
                                                                      aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTTGCCTCAAAAAAAAAAAAAAAATCCCAATTTTTTTCTTTACCAAAAAATAACaa                                                                                >  1:1130429/1‑147 (MQ=255)
                                                                           ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTTTCCCCAAAAAAAAAACAAAATTTCCAATTTTTTTTTTTAACAATGAATAaaaaaaataaaa                                                                        >  2:946702/1‑142 (MQ=255)
                                                                              tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTAGTCATCAACAGCGTAGAGAGGCGGTTATGTAAAACCACTCATTAGCGTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAat                                                                     <  2:510375/150‑1 (MQ=255)
                                                                               cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACACATTTACCTCAAAAAAAAAAAAAAAATCCAGATTTTTTTATTTACCAAAGAATAAAAAAATTaaaaaaaa                                                                    >  2:555877/1‑148 (MQ=255)
                                                                                aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTATAAAAACAATCAATATCCCCCAAAAAAAAAACAAATTCTTGATTTTTTTTTTAACCAATGAAAAAAAAAACTaaaaaaaa                                                                    >  1:140680/1‑147 (MQ=255)
                                                                                aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAACCTCAAAAAAAAAAAAAAACTTCAGATTTTTTACTTTACCAATAAATAACAAACTTaaaaaaaa                                                                    >  1:918033/1‑147 (MQ=255)
                                                                                    aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTTTAAAAACACCCCTTATCCTCCAAACAAAAACAAAACTTCAAGTTTTTTTTTTTAACCAAGAAAAAAAAACCTAaaaaaaaaaatgt                                                               >  1:503950/1‑144 (MQ=255)
                                                                                      aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTGTTTAAAACCACCCATCAGTCTCTAAAAAAAAACCAACCTCCAGCTTTTTTACTTAACAAATAAAAAAAAAGAATAAAAAATATACttttt                                                            >  1:54770/1‑150 (MQ=255)
                                                                                      aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTATGTTAAAAAAAACATTTTACCCAAAAAAAAAAAAAAATTTCCGATTTATTATTTTAACAAAGAATAAAAAAATTAAAAAAAATAATGtt                                                             >  2:116137/1‑150 (MQ=255)
                                                                                       gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAt                                                           <  1:620306/151‑1 (MQ=255)
                                                                                        tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACCCCCATTATCCTCAAAAAAAAAAAAAACATCCCGAATTTTTATTTTAACAAAAAATAAAAAACTTAAAAAAAATACTGTtttt                                                          >  1:974812/1‑151 (MQ=255)
                                                                                          tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCATTAACCTCAAAAAAAAAAAAAAAATCCAAATTTTTTAATTTACCAATGAATAAAAAACTTAAAAAAAATAATTTTAtttt                                                        >  1:3447/1‑149 (MQ=255)
                                                                                             aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                     <  2:90732/151‑1 (MQ=255)
                                                                                                ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                                                    <  1:908765/149‑1 (MQ=255)
                                                                                                  atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaa                                                <  2:745405/151‑1 (MQ=255)
                                                                                                       tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCATTATTCTCAAAAAAAAAAAAAAAATCCAGAATTTTTACTTAAACAAAAAAAAAAAAAATTAAAAAAAAAACTTTTTTTTTTAaaaaaaaa                                              >  2:74794/1‑148 (MQ=255)
                                                                                                              aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCa                                    <  1:245009/151‑1 (MQ=255)
                                                                                                              aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCa                                    <  1:158984/151‑1 (MQ=255)
                                                                                                               tAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTTACCTCCAAAAAAAAACAAAAATTCAGATTTTTTTTTTTAACAAATAAAAAAAAACTTTAAAAAAAAAATTTTTTTTTTaaaaaaaaaa                                              >  1:157028/1‑141 (MQ=255)
                                                                                                                aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAt                                   <  2:1052371/150‑1 (MQ=255)
                                                                                                                aaGAACTATTCATCAACAGCGTATAGAGGCGGGTTTGTAAAAACACTCATTAGCCCCAAAAAAAAAAAAAAAATTCAGATTTATTTATTTAACAATGAAAAAAAAAAATAAAAAAATAAATTTTTTTTTAAAAAAAAAATAAAAAgcatt                                   >  1:677540/1‑146 (MQ=255)
                                                                                                                   aaCTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACCCCCCTTATACTCCAAAAAAAAAAAAAATTCCCAATTTTTTACTTAAACAAAAAATAAAAAAATTTAAAAAAAAATTGTTTTTTTAGAGAAAAAAATAAAAAGCATTaa                                >  1:972823/1‑150 (MQ=255)
                                                                                                                       aTTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAACCATTTACCTCAAAAAAAAAAAAAAACTTCAGATTTATTTCTTAACCAAATAAAAAAAAAATTAAAAAAAAAAATGTTTTTTTAAAAAAAAAAATAAAAAGCAATAAtt                              >  1:624543/1‑148 (MQ=255)
                                                                                                                          catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                        <  2:949797/151‑1 (MQ=255)
                                                                                                                          catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                        <  1:757165/151‑1 (MQ=255)
                                                                                                                               aCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaa                   <  1:314880/151‑1 (MQ=255)
                                                                                                                                 aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                 <  1:1129747/151‑1 (MQ=255)
                                                                                                                                      ataGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataa             <  2:30945/150‑1 (MQ=255)
                                                                                                                                       taGAGGCGGTTATTTAAAAACACTCATTAGCCCCAAAAAAAAAAAAAAAATCCAGAATTATTACTTAAACAAAGAATAAAAAAAATAAAAAAAATAATGTTATTTTTAAGAAAAAAATAAAAAGAATTAATTCAATATTAAAAAaataa             >  2:502909/1‑149 (MQ=255)
                                                                                                                                         gagGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa         <  2:47406/151‑1 (MQ=255)
                                                                                                                                          agGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa         <  2:47398/150‑1 (MQ=255)
                                                                                                                                            gCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa      <  1:57990/151‑1 (MQ=255)
                                                                                                                                             cGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa      >  1:1109039/1‑150 (MQ=255)
                                                                                                                                                 tATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAt  <  1:576527/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCTGTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAT  >  NC_000913/2903538‑2903832

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