Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,6920CA45.2% 110.3 / 70.1 84intergenic (+318/+1055)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (2/44);  new base A (38/0);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGA  >  NC_000913/2903542‑2903842
                                                                                                                                                      |                                                                                                                                                          
gTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAc                                                                                                                                                            <  1:546755/151‑1 (MQ=255)
   tttACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACa                                                                                                                                                          >  2:901117/1‑150 (MQ=255)
      aCAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTc                                                                                                                                                       <  1:1181176/150‑1 (MQ=255)
       cAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCAt                                                                                                                                                     <  1:301197/151‑1 (MQ=255)
        aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAACAAATCAtt                                                                                                                                                    >  1:1295106/1‑151 (MQ=255)
        aaGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCAtt                                                                                                                                                    <  2:309246/151‑1 (MQ=255)
         aGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTa                                                                                                                                                   <  2:629994/151‑1 (MQ=255)
          ggAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTa                                                                                                                                                   >  1:976878/1‑150 (MQ=255)
           gAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTTAAAAAACCCATTaac                                                                                                                                                 >  1:576107/1‑149 (MQ=255)
                ttAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACTCATTAACCTCa                                                                                                                                            >  2:88271/1‑150 (MQ=255)
                 tAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCCTTAACCTcca                                                                                                                                           >  1:1115829/1‑148 (MQ=255)
                  aaCTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCaaa                                                                                                                                          <  1:478617/150‑1 (MQ=255)
                     tACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAAACATTTATCTCaaaaaa                                                                                                                                       >  2:417806/1‑148 (MQ=255)
                      acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACACATTATCCTCAaaaaaa                                                                                                                                      >  1:791970/1‑150 (MQ=255)
                        acTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGa                                                                                                                                   <  1:1174818/151‑1 (MQ=255)
                           cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGAc                                                                                                                                  <  2:1300261/149‑1 (MQ=255)
                           cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACaa                                                                                                                                <  2:995541/151‑1 (MQ=255)
                             aTGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACACATTAACCTCaaaaaaaaaaaaaa                                                                                                                               >  1:1082695/1‑150 (MQ=255)
                              tGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACAACCCTTAACCTCCAaaaaaaaaaaaaa                                                                                                                              >  1:936876/1‑149 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTTTCCTCaaaaaaaaaaaaaaaat                                                                                                                            >  1:288670/1‑148 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACt                                                                                                                             <  1:1135864/150‑1 (MQ=255)
                               gTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATACATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACt                                                                                                                             <  2:818437/150‑1 (MQ=255)
                                     aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAATCAATATAATCAAAAAAAAAAAAAAACTCCaca                                                                                                                      >  2:1073623/1‑148 (MQ=255)
                                       aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACTCATTATCCTCAAAAAAAAAAAAAAATTTCAGAttt                                                                                                                   >  2:47996/1‑151 (MQ=255)
                                        aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACCCATTATCCTCCAAAAAAAAAAAAAAATCCAAAtttt                                                                                                                  >  2:1051186/1‑150 (MQ=255)
                                        aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATtta                                                                                                                  <  1:865203/151‑1 (MQ=255)
                                         tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTTAAAACACTCAATATACTCCAAAAAAAAAAAAAAATCCAGAttttt                                                                                                                 >  1:1172223/1‑149 (MQ=255)
                                         tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACACATTTACCTCCAAAAAAAAAAAAAAATTCACATttat                                                                                                                 >  2:20064/1‑151 (MQ=255)
                                         tAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAAATCATTATTCTCAAAAAAAAAAAAAAAATACAGATttat                                                                                                                 >  2:636498/1‑151 (MQ=255)
                                               aCAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                                                          <  1:831372/151‑1 (MQ=255)
                                                cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAACAAACATTTTCCCCAAAAAAAAAAAAAAAATCCCAATTTTTTAtttt                                                                                                          >  2:795234/1‑149 (MQ=255)
                                                cAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTAACCCCCAAAAAAAAAAAAAAATCCAGAAttttttattt                                                                                                           >  1:270406/1‑145 (MQ=255)
                                                   tACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAATCATTATCCTCAAAAAAAAAAAAAAATTCCTAAATTTtttatttacc                                                                                                        >  2:522302/1‑148 (MQ=255)
                                                    aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCa                                                                                                      <  1:1004630/150‑1 (MQ=255)
                                                    aCGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCCATTAACCTCAAAAAAAAAAAAAAATTCCAGATTTTTTTTTTAACCa                                                                                                      >  1:904210/1‑150 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACACATTATCCTCCAAAAAAAAAAAAAAATTCAGATTTATTATTTTAACAAta                                                                                                   >  2:1069684/1‑149 (MQ=255)
                                                       aaTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAACGa                                                                                                  <  2:744645/151‑1 (MQ=255)
                                                        aTCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGaa                                                                                                 <  1:734572/151‑1 (MQ=255)
                                                         tCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCATTAACCTCAAAAAAAAAAAAAAATTTCACAtttttttttttaacaataaa                                                                                                  >  2:661900/1‑136 (MQ=255)
                                                                ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                          <  1:145902/149‑1 (MQ=255)
                                                                ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                          <  1:315273/149‑1 (MQ=255)
                                                                ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAg                                                                                         <  2:58358/150‑1 (MQ=255)
                                                                ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAg                                                                                         <  1:1125323/150‑1 (MQ=255)
                                                                  aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCtt                                                                                      <  2:621979/151‑1 (MQ=255)
                                                                   tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCATTTTCCTCAAAAAAAAAAAAAAAATTCAGATTTTTTATTTTACaaaaaaaaaaaaaaact                                                                                       >  2:830156/1‑136 (MQ=255)
                                                                    cAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaa                                                                                    <  2:1358324/151‑1 (MQ=255)
                                                                     aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTAGACTCAAAAAAAAAAAAAACTTTCCGATTTTTTTTTTTAACAAAGAATAAAAAACTTaa                                                                                    >  1:1195341/1‑150 (MQ=255)
                                                                        tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaaaa                                                                                 <  2:397446/150‑1 (MQ=255)
                                                                          tCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTCCTTAACCTCCAAAAAAAAAAAAAACTTCAGATTTTTTTTTTTACCCAATAAAAAAAAACTTaaaaaaa                                                                               >  1:524029/1‑149 (MQ=255)
                                                                           cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACACATTATCCTCCAAAAAAAAAAAAAATTTCCGATTTTTTATTTTAACAAAGAAAAAAAAAATTaaaaaaaa                                                                              >  2:594009/1‑148 (MQ=255)
                                                                           cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACTCATTATCCTCAAAAAAAAAAAAAAAATCCAGATTTATTTATTTAACAAAGAAAAAAAAAAATaaaaaaaa                                                                              >  2:963317/1‑148 (MQ=255)
                                                                           cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAata                                                                              <  1:88274/150‑1 (MQ=255)
                                                                               tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCCTTATCCTCAAAAAAAAAAAAAAAATCCCGAATTTTTTTTTTACAAAAAAAAAAAAAACATAAAAAAaataatt                                                                         >  2:524523/1‑148 (MQ=255)
                                                                                  aGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTa                                                                      <  1:672327/151‑1 (MQ=255)
                                                                                    tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACTCATTAACCTCAAAAAAAAAACAAAAATCCAAATTTTTTTCTTTAACAAAGAAAAAAAAAACTTAAAAAAATAATGTtttt                                                                    >  1:876966/1‑151 (MQ=255)
                                                                                      tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTTAAAAAACCCATTATTCTCCAAAAAAAAACAAAATTCCCAATTTTTTTTTTTAACAACTAATAAAAAACTTAAAAAAaataattttttttt                                                                  >  1:27247/1‑141 (MQ=255)
                                                                                         aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAGTGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                               <  1:1035058/151‑1 (MQ=255)
                                                                                         aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                                                               <  2:1316996/151‑1 (MQ=255)
                                                                                           tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTCCTTTACCTCCAAAAAAAAAAAAAAATCCACAATTTTTAATTTACCCAATAAAAACAAAATTAAAAAAAAAAATTTTAtttttaaa                                                              >  2:1132104/1‑144 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGa                                                             <  1:714521/150‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGa                                                             <  2:162014/150‑1 (MQ=255)
                                                                                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGa                                                             <  2:1256875/150‑1 (MQ=255)
                                                                                              atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACTCATTAACCCCCAAAAAAAAAAAAAAATCCCGATTTATTACTTAAACAAAGAAAAAAAAACTTAAAAAAAAAAATGTTTTTTTTaaaaaa                                                           >  1:1188852/1‑150 (MQ=255)
                                                                                                    ttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAt                                                       <  2:1015812/148‑1 (MQ=255)
                                                                                                       ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCCATTATCCTCAAAAAAAAAAAAAAAATCCAGATTTATTACTTTACCAATGAATAAAAAACTTAAAAAAAAAAATGTTATTTTTAAAAAAAAAATaaaa                                                  >  2:408045/1‑148 (MQ=255)
                                                                                                         cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAgg                                                <  1:806755/150‑1 (MQ=255)
                                                                                                         cATAAGAACTATTCATCAACAGCGTATAGAGGCGGCGTTTTAAAAACAATCTTTATCCTCAAAAAAAAAACAAACTTACAGGTTTATTACTTTAACCCATAAAAAAAAAAATAAAAAAAATAAAGTTTTTTTTAATAAAAAAATaaaaagc                                               >  2:177994/1‑146 (MQ=255)
                                                                                                           tAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAAACCCACTAGCCTCCAAAAAAAAAAAAACTTTCAGCTTTTTTTCTTAAACAATTAATAAAAAAACTTAAAAAAATAATTTTTTTTTAAAAAAAAAAATAACCGCCAt                                             >  1:1018920/1‑151 (MQ=255)
                                                                                                             aGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTGCAAACCCACCTTATTCACAAAAAAAAAACAAACCTTCCCATTTATTTCTTTAACAAAGAATAAAAAACCTTAAAAAAATATTGTTTTTCTTAAGAAAAAAATaaaaaggatt                                            >  1:338703/1‑142 (MQ=255)
                                                                                                                aCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAAtt                                        <  2:1089924/151‑1 (MQ=255)
                                                                                                                   aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTc                                       <  1:64575/149‑1 (MQ=255)
                                                                                                                   aTTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCaa                                     <  2:711150/151‑1 (MQ=255)
                                                                                                                       atcaACAGCGTATAGAGGCGGTTATTTTAAAACCCTCATTAACCCCAAAAAAAAAAAAAAAATTCCGATTTTTTTTTTTAACAAAAAAAAAAAAAATTAAAAAAAATAATTTTTTTTTTAAAAAAAAAATAAAAAGCATTAAATCAAAAtt                                 >  1:969473/1‑151 (MQ=255)
                                                                                                                       atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATAtt                                 <  2:1298090/151‑1 (MQ=255)
                                                                                                                         caACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaa                               <  1:268500/151‑1 (MQ=255)
                                                                                                                            cAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaa                            <  2:559114/151‑1 (MQ=255)
                                                                                                                            cAGCGTATAGAGGCGGTTATGTAAAAACACCCATTAGCCTCCAAAAAAAAAAAAAATTCCAGATTTATTACTTAACCAAAGAAAAAAAAACTTAAAAAAAATAATGTTATTTTAAAGAAAAAAATAAAAAGCATTAATTTAATATTaaaa                             >  1:1326920/1‑150 (MQ=255)
                                                                                                                                      agGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAt                  <  1:430148/151‑1 (MQ=255)
                                                                                                                                      agGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAt                  <  2:119088/151‑1 (MQ=255)
                                                                                                                                           gTTATGTAAAACGACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaaa              >  2:787721/1‑150 (MQ=255)
                                                                                                                                              aTGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt          <  1:1363605/151‑1 (MQ=255)
                                                                                                                                              aTGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAAttt          <  1:267981/151‑1 (MQ=255)
                                                                                                                                                  aaaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt       <  1:540608/150‑1 (MQ=255)
                                                                                                                                                  aaaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt       <  1:1216811/150‑1 (MQ=255)
                                                                                                                                                  aaaaCCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCtt       <  2:529314/150‑1 (MQ=255)
                                                                                                                                                      ccACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGa  <  1:289048/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                          
GTTTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGA  >  NC_000913/2903542‑2903842

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