Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6920CA42.2% 117.6 / 48.4 71intergenic (+318/+1055)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (3/38);  new base A (30/0);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 5.69e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903554‑2903841
                                                                                                                                           |                                                                                                                                                        
aaTATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGcc                                                                                                                                             <  2:149260/151‑1 (MQ=255)
        cTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACCCATTAACCTCaaaa                                                                                                                                      >  1:674851/1‑149 (MQ=255)
         tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAAAACCTTAACCTCaaaaaa                                                                                                                                    >  1:637913/1‑148 (MQ=255)
          acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAACACATTAACCTCCaaaaaa                                                                                                                                   >  1:412791/1‑150 (MQ=255)
          acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCCATCATTAACCTCAaaaaaa                                                                                                                                   >  1:643337/1‑150 (MQ=255)
          acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACaaa                                                                                                                                  <  1:527013/151‑1 (MQ=255)
          acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACaaa                                                                                                                                  <  1:24125/151‑1 (MQ=255)
          acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACaa                                                                                                                                   <  1:982160/150‑1 (MQ=255)
          acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACa                                                                                                                                    <  1:668419/149‑1 (MQ=255)
            acTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACACATTAGCCTCCAaaaaaa                                                                                                                                  >  2:495934/1‑149 (MQ=255)
              tCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCATTTTCCTCCAaaaaaaaa                                                                                                                                >  1:454229/1‑147 (MQ=255)
               cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaa                                                                                                                             <  1:593028/151‑1 (MQ=255)
                 aTGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAAACATTAGCCTCaaaaaaaaaaaaaa                                                                                                                            >  1:408605/1‑150 (MQ=255)
                   gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCCCTAGCCTCaaaaaaaaaaaaaa                                                                                                                            >  2:405636/1‑148 (MQ=255)
                   gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACt                                                                                                                          <  1:349591/150‑1 (MQ=255)
                      tttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAATCATTTACCTCCAAAAAAAAAAAAAAATTcc                                                                                                                       >  1:717584/1‑151 (MQ=255)
                       ttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTcc                                                                                                                       <  1:594845/149‑1 (MQ=255)
                              aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCCTTTAACTCCAAAAAAAAAAAAAATTCCAAAtttttt                                                                                                             >  2:28984/1‑151 (MQ=255)
                              aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATttat                                                                                                              <  2:497990/150‑1 (MQ=255)
                                  gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAACTCATTAACCTCAAAAAAAAAAAAAAACTCCAGATTTTTTTCtt                                                                                                         >  2:667356/1‑151 (MQ=255)
                                   aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                                                        <  1:154833/151‑1 (MQ=255)
                                      aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAAcc                                                                                                     <  2:920945/151‑1 (MQ=255)
                                          gAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTATCCTCAAAAAAAAAAAAAAAATCCAGATTTATTTCTTAATCaa                                                                                                   >  1:490081/1‑149 (MQ=255)
                                           aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTAGCCTCAAAACAAAAACAAAAATCCAGATTTATTACTTTAACaa                                                                                                   >  2:1018674/1‑148 (MQ=255)
                                            aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGaa                                                                                               <  2:1001848/151‑1 (MQ=255)
                                              cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAAt                                                                                              <  2:698981/150‑1 (MQ=255)
                                                cACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAc                                                                                           <  2:87975/151‑1 (MQ=255)
                                                 aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAACTCTTTATTCTCAAAACAAAAACAAAATTTCCGATTTTTTTCTTTAACaaataa                                                                                               >  1:428989/1‑142 (MQ=255)
                                                   tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACAACCCTTTAACTCAAAAAAAAGAAAAAAATCCACATTTTTTTATTaaaccaaaaaaaaaaa                                                                                         >  1:563635/1‑136 (MQ=255)
                                                    ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACTCATTTTCCTCAAAAAAAAAAAAAAATTTCAGATTTTTTTTTTTAACAAAGAATAaaaa                                                                                         >  2:102239/1‑149 (MQ=255)
                                                    ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                         <  1:571245/149‑1 (MQ=255)
                                                         aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaa                                                                                  <  2:403643/151‑1 (MQ=255)
                                                          aTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAAAACAATAACCTCAAAAAAAAAAAAAAATTTCAGAATTATTAATTTAACAAAGAAAAAAAAACCTaaa                                                                                  >  1:483104/1‑150 (MQ=255)
                                                           ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTTGCCTCAAAACAAAAACAAAATTCCCGATTTATTTATTAAACAAAAAAAAAAAAACTTtaaa                                                                                 >  2:791612/1‑150 (MQ=255)
                                                               cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                                                            <  1:108742/151‑1 (MQ=255)
                                                               cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAata                                                                             <  2:267244/150‑1 (MQ=255)
                                                                aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAAtata                                                                           <  1:301800/151‑1 (MQ=255)
                                                                  ttAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTATAAAAACACCCATTTTCCCCTAAACAAAAAAAAAATTCTCGATTTTTTTTTTTTACAAAGAAAAAAAAACTTTAAAAAtataa                                                                          >  1:933322/1‑149 (MQ=255)
                                                                       gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAATATACTGTTAtt                                                                   <  2:873743/151‑1 (MQ=255)
                                                                        tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTTTAAAAAAAAACCATTAACCCCAAAAAAAAAAAAAACATCCATATTTTTTACTTAAACAAAAAATAAAAAAACTAAAAAAAATAATTTTAtt                                                                   >  2:333792/1‑151 (MQ=255)
                                                                           gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTATACCCCAAAAAAAAAAAAAAATTCCGATTTATTTATTTAACAAAAAAAAAAAAAATTAAAAAAaataaatttttttt                                                                 >  1:47888/1‑140 (MQ=255)
                                                                           gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCt                                                                 <  1:702384/150‑1 (MQ=255)
                                                                            aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTAAGCTCAAAAAAAAAAAAAAAATCCCAATTTTTTTTTTAAACAATGAAAAAAAAAATTTAAAAAAATAATTTTTTTCt                                                                 >  1:394896/1‑149 (MQ=255)
                                                                            aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAACCCTTAGGCCCAAAAAAAAAAAAAAAATCCAGATTTTTTTTTTAAACCAAGAAAAAAAAACTTTAAAAAaataatttttt                                                                    >  2:410868/1‑139 (MQ=255)
                                                                             aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATATGTAAAAACACTCAATAGCCTCAAAACAAAAAAAAAATTTCCGATTTATTACTTTACCAATGAATAACAAACTTTAAAAAAATAATTTTATTCtta                                                               >  2:191544/1‑148 (MQ=255)
                                                                               tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAGCCTCCAAAAAAAAAAAAACTTCCAGATTTTTTTTTTTAACAAAGAAAAAAAAACTTAAAAAAAATACTGTTATTTTaaa                                                              >  2:971044/1‑149 (MQ=255)
                                                                                  atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaa                                                         <  1:600198/151‑1 (MQ=255)
                                                                                     tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTAGCCTCAAAAAAAAAAAAAAATTTCAGATTTATTTATTTAACAAATAAAAAAAAAATTAAAAAAAATAATGTTATTTTTAAGaaaaaa                                                       >  1:409654/1‑150 (MQ=255)
                                                                                        ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAt                                                      <  2:1019940/148‑1 (MQ=255)
                                                                                        ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAt                                                      <  1:783207/148‑1 (MQ=255)
                                                                                          aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAAc                                                  <  2:537364/150‑1 (MQ=255)
                                                                                              aTAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTTGCCTCCAAAAAAAAAAAAAAATTCAGAATTTTTTCTTTAACAAAGAAAAAAAAAATTTAAAAAAATAATTTTTTTTTAAAGAAAAAAATAAAAgg                                               >  2:233298/1‑149 (MQ=255)
                                                                                                aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAt                                            <  2:551619/150‑1 (MQ=255)
                                                                                                 aGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACAACCCTTTGCCTCAAAACAAAAAAAAAATTCCAGATTTTTTATTTTACCAATAAATAAAAAGCTTAAAAAAAATAATGGTATTTTTAAAAAAAAAATAAAAAGCAtt                                           >  2:971034/1‑150 (MQ=255)
                                                                                                 aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTa                                          <  1:204801/151‑1 (MQ=255)
                                                                                                   aaCTATTCATCAACAGCGTATAGAGGCGGCTCTATGAAAACAATCTTTTACCTCAAAACAAAAATAAACTTCCAGACTTTTTTTTTAACCCATTAAAAACAAGCCTTAAAAAAATAAAGTTTTTCTTAAGAAAAAAATAAAAAGCAtt                                           >  2:931169/1‑148 (MQ=255)
                                                                                                      tATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCCTTATCCTCCAAAAAAAAAAAAAATTTCCGATTTTTTTATTAAACAAAAAATAAAAAACCTTAAAAAAATAATGTTTTTTTTAAAAAAAAAATAAAAAGCATTaaat                                       >  1:861614/1‑147 (MQ=255)
                                                                                                         tcatcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAat                                   <  1:325601/150‑1 (MQ=255)
                                                                                                            tcaACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTAGACTCCAAACAAAAAAAAAATTTCAGATTTTTTAATTTAACAAAGAATAACAAACTTAAAAAATATACTGTTTTTTTTAAAAAAAAAATAAAAAGCATTTATTTAATAtt                                >  2:268062/1‑150 (MQ=255)
                                                                                                             caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                               <  1:74351/150‑1 (MQ=255)
                                                                                                              aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaa                             <  1:207148/151‑1 (MQ=255)
                                                                                                                        agagGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAataata                    <  2:846708/150‑1 (MQ=255)
                                                                                                                         gcggCGGTTATGTAAAAACACTCCGTAGCCTCAAAAAAAAGACAAACTTCCGGATTTATTACTTAACCAATTAATAAAAAGCTTAAAAAAAATACTGTTATTCTAAAGAAAAAAATAAAAAGCATTAATTCAATATTAAAAAaataa                      >  1:595519/3‑147 (MQ=255)
                                                                                                                          agGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGa                  <  2:926480/150‑1 (MQ=255)
                                                                                                                          agGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAt                 <  1:489068/151‑1 (MQ=255)
                                                                                                                               gTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaaaa            <  2:762890/151‑1 (MQ=255)
                                                                                                                                 tATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATAAATTAAATATTAAAAAAATAATAGATTAAAAtt          >  1:338224/1‑151 (MQ=255)
                                                                                                                                     gTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTc        <  1:444482/150‑1 (MQ=255)
                                                                                                                                     gTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt       <  2:832517/151‑1 (MQ=255)
                                                                                                                                        aaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa    <  2:1004683/151‑1 (MQ=255)
                                                                                                                                          aCCACTCATTAGCCTCAAAACAAAGACAAACCTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACg  <  2:461577/151‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                        
AATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903554‑2903841

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