Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA43.9% 83.6 / 34.3 66E384D (GAA→GATyhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (9/28);  new base A (29/0);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 3.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAT  >  NC_000913/3508518‑3508808
                                                                                                                                                  |                                                                                                                                                  
tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAAttaaa                                                                                                                                                >  1:1132320/1‑147 (MQ=255)
tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATcaaa                                                                                                                                               >  2:716774/1‑148 (MQ=255)
  tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCAACCATTTTGTGGCGGTGGCGCGCGGGAAAAATcacaca                                                                                                                                             >  2:852480/1‑150 (MQ=255)
      tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACgg                                                                                                                                          <  1:258349/149‑1 (MQ=255)
      tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGt                                                                                                                                         <  1:719665/150‑1 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                        <  2:609867/150‑1 (MQ=255)
       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                       <  2:386976/151‑1 (MQ=255)
        aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGt                                                                                                                                         >  1:729800/1‑148 (MQ=255)
             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGccccc                                                                                                                                  >  2:952005/1‑150 (MQ=255)
              gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcac                                                                                                                                <  2:208146/151‑1 (MQ=255)
                  gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                            <  1:1197623/151‑1 (MQ=255)
                   aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAACACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGCAAAATTACCCACAGTCGGCCCCCCAACa                                                                                                                            >  2:703164/1‑150 (MQ=255)
                   aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGGAGccccccc                                                                                                                                >  2:573231/1‑144 (MQ=255)
                      aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTACACACGGTAGCCCCCCCATCaaaa                                                                                                                         >  2:478469/1‑150 (MQ=255)
                       atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCCTCAGcaa                                                                                                                        >  1:485691/1‑150 (MQ=255)
                         atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                      <  1:101970/150‑1 (MQ=255)
                         atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                      <  1:57481/150‑1 (MQ=255)
                                 gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                              <  1:286206/150‑1 (MQ=255)
                                   tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACGGGAGCCCCCCCATCAGAAAACgtgggggg                                                                                                              >  2:1161258/1‑144 (MQ=255)
                                    cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGAGGGGGGGGcac                                                                                                           >  2:808937/1‑150 (MQ=255)
                                     aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGCGGGgcgcgc                                                                                                           >  1:982292/1‑147 (MQ=255)
                                      tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACGGGAGCCCCCACATCAGAAAACGGGGGGGGGACAcc                                                                                                        >  2:96865/1‑151 (MQ=255)
                                          cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                    <  2:429779/151‑1 (MQ=255)
                                          cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                    <  1:154327/151‑1 (MQ=255)
                                              ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGCACACCGAGCTg                                                                                                  >  2:113218/1‑149 (MQ=255)
                                              ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGc                                                                                                 >  2:1178490/1‑150 (MQ=255)
                                                     cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAAACGTAGCCCCCACCACAGAAAACGTGGTGCGCACCCCGGGCTg                                                                                                  >  2:721341/1‑142 (MQ=255)
                                                        gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTAAACACAGTAGCCCCCACACCAGAAAACGTGGGGCGCCCACCCAgcggcggggcg                                                                                           >  1:618310/1‑141 (MQ=255)
                                                         tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGCGGGAAAAATTAAACACGGGAGCCCCCCCACCAAAAAACGGGGGGGGCCCACCCAGGCGCGGAGCGatc                                                                                        >  1:1196059/1‑148 (MQ=255)
                                                             aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCAGCCGCGGAGCGGTCATcccct                                                                                 >  1:393395/1‑149 (MQ=255)
                                                             aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                 <  1:524566/151‑1 (MQ=255)
                                                             aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATTAAACACGGGAGCCCCCCCCACAGAAAACGGGGTGGGCACACCCAGCCGCGTAGCGATCAAccc                                                                                   >  1:1035118/1‑149 (MQ=255)
                                                               cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                               <  2:677719/151‑1 (MQ=255)
                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCATCAGAAAGCGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCt                                                                              >  2:760983/1‑151 (MQ=255)
                                                                    cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAACGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCCGCGGAGCGATCATCCCTTTCTGAt                                                                           >  1:1043559/1‑150 (MQ=255)
                                                                       aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATTAAACACGGGAGCCCCCCCAACAGAAAACGGGGGGGGGGCAACCAGCTGCGGAGCGAGCATCCCCTTCTGATaaa                                                                        >  1:38048/1‑150 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCCTCAAAAAACGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATaaaac                                                                      >  1:485623/1‑148 (MQ=255)
                                                                         aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                     <  1:737835/151‑1 (MQ=255)
                                                                             tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACAACGGAGCCCCCCCACCAGAAAAGGGGGGGCGCACACCCAGGTGCGGAGCGATCAACCCTTGCTGATaaa                                                                        >  1:1171905/1‑144 (MQ=255)
                                                                               tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTAAACACGGGAGCCCCCCAATAAGAAAACGTGGTGGGCCCACCGAGCTGCGTAGCGATCACCCCTTGCTGTTAAAACa                                                                     >  2:108362/1‑145 (MQ=255)
                                                                                 tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTgg                                                              >  2:1203341/1‑150 (MQ=255)
                                                                                  ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGc                                                             >  2:340319/1‑150 (MQ=255)
                                                                                  ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                              <  1:33938/149‑1 (MQ=255)
                                                                                  ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                               <  2:233266/148‑1 (MQ=255)
                                                                                             gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                 <  2:309081/151‑1 (MQ=255)
                                                                                             gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                 <  1:197749/151‑1 (MQ=255)
                                                                                              agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAACCAACACTGGCACCACTTCGCggg                                                >  1:1112258/1‑151 (MQ=255)
                                                                                                ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCCGGAAAAATTACACACCGGAGCCCCCCCATACGAAAGAGTGGGGGGCACACCCCGCCGCGGAGCGATCATCCCCTTCTGGTAAAACAACACTGGgacccccttccgggg                                               >  1:134578/1‑139 (MQ=255)
                                                                                                  aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACACCCGACTGCGTGGCGGTCATCCCCTGCTGATAAAACAGCACCGGCACCACCTCGCggg                                                >  1:592205/1‑147 (MQ=255)
                                                                                                    gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                          <  1:1090975/151‑1 (MQ=255)
                                                                                                     aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATAAAACACGGTAGGCCCCCCATTAGAAAACATGGGGCGCACACCAAGATGTGGAGCGAATATGTCTTGCTGGTTAACGACAACCGGCACCCCCTCGCGGGTAtgggg                                         >  2:624908/1‑149 (MQ=255)
                                                                                                      aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAACCCCCCCATCAGAAAACGGGGTGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAACACTGGCACCCCCTCGCGGGGAtggggg                                        >  2:609851/1‑151 (MQ=255)
                                                                                                       agaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGCCCCCACATCAGAAAACGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGAtggggg                                        >  2:39926/1‑150 (MQ=255)
                                                                                                                   ccATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                           <  2:276374/151‑1 (MQ=255)
                                                                                                                       tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                        <  2:1115736/150‑1 (MQ=255)
                                                                                                                        ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                       <  2:894108/150‑1 (MQ=255)
                                                                                                                         tgtgGCGGTGGGGCGCGGAAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGTGGGCACAACCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACCTCGCGGGGATGGTGGGTGTGACCAATGGTCgg                       >  1:663218/1‑149 (MQ=255)
                                                                                                                          gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGTTGTGACCAATGGTCGGATc                    <  2:986414/151‑1 (MQ=255)
                                                                                                                              cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                <  1:156049/151‑1 (MQ=255)
                                                                                                                                  ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAATCCAGCACTGGCAACACTTCGCGGGTATGGTGGCTGTGACCAATGGGCGGATCGTTGCCGt            <  1:777475/151‑1 (MQ=255)
                                                                                                                                  ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg             <  1:654602/150‑1 (MQ=255)
                                                                                                                                    cgcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg           >  2:530265/1‑150 (MQ=255)
                                                                                                                                         cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg      >  2:1168696/1‑151 (MQ=255)
                                                                                                                                           ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc     <  2:29960/151‑1 (MQ=255)
                                                                                                                                            gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc    >  2:171749/1‑151 (MQ=255)
                                                                                                                                             aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc    >  2:555071/1‑150 (MQ=255)
                                                                                                                                              aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt  <  1:459023/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAT  >  NC_000913/3508518‑3508808

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: