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

GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903542‑2903841
                                                                                                                                                        |                                                                                                                                                          
gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAc                                                                                                                                                            <  1:9449/151‑1 (MQ=255)
              tatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACCCATTAAcc                                                                                                                                               >  2:282599/1‑149 (MQ=255)
                     tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACTCATTAACCTCaaaaaa                                                                                                                                      >  2:612232/1‑148 (MQ=255)
                          tCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTAACCTCAAAACaaaaa                                                                                                                                  >  2:795674/1‑147 (MQ=255)
                           cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACa                                                                                                                                <  1:90445/150‑1 (MQ=255)
                              tGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGATTTTTAAAAACACACATTATTCCCCACACAAAAACaaaa                                                                                                                             >  1:877110/1‑149 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACt                                                                                                                            <  2:519006/150‑1 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACt                                                                                                                            <  1:66907/150‑1 (MQ=255)
                                   ttAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTc                                                                                                                          <  2:127303/148‑1 (MQ=255)
                                     aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACACATTAAACTCaaaaaaaaaaaaaaaattcaaa                                                                                                                      >  1:725165/1‑142 (MQ=255)
                                     aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGAt                                                                                                                   <  1:735981/151‑1 (MQ=255)
                                      gAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGAtt                                                                                                                  <  2:103100/151‑1 (MQ=255)
                                       aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCACTATGCTCCAAAAAAAAAAAAACATTCCGAttt                                                                                                                 >  2:427793/1‑151 (MQ=255)
                                         tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACACCTTAAACTCAAAAAAAAAAAAAAAATCCCGAttt                                                                                                                 >  1:564201/1‑149 (MQ=255)
                                              gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAAAAACTCATTAACTCAAAACAAAAAAAAAccttcaaacttt                                                                                                                  >  1:674713/1‑134 (MQ=255)
                                              gACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTAACCTCAAAAAAAAAAAAAAATTCCACATttatt                                                                                                              >  2:851488/1‑147 (MQ=255)
                                                 aaTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTTTCCTCTAAAAAAAAAAAAAACTTcctagtttttaattaa                                                                                                           >  1:341404/1‑136 (MQ=255)
                                                 aaTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAAc                                                                                                       <  2:674335/151‑1 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAAAACATTATCCTCAAAAAAAAAAAAAAAATTCAAATTTATTTCTTAAACa                                                                                                     >  1:387590/1‑147 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                                                                 <  2:373648/151‑1 (MQ=255)
                                                         tCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTATCCTCAAAAAAAAAACAAAAATTCAGATTTTTTTATTAAACAATGaa                                                                                                >  2:431305/1‑150 (MQ=255)
                                                          cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACCCCTTAACCCCAAAAAAAAAAAAAAAATCCCGAATTTTTTATTTAACCaaaaa                                                                                                 >  2:821272/1‑146 (MQ=255)
                                                                    cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTa                                                                                     <  2:698205/150‑1 (MQ=255)
                                                                      aTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaaa                                                                                  <  1:478618/151‑1 (MQ=255)
                                                                        tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACACACATTTACCTCAAAACAAAAAAAAAAATCCAAAATTTTTTATTTAACAAATAaaaaaaaaaataaaaa                                                                                  >  1:488503/1‑142 (MQ=255)
                                                                           cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                                                             <  2:903321/151‑1 (MQ=255)
                                                                               tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACTCATTAACCTCAAAAAAAAAAAAAAAATCCCAATTTTTTTCTTTAACAAAAAATAAAAAAATTaaaaaaaaaaat                                                                          >  1:885437/1‑144 (MQ=255)
                                                                                 aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCCATTAATCTCAAAAAAAAAAAAAAACTCCCGATTTTTTATTTTACCAATGAAAAAAAAACTTtaaaaaaaaaatttt                                                                       >  2:207427/1‑142 (MQ=255)
                                                                                  aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCATTATCCTCAAAAAAAAAAAAAAAATCCAGATTTTTTTCTTTACCAAAGAAAAAAAAACTTAAAAAAAATACTGtt                                                                      >  2:468960/1‑150 (MQ=255)
                                                                                    tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCCTTAGACTCAAAAAAAAAACAAAATTCCAGATTTTTTTCTTTACCCATGAAAAAAAAAAATAAAAAAAAAAATGTTAtt                                                                   >  1:494512/1‑151 (MQ=255)
                                                                                      tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGCTTTATGAAAACACTCCTTCATGTTCTAAAAAATGAAAAAGTTCCAACTTTTTTAATTAACTACTGATGACAAGGTTAAAAAAAATACTattatt                                                                   >  2:425917/1‑149 (MQ=255)
                                                                                       gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTa                                                                <  1:22406/151‑1 (MQ=255)
                                                                                        aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaa                                                               <  1:112988/151‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACACATTTTCCTCCAAAAAAAAAAAAAAATCCACATTTTTTTTTTTAACAAAAAAAAAAAAACTTAAAAAAAATAATGTTATTTTaaaaa                                                            >  2:732619/1‑148 (MQ=255)
                                                                                             gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAAAACATTTTCCTCAAAAAAAAAAAAAACTTTCCAATTTTTTTCTTTAACAAAAAAAAAAAAACTTAAAAAAAATAATGTtttttttaaaaa                                                           >  1:400188/1‑142 (MQ=255)
                                                                                             gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaa                                                          <  1:584951/151‑1 (MQ=255)
                                                                                              atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAAAACATTTTACTCAAAAAAAAAAAAAAAATTCCGATTTTTTTTTTAAACAAAAAAAAAAAAAAATAAAAAAAAAACTGTTTTTTTAaaaaaa                                                          >  1:845219/1‑150 (MQ=255)
                                                                                              atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaa                                                         <  1:427782/151‑1 (MQ=255)
                                                                                                    ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTa                                                    <  1:178917/150‑1 (MQ=255)
                                                                                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTAGCCCCCAAAAAAAAACAAACTTTCAGAATTTTTTCTTTAACAATGAAAAAAAAACTTAAAAAAAAAAATGTTTTTTTTAAAAAAAAAATAAAAAGCAtt                                           >  2:47817/1‑151 (MQ=255)
                                                                                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACACCATAACCCCAAAAAAAAAAAAAAAATTCCGAATTTTTACTTTAACAAAGAAAAAAAAAACTAAAAAAAAAAAAATTTTTTTAAAAAAAAAAATAAAAAGCAtt                                           >  1:830725/1‑151 (MQ=255)
                                                                                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAt                                            <  1:795658/150‑1 (MQ=255)
                                                                                                             aGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAAAACATTTACCTCAAAAAAAAAAAAAACTTCCAGATTTATTAATTTAACAAAGAAAAAAAAAATTAAAAAAAAAAATGTTTTTTTAAAGAAAAAAATAAAAAGCAtt                                           >  1:235485/1‑150 (MQ=255)
                                                                                                             aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTa                                          <  2:459081/151‑1 (MQ=255)
                                                                                                              gAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTaa                                         <  2:877128/151‑1 (MQ=255)
                                                                                                               aaCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCCTTTTCCTCCAAAAAAAAAAAAACTTCCCGATTTTTTTATTTAACAAAGAATAAAAAACTTTAAAAATATAATGTTTTTTTAAAAAAAAAAATAAAAAGCAtt                                           >  2:755190/1‑148 (MQ=255)
                                                                                                                 cTATTCATCAACAGCGTATAGAGGCGGATTTTTAAAAACACCCATTAGCCTCCAAACAAAAAAAAAATTCCACATTTTTTTCTTTAACAATGAAAAAAAAACTTAAAAAATATAATTTTTTTTTTAAAAAAAAAATAAAAAGCATTAAtt                                       >  2:776534/1‑150 (MQ=255)
                                                                                                                   aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTc                                      <  2:471941/149‑1 (MQ=255)
                                                                                                                   aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCaa                                    <  2:510039/151‑1 (MQ=255)
                                                                                                                     tcatcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCCTTTAACTCAAAACAAAAAAAAAATTTCCGATTTTTTACTTTAACAAAAAAAAAAAAACTTAAAAAAAATAATTTTTTTTTAAAAAAAAAAATAACAAGCATTTATTTaa                                    >  2:148809/1‑149 (MQ=255)
                                                                                                                     tcatcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAat                                   <  2:281932/150‑1 (MQ=255)
                                                                                                                       atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAtt                                <  2:351213/151‑1 (MQ=255)
                                                                                                                       atcaACAGCGTATAGAGGCGGGTAAGTTAAAACCATCCATATCCTCCAAAAAAAAAAAAAATTCCAGATTTTTTTATTTAACAACGAAAAAAAAACCTAAAAAAAATAATGTTTTTTTTAAAAAAAAAATAAAAAGGATTTAAttaataat                                 >  2:471608/1‑145 (MQ=255)
                                                                                                                            cAGCGTATAGAGGCGGTTATTTAAAAACCCCCATTAACCCCCCAAAAAAAAAAAACTTCCAGAATTTTTAATTAAACAAAAAAAAAAAAAACTAAAAAAAAAACTTTTATTTTTAAGAAAAAAATAAAAAGCAATAATTCAAAATTaaaa                            >  1:248287/1‑150 (MQ=255)
                                                                                                                               cGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaaa                         <  2:674729/150‑1 (MQ=255)
                                                                                                                               cGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaat                        <  1:429740/151‑1 (MQ=255)
                                                                                                                                  ataGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaataat                     <  2:152621/151‑1 (MQ=255)
                                                                                                                                        gCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa               <  2:769811/151‑1 (MQ=255)
                                                                                                                                         cGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa              <  1:821256/151‑1 (MQ=255)
                                                                                                                                          ggTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaaa             <  2:387601/151‑1 (MQ=255)
                                                                                                                                              aTGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAtt          <  1:851471/150‑1 (MQ=255)
                                                                                                                                                    aaaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt      <  1:612217/150‑1 (MQ=255)
                                                                                                                                                    aaaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa     <  2:341411/151‑1 (MQ=255)
                                                                                                                                                     aaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTa     >  1:396297/1‑150 (MQ=255)
                                                                                                                                                       aCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACg  <  2:845236/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                          
GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACG  >  NC_000913/2903542‑2903841

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: