Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6920CA43.1% 101.7 / 45.7 72intergenic (+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 (31/0);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 2.69e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAA  >  NC_000913/2903542‑2903839
                                                                                                                                                        |                                                                                                                                                        
gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAc                                                                                                                                                          <  1:55586/151‑1 (MQ=255)
      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAGAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTc                                                                                                                                                    <  2:1049463/150‑1 (MQ=255)
      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACa                                                                                                                                                        >  2:315215/1‑147 (MQ=255)
        aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAAACAtt                                                                                                                                                 >  1:1070307/1‑151 (MQ=255)
          ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTGAAAACACTc                                                                                                                                                    >  2:221411/1‑146 (MQ=255)
          ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAg                                                                                                                                               <  1:599842/151‑1 (MQ=255)
            aaTATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAAAACTCAtt                                                                                                                                                 >  1:848967/1‑147 (MQ=255)
               atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAACCATTAAcccc                                                                                                                                           >  2:139306/1‑148 (MQ=255)
                  aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCAtttaactccaa                                                                                                                                        >  1:1187721/1‑141 (MQ=255)
                   aCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACAAACATTTACCTCaaaa                                                                                                                                    >  1:760163/1‑150 (MQ=255)
                       cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAg                                                                                                                               <  1:1170316/151‑1 (MQ=255)
                          tCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAAAACACTAACCTCcaaacaaaaaa                                                                                                                             >  1:1212981/1‑147 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACtt                                                                                                                       <  2:989603/151‑1 (MQ=255)
                                     aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTCTAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGa                                                                                                                  <  2:1060399/150‑1 (MQ=255)
                                        aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATtta                                                                                                             <  1:139300/151‑1 (MQ=255)
                                           aaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAACACCCCTTTATCTCAAAACAAAAAAAAAAttttcgattttttt                                                                                                           >  2:150481/1‑139 (MQ=255)
                                           aaaCACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACCCCTTATTCTCAAAAAAAAAAAAAAACTACAGAtttttt                                                                                                           >  1:579555/1‑150 (MQ=255)
                                            caGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCAATTTCCTCCAAAAAAAAAAAAAAATTCCACAttattt                                                                                                           >  1:1050447/2‑146 (MQ=255)
                                            aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTTTCCCCAAAAAAAAAAAAAAAATTCCAAAttttt                                                                                                            >  2:779913/1‑147 (MQ=255)
                                              gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCCATTAAACTCAAAAAAAAAAAAAAATTCCCGAAttttttttt                                                                                                       >  2:337586/1‑147 (MQ=255)
                                               aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                                                      <  2:372090/151‑1 (MQ=255)
                                                    aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCaa                                                                                                 <  2:635516/151‑1 (MQ=255)
                                                    aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACCTAACCa                                                                                                  <  1:31191/150‑1 (MQ=255)
                                                     cGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTATTTAAAAACCCCCAATAGCCTCAAAAAAAAAAAAAAAATCCCAATTTTTTACTTAACCaa                                                                                                 >  2:907516/1‑150 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                                                              <  2:511621/151‑1 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                                                              <  2:597128/151‑1 (MQ=255)
                                                        aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTTTTTAAAAAAACCCATAACCTCCAAAAAAAAACAAAAATTCACATTTTTTAATTTACCaa                                                                                                 >  2:1054214/1‑147 (MQ=255)
                                                           gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAAACAATTTACTCAAAAAAAAAAAAAAAATTCCGATttatttttttaacaaaaaa                                                                                             >  1:1104550/1‑134 (MQ=255)
                                                           gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTAGATTAAAAAAACTCATTATCCACCAAAAAAAAAAAATTTCCCAGATtaatatttgtac                                                                                                    >  1:652997/1‑134 (MQ=255)
                                                                     aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaa                                                                                <  2:902986/151‑1 (MQ=255)
                                                                      aTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAAATCATTTGACTCAAAACAAAAAAAAAATTCCAGATTTTTTTTTTAACAAAAGAATAAAAAAATTaaa                                                                                >  1:818732/1‑150 (MQ=255)
                                                                              ttAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAACACCCCATTTAACTCCAAACAAAAACAAACCTCCCGATTTATTATTTTATCCAAGAATAAAAAACTTAAAAAAtata                                                                         >  1:435049/1‑149 (MQ=255)
                                                                                 aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGAATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                                                                    <  2:490069/151‑1 (MQ=255)
                                                                                   gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAACCCCAAAAAAAAAACAAAATTCCAGATTTTTTTATTTACCAAAAAAAAAAAAACTTaaaaaaaaaaatttttt                                                                  >  2:891373/1‑140 (MQ=255)
                                                                                     aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAAAAACCATTATCCTCACAAAAAAAAAAAAAATCCCAAATTTTTTATTTAAAAATAAATAAAAAAATTaaaaaaaaaaatattttt                                                                 >  1:1183035/1‑138 (MQ=255)
                                                                                     aTGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTc                                                                <  1:821572/151‑1 (MQ=255)
                                                                                         aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                            <  2:919783/151‑1 (MQ=255)
                                                                                           tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTTAAAACACCCCTTATCCCCTAAAAAAAAAAAAAACTTCAACTTTTTTTTTTTACCAAAGAATAAAAAAATTAAAAAAAAAAATGTtttt                                                                 >  2:1025060/1‑144 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                                                           <  1:273870/149‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                                                           <  1:295981/149‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGa                                                          <  2:879707/150‑1 (MQ=255)
                                                                                             gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATATTTAAAAAAACTCCCTATCCTCCAAAAAAAAAAAAAAATTCCGATTTTTTATTTTAACAAAGAAAAAAAAAATTAAAAAAAATAATGTtttttttaaaaa                                                         >  2:601923/1‑142 (MQ=255)
                                                                                              atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAACCTCAAAAAAAAAACAAAAATCCAAATTTTTTATTTTAACAATAAAAAAAAAAATTAAAAAAAAAACTTTTATTTTTAAGaaa                                                        >  2:524703/1‑150 (MQ=255)
                                                                                                atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAAACACATTTCCTTCAAAAAAAAACAAAAATCCACAGTTTTTTTTTTACCAAAAAATAAAAAACTTAAAAAAAAAAATGTTATTTTTaaaaaaaa                                                      >  1:1197616/1‑150 (MQ=255)
                                                                                                atCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaa                                                      <  2:1117644/150‑1 (MQ=255)
                                                                                                   tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTTAAAAACCCATTATACCCACAAAAAAAAAAAAATTCCCGATTTTTTACATTAACAAAAAATAACAAACCTAAAAAAAAAAATTTTTTTTTTAaaaaaaaa                                                     >  1:920812/1‑148 (MQ=255)
                                                                                                    ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCCTTATCCCCCAAAAAAAAAAAAAAATCCAGATTTTTTTATTTAACAAAGAAAAAAAAACTTAAAAAAAAAAATTTTTTTTTAAAAAAAAAAATa                                                  >  2:738080/1‑150 (MQ=255)
                                                                                                       ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAAc                                                <  2:124028/149‑1 (MQ=255)
                                                                                                        tCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAg                                              <  1:961060/150‑1 (MQ=255)
                                                                                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGATATTTAAAAACACACATTTTTCCCAAAAAAAAAAAAAACAACCTTATTTTTTAATTATCCAAAGAAAAAAAAAACTAAAAAAAATTATTTTATTTTTAaaaaaaaaaaaaaaagca                                           >  1:831288/1‑139 (MQ=255)
                                                                                                             aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTGAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTa                                        <  1:533765/151‑1 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAACTCCATTATCTCAAAAAAAAAAAAAAAATTCAGATTTTTTTTTTTAACAAATAAAAAAAAAAATAAAAAAAAAAATGTTTTTTTTAAGAAAAAAATaaa                                                >  1:918183/1‑141 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTaa                                       <  1:81907/151‑1 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAATCCCTATCCTCCAAAAAAAAACAAAAATCCAAATTTTTTTATTTAACAAATAAAAAAAAAATTAAAAAAAAAAATGTTTTTTTTAAAAAAAAAATAAAAAGAATTaa                                       >  2:646684/1‑150 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCCATTAGCCCCCAAAAAAAAAAAAAAATCCCAATTTTTTTCTTTAACAAAGAAAAAAAAACTTAAAAAAAATAATTTTTTTTTTAAAAAAAAAATAAAAAGGATTaa                                       >  2:548083/1‑150 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAAt                                      <  2:238250/151‑1 (MQ=255)
                                                                                                                    ttCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCCTTAAACTCCAAAAAAAAAAAAAAATCCCCATTTTTTTCTTTAACAAAGAATAAAAAAATTTAAAAAAATACTGTTTTTTTAAAAAAAAAAATAAAAAGCATTaa                                       >  2:1049/1‑145 (MQ=255)
                                                                                                                      catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                                <  1:192730/150‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaa                            <  2:638627/151‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                             <  1:1037498/150‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                             <  1:1483/150‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                             <  2:859331/150‑1 (MQ=255)
                                                                                                                             aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                        <  1:253421/151‑1 (MQ=255)
                                                                                                                             aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                        <  1:987429/151‑1 (MQ=255)
                                                                                                                                gTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACATCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaata                     <  1:868093/151‑1 (MQ=255)
                                                                                                                                  ataGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataa                    <  1:928932/150‑1 (MQ=255)
                                                                                                                                          ggTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAGTAAAAAAATAATAGATTaaa           <  1:614657/151‑1 (MQ=255)
                                                                                                                                           gTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaaa           >  2:875738/1‑150 (MQ=255)
                                                                                                                                              tATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAtt        >  2:175054/1‑151 (MQ=255)
                                                                                                                                                   tAAAACCACTCATTAGCCTCAAAACAAAGACAAAAATCCAGAGTTATTACTTAACCAATGAATAACAAGATTAAAAAATAAACTGTAATTATAAAGAAAAAATTAACAGGCAGTAATTCAATATTAAAAAAATAATAGATTAAAATTTCt     <  1:210914/150‑1 (MQ=255)
                                                                                                                                                    aaaaCCACTCATTAGCCTCAAAACAAAGATAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt    <  2:667626/150‑1 (MQ=255)
                                                                                                                                                     aaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa  <  2:877801/151‑1 (MQ=255)
                                                                                                                                                     aaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTaa  <  2:1021269/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                        
GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAA  >  NC_000913/2903542‑2903839

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