Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG40.6% 100.8 / 10.6 65F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (13/25);  new base G (26/0);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 2.59e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508509‑3508805
                                                                                                                                                    |                                                                                                                                                    
aGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaa                                                                                                                                                    <  2:870633/151‑1 (MQ=255)
aGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaaa                                                                                                                                                    <  1:422900/151‑1 (MQ=255)
  tCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCCGgaaa                                                                                                                                                   >  2:192932/1‑150 (MQ=255)
   ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGgaaaa                                                                                                                                                  >  2:396157/1‑150 (MQ=255)
    cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAtt                                                                                                                                                <  1:978916/151‑1 (MQ=255)
       gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcac                                                                                                                                             <  1:741964/151‑1 (MQ=255)
         tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcaca                                                                                                                                            <  2:253232/150‑1 (MQ=255)
          gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAttta                                                                                                                                              >  2:48268/1‑145 (MQ=255)
          gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacac                                                                                                                                           <  1:45150/150‑1 (MQ=255)
              ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAAgg                                                                                                                                       >  2:1038492/1‑150 (MQ=255)
               tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGt                                                                                                                                      >  2:515896/1‑150 (MQ=255)
                  aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACAAACAGGAGcc                                                                                                                                  >  1:312172/1‑151 (MQ=255)
                  aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGcc                                                                                                                                  <  1:944357/151‑1 (MQ=255)
                    ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGTAGcccc                                                                                                                                >  2:721560/1‑151 (MQ=255)
                     cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGccccc                                                                                                                               >  1:667689/1‑151 (MQ=255)
                     cGGAGAGATCAATATGCAGATGCTCATCGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAATCCGGAGCGAGGAAAGAAAAGAGCGTCCCGTTTGTGGCGGTGGCGCGCGGAAAAAGTCACACACGGTAGccccc                                                                                                                               <  1:428666/151‑1 (MQ=255)
                      ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGCAGcccccc                                                                                                                              >  1:1041344/1‑150 (MQ=255)
                           gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGCAAAAATCACACACGGGAGCCCCCACATc                                                                                                                          >  1:371919/1‑150 (MQ=255)
                           gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACAAAAGGTAGCCCCCACATc                                                                                                                          >  2:645603/1‑150 (MQ=255)
                             tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGGGGGAAAAATTACACAACGGAGccccccc                                                                                                                             >  2:105993/1‑143 (MQ=255)
                                atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCACATCAGaaa                                                                                                                     >  2:52548/1‑150 (MQ=255)
                                 tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCACATCAGaaa                                                                                                                     >  2:56294/1‑149 (MQ=255)
                                      aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGTGGGCCGAAAAATTCACACAAGGTAGTCCCCCCATCAGAAAGCgtggg                                                                                                              >  2:140428/1‑150 (MQ=255)
                                       gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGACCCCACATCAGAAAGCggggtg                                                                                                             >  1:32730/1‑151 (MQ=255)
                                          gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTACACACGGGAGCCCCCACCCCCGAAAACGGGGGgcg                                                                                                           >  1:13653/1‑150 (MQ=255)
                                          gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATTAAAAAAGATAACCCCCACAACCGAAACGgtggg                                                                                                              >  2:674983/1‑146 (MQ=255)
                                           cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGGGCGCGGGAAAAATCACACAACGGAGCCCCCACCTCCAAAAGCgggg                                                                                                               >  1:980317/1‑145 (MQ=255)
                                             cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCACAACCGAAAACGGGGGGCGcac                                                                                                        >  1:587089/1‑150 (MQ=255)
                                             cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcac                                                                                                        <  1:477999/150‑1 (MQ=255)
                                              aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCCGAAAACGTGGGGCGcaca                                                                                                       >  2:894043/1‑150 (MQ=255)
                                                   cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                 <  2:1006010/151‑1 (MQ=255)
                                                    gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                <  2:37409/151‑1 (MQ=255)
                                                    gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGACGCGCGGGAAAATTTACAAACCGTGACCCCCCCATCAGGAAACGCGGGgcgtgcccac                                                                                                    >  1:406028/1‑140 (MQ=255)
                                                             aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                       <  1:936366/151‑1 (MQ=255)
                                                             aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                                       <  2:628103/151‑1 (MQ=255)
                                                              cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                      <  2:357862/151‑1 (MQ=255)
                                                                cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAAGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCCGCGTAGCGatc                                                                                     >  1:286790/1‑150 (MQ=255)
                                                                cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                                                                    <  2:684523/151‑1 (MQ=255)
                                                                    ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATccc                                                                                >  2:862029/1‑151 (MQ=255)
                                                                    ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCCTCAGAAAACGGGGGGGGGACACCGAGCTGCGTAGCGATCATccc                                                                                >  1:226636/1‑151 (MQ=255)
                                                                        cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCCCACCCAGCTGCGGAGCGATCATccc                                                                                >  2:409283/1‑147 (MQ=255)
                                                                              tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACAAACGGTAGCCCCCACCTCAGAAAGCGGGGTGCGCACAACGCGCCGCGTGGCGAGCATccccttct                                                                           >  1:344960/1‑141 (MQ=255)
                                                                               cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGTAGCCCCCACATCAGAAAACGTGGGGGGCACACCCAGCCGCGTAGCGATCATCCCCTGCTGATaa                                                                      >  1:6044/1‑150 (MQ=255)
                                                                                    ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                                <  1:752251/151‑1 (MQ=255)
                                                                                       gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                              <  1:312624/150‑1 (MQ=255)
                                                                                       gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                             <  2:1058497/151‑1 (MQ=255)
                                                                                       gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTAAACACGGTAGCCCCCACAACAAAAAAAGGTGGGCGGACACCCAACCGCGGAGCGATCATCCCTTTCCGATTAAACAGCAc                                                              >  2:651139/1‑150 (MQ=255)
                                                                                         ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGTAGCCCCCACACCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACTgg                                                           >  2:762271/1‑151 (MQ=255)
                                                                                           ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                            <  1:120555/148‑1 (MQ=255)
                                                                                           ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca                                                         <  1:928900/151‑1 (MQ=255)
                                                                                             tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAACGGGGTGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAGCACTGGgac                                                        >  2:481845/1‑150 (MQ=255)
                                                                                                aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACAGGAGCCCCCACAACAGAAAACGGGGGGCGCCCACCCAGCTGGGGAGCGATCATCCCTTGCTGATAAACCACCACCGGcacccc                                                     >  2:196828/1‑148 (MQ=255)
                                                                                                    cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCCGCGTAGCGATCATCCCTTGCTCATAAAACAGCACTGGCCCCACTTcgc                                                >  1:34177/1‑151 (MQ=255)
                                                                                                    cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGGAGCCCCCCCATCAGAAAGCGGGGTGGGCACACCGAGCTGCGGAGCGGTCATCCCTTTCTGGTAAAACAACACCGGcaccac                                                     >  1:709105/1‑146 (MQ=255)
                                                                                                      gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAACACGGTAGCCCCCCCATCAGAAAACGGGGGGGGCACACCCCGCTGCGTAGGGATCATCCCTTGCTGATAAAACAGCACCGGCACCACTTCGCgg                                              >  2:83596/1‑151 (MQ=255)
                                                                                                      gagggaggAAAGAAAAGAGCGACCATTTGGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                             <  1:396093/151‑1 (MQ=255)
                                                                                                         ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACGCCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTa                                           <  1:167386/151‑1 (MQ=255)
                                                                                                         ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGAGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTa                                           <  1:674875/151‑1 (MQ=255)
                                                                                                             gaaagaaaAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTAAACACCGTAGGCCCCACATCAGAAAGAGTGGGGGGCACACCGAGATGCCGAGCGATCATCCCTTGTTGATAAACCAACACTGGATACACTTCTCGGGTagg                                         >  2:625968/1‑147 (MQ=255)
                                                                                                              aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGAtggtg                                      >  2:684286/1‑151 (MQ=255)
                                                                                                                       agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                             <  2:750311/151‑1 (MQ=255)
                                                                                                                                   gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                 <  1:544350/151‑1 (MQ=255)
                                                                                                                                       cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg             <  2:313965/151‑1 (MQ=255)
                                                                                                                                        ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc            <  1:1038360/151‑1 (MQ=255)
                                                                                                                                         gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc           >  1:37026/1‑151 (MQ=255)
                                                                                                                                                gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc    >  2:558810/1‑151 (MQ=255)
                                                                                                                                                 cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg   >  1:821152/1‑151 (MQ=255)
                                                                                                                                                  ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc  >  2:1009451/1‑151 (MQ=255)
                                                                                                                                                  ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc  >  1:994861/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGC  >  NC_000913/3508509‑3508805

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