Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6920CA48.6% 65.2 / 58.6 70intergenic (+318/+1055)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (2/34);  new base A (33/1);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 3.83e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATT  >  NC_000913/2903553‑2903827
                                                                                                                                           |                                                                                                                                               
gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGc                                                                                                                                   <  1:483604/151‑1 (MQ=255)
  atatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCATTATcc                                                                                                                                  >  2:267950/1‑150 (MQ=255)
   tatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACGCATTattt                                                                                                                                  >  1:527797/1‑146 (MQ=255)
    atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTc                                                                                                                                <  2:425392/150‑1 (MQ=255)
    atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCa                                                                                                                               <  2:864925/151‑1 (MQ=255)
       aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaaaa                                                                                                                            <  1:381388/151‑1 (MQ=255)
          tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTTACCTCaaaaaa                                                                                                                          >  2:966505/1‑148 (MQ=255)
            cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAg                                                                                                                       <  2:342948/151‑1 (MQ=255)
             acTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAg                                                                                                                       <  2:765794/150‑1 (MQ=255)
              cTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACACCTTTACCTCCAaaaaaaaa                                                                                                                      >  2:129677/1‑148 (MQ=255)
                cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACAACCATTAACCTCaaaaaaaaaaaa                                                                                                                    >  2:701821/1‑146 (MQ=255)
                cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCCATTAACCTCAAAACaaaaa                                                                                                                      >  2:1008803/1‑146 (MQ=255)
                 gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCAATAACCCCaaaaaaaaaaaaa                                                                                                                   >  1:651453/1‑145 (MQ=255)
                 gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTAGCCTCCAaaaaaaaaaaa                                                                                                                   >  2:152364/1‑145 (MQ=255)
                 gATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaaa                                                                                                                 <  1:966026/151‑1 (MQ=255)
                  aTGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAACACTCATTAACCTCAAAACAaaaaaaaa                                                                                                                 >  2:60256/1‑150 (MQ=255)
                      cTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTTTTTATAAAAAATCATTAATCTCAAAAAAAAAAAAAAcattc                                                                                                           >  1:1058679/1‑147 (MQ=255)
                         tAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAg                                                                                                         <  2:365313/150‑1 (MQ=255)
                          aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGAt                                                                                                      <  1:265037/151‑1 (MQ=255)
                             aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGGAAAAACACTCCTTAACCTCAAAACAAAAACAAAAATCCCGAttt                                                                                                    >  2:1032753/1‑150 (MQ=255)
                                    aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCATTTTACTCAAAAAAAAAACAAACTTCCAGATttatt                                                                                                 >  1:516802/1‑146 (MQ=255)
                                    aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAAAAACCATTAAACTCAAAACAAAAAAAAAATTCCAGATTTATTAAtt                                                                                             >  1:1013158/1‑150 (MQ=255)
                                    aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                                            <  2:447626/151‑1 (MQ=255)
                                     cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTaa                                                                                           <  2:217042/151‑1 (MQ=255)
                                        tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAAAAATCATTAGCCCCCAAAAAAAAAAAAACATCCCGAtttttttttt                                                                                             >  2:718262/1‑142 (MQ=255)
                                            aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                                                   <  2:641685/151‑1 (MQ=255)
                                               cGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAAt                                                                                 <  1:251844/150‑1 (MQ=255)
                                                  aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACCCAATTACCCCAAAAAAAAAAAAAAACTTCAAATTTATTACTTTACCAAAGAaaaa                                                                               >  2:41827/1‑149 (MQ=255)
                                                  aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACa                                                                             <  2:1019180/151‑1 (MQ=255)
                                                     ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCCTTATACTCAAAAAAAAAAAAAACATTCAGATTTTTTACTTAACaaaaacaaaaaaa                                                                            >  1:409095/1‑139 (MQ=255)
                                                     ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAg                                                                           <  2:884640/150‑1 (MQ=255)
                                                      tATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGc                                                                          <  2:74472/150‑1 (MQ=255)
                                                       aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCtt                                                                        <  2:458259/151‑1 (MQ=255)
                                                        tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTAACCTCAAAAAAAAAAAAAAAATCCCAAATTATTATTTAAACaaaaaaaaaaaaaaat                                                                        >  2:285143/1‑136 (MQ=255)
                                                         cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaa                                                                      <  2:517095/151‑1 (MQ=255)
                                                            ttttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAAAACTCATTACCATCAAAAAAAAAAAAAACCTTCTTACTTAATTACTTCAACAATAATAAAAAACTTaaa                                                                     >  1:205380/1‑149 (MQ=255)
                                                             tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaaaa                                                                   <  2:1135686/150‑1 (MQ=255)
                                                                cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAata                                                                <  2:174172/150‑1 (MQ=255)
                                                                 aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAAtata                                                              <  2:51274/151‑1 (MQ=255)
                                                                    tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTTACCTCCAAAAAAAAAAAAACATCCCGATTTTTTATTTAAACAAATAAAAAAAAAATTaaaaaaaa                                                                >  2:252929/1‑144 (MQ=255)
                                                                     aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAAATCATTAGTCTCAAAAAAAAAAAAAAATTCCAAATTTATTATTTTACCAATGAATAACAAGATTAAAAAAaataat                                                            >  1:428552/1‑147 (MQ=255)
                                                                      aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                                                        <  1:930409/151‑1 (MQ=255)
                                                                      aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTGCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTg                                                          <  2:870134/149‑1 (MQ=255)
                                                                        gTATGAACTGGATATCTTTATCCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAt                                                      <  1:756091/151‑1 (MQ=255)
                                                                              aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTATCCTCAAAAAAAAAACAAAAATTCAGAATTATTACTTTACCAAAAAATAAAAAACTTAAAAAAAAAAATTTTAtttt                                                   >  2:951642/1‑146 (MQ=255)
                                                                                tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTTTCCTCCAAAAAAAAAAAAAATTCCAAATTTATTTATTAAACCAAAAAAAAAAAACTTTAAAAAAAAAATGTTTTTTTaaa                                                >  2:750605/1‑149 (MQ=255)
                                                                                      tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCATTAGCCCCAAAACAAAGAGAAACTTCACGGTTTTTTTATTAACCAAATAAAAAACAAATTAAAAAAAAAACTTTTATTTTTCaaaaaaaa                                         >  2:328399/1‑150 (MQ=255)
                                                                                        tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACAAATTATCACAAAAAAAAAAAAAAAATTCCCGATTTATTTTTTTAACAAATAAAAACAAAATTAAAAAATAAAATTTTATTTTTAAGaaaaaa                                         >  1:562897/1‑148 (MQ=255)
                                                                                         ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATCCACCTCAAAAAAAAAAAAAAAATTCCGAATTATTTATTAAACAAAGAATAAAAAAATTTAAAAAAATAATTTTATTTTaaaaaaaaaaat                                       >  2:59346/1‑147 (MQ=255)
                                                                                         ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaaa                                         <  2:514613/147‑1 (MQ=255)
                                                                                         ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaaa                                         <  2:554646/147‑1 (MQ=255)
                                                                                          tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTaa                                     <  1:206392/150‑1 (MQ=255)
                                                                                           aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAAc                                    <  1:397406/150‑1 (MQ=255)
                                                                                             tCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCCATTAAACTCCAAAAAAAAAAAAAATTCCCGAATTTTTTTTTTACCAAATAAAAAAAAAACTTAAAAAAATACTTCTTTTTTTAAAAAAAAAATaaaa                                    >  2:211720/1‑148 (MQ=255)
                                                                                              cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAGCCTCCAAAAAAAAAAAAAATTCCAGATTTTTTTCTTTAACAAAGAAAAAAAAACTTAAAAAAAATAATTTTATTTTTAAAAAAAAAATaa                                     >  2:965115/1‑146 (MQ=255)
                                                                                                   gAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTaa                           <  2:454204/151‑1 (MQ=255)
                                                                                                   gAACTATTCATCAACAGCGTATAGAGGCGGATATGTAAAAAAACACCATAAACCCCAAAAAAAAAAAAACTTTCCGAGTTATTAATTAACCAAATTATAACAAAATTAAAAAAAATAATGTTATTTTTAAGAAAAAAATAACaaagaata                            >  2:972184/1‑143 (MQ=255)
                                                                                                    aaCTATTCATCAACAGCGTATAGAGGCGGATATTAAAACACACCCCTTAAACCCAAAATAAAAACAAAACTACCAACTTTTATCTTTACCCAAAAAAAAAAAAAATTAAAAAAAATACTTTTTTTTTAAAAAAACAATAACAAGCAttta                           >  1:492002/1‑148 (MQ=255)
                                                                                                     aCTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTAACCCCAAAAAAAAAAAAAAAATCCAGATTTTTTTATTAAACAAAAAAAAAAAAACTTAAAAAAAATACTGTTATTTTTAAAAAAAAAATAAAAAGCATTAAtt                         >  2:669575/1‑151 (MQ=255)
                                                                                                        aTTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTAACCTCAAAAAAAAAAAAAAAATCCAAAATTTTTTTTTTAACAAAGAAAAAAAAACTTAAAAAAAATACTTTTTTTTTTAAAAAAAAAATAAAAAGCAATAAtt                         >  2:462706/1‑148 (MQ=255)
                                                                                                        aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCa                       <  2:503374/150‑1 (MQ=255)
                                                                                                         ttCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCACTAAACTCAAAAAAAAAAAAAAAATCCAAAATTTTTTTTTAAACCATAAAAAAAAAACTTAAAAAAAAAAATTTTTTTTTTAAAAAAAAAATAAAAGGGAATaa                           >  2:188083/1‑145 (MQ=255)
                                                                                                         ttCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACACCCTTGACCCAAAAAAAAAAAAAAATTCCAGATTTTTTTATTAAACAAAAAATAAAAAACATTAAAAAAATAATATTTTTTTTAAAAAAAAAATAAAAAGCAATAATTTaa                      >  2:743731/1‑150 (MQ=255)
                                                                                                          tcatcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAat                     <  2:510356/150‑1 (MQ=255)
                                                                                                           catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                   <  1:152316/151‑1 (MQ=255)
                                                                                                             tcaACAGCGTATAGAGGCGGTTATGTAAAAAAACCCATTAGCCTCAAAAAAAAAAAAAAATTCCAGGTTTTTTACTTTAACAAAGAATAACAAACTTAAAAAAAATAATGTTATTTTATAGAAAAAAATAAAAAGAATTaaataaa                      >  1:1133293/1‑141 (MQ=255)
                                                                                                              caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTa                 <  2:697448/150‑1 (MQ=255)
                                                                                                                  aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa            <  1:615648/151‑1 (MQ=255)
                                                                                                                             gCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt  >  1:963047/1‑150 (MQ=255)
                                                                                                                                 ttATGTAAAAACACTCATTAGCCTCAAAACAAAGAAAAACATCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt  <  1:325805/146‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATT  >  NC_000913/2903553‑2903827

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