Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,619,2300TC40.9% 147.5 / 13.6 88I311V (ATT→GTT) uraAuracil permease
Reads supporting (aligned to +/- strand):  ref base T (41/11);  new base C (0/36);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 2.85e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGA  >  NC_000913/2619081‑2619375
                                                                                                                                                       |                                                                                                                                                 
gtAATCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAAt                                                                                                                                                   <  2:343375/150‑1 (MQ=255)
gtAATCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATc                                                                                                                                                  >  2:170218/1‑151 (MQ=255)
gtAATCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATc                                                                                                                                                  <  1:1469889/151‑1 (MQ=255)
  aaTCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATcgc                                                                                                                                                >  2:810148/1‑151 (MQ=255)
  aaTCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATcg                                                                                                                                                 <  2:1468454/150‑1 (MQ=255)
   aTCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATcgcc                                                                                                                                               <  2:205517/151‑1 (MQ=255)
   aTCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATcgcc                                                                                                                                               <  1:992193/151‑1 (MQ=255)
          ttcgattcgaTCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGccgccccg                                                                                                                                         <  2:484477/150‑1 (MQ=255)
           tcgattcgatCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGccgccccgc                                                                                                                                        <  2:1193429/150‑1 (MQ=255)
              attcgatCAAAACACGAATACTGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGCGATCATCTGGATA‑GCGGCAGACAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGGCGCCTCGCCGa                                                                                                                                     >  1:302647/1‑150 (MQ=255)
              attcgatCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGAt                                                                                                                                    >  2:123665/1‑151 (MQ=255)
              attcgatCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGAt                                                                                                                                    >  1:710419/1‑151 (MQ=255)
                 cgatCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAAc                                                                                                                                 >  2:381152/1‑151 (MQ=255)
                  gatCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAAcc                                                                                                                                >  1:1026049/1‑151 (MQ=255)
                     cAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAAccc                                                                                                                               >  2:283551/1‑149 (MQ=255)
                     cAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAAccc                                                                                                                               >  2:1210924/1‑149 (MQ=255)
                      aaaaCACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAgg                                                                                                                            >  1:584327/1‑151 (MQ=255)
                       aaaCACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGt                                                                                                                           >  2:310611/1‑151 (MQ=255)
                            cGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGt                                                                                                                      >  2:126470/1‑151 (MQ=255)
                                 aCCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAac                                                                                                                  >  1:1130763/1‑150 (MQ=255)
                                  ccGAAAGCCCCGATGACACAATAAAGCAGCAGCGAACCCCCGCCCACCCCCGACAATGGGATCTTCTGTATA‑GCGGCACCCAGTTT‑CCCGACGCAGGAAAGGAGGAAAGCGAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAac                                                                                                                  <  2:1106311/149‑1 (MQ=255)
                                    ggAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAAGCCAGGTACTGTAAACACgg                                                                                                              >  1:64946/1‑151 (MQ=255)
                                       aGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTg                                                                                                           >  2:617774/1‑151 (MQ=255)
                                           ccGATGACACCAAAAAGCGGCGGCAAACCCCCCCCCACCCCCGGCAATGGGATCATGGGGATA‑GCGGCAGCCGGTTC‑CCCGACGCAGGAAAGCAGGATAGCAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCg                                                                                                       <  2:601780/151‑1 (MQ=255)
                                             gaagccacCATAAAGCACCAGCGAACCCCCCCCCCCCCCCGGAAGTGGTATCTTCTGGATA‑GCGGCACCCATTTA‑ACCGACGCAGAAAAGCAGAAAGGCAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGa                                                                                                          <  1:1305116/141‑1 (MQ=255)
                                             gATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGcc                                                                                                     >  1:1060249/1‑151 (MQ=255)
                                              aTGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGcc                                                                                                     >  2:1380208/1‑150 (MQ=255)
                                                gACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATc                                                                                                  >  1:785714/1‑151 (MQ=255)
                                                aacacAAAAAAGCAGCAGAAAAACCCCCCCCATCACGGGAAAGGGATTCTTCTGGATG‑GCGGCAGCAATTTC‑CCCGACCCAGAAAAGCAGGATAGGAAAAACCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCAt                                                                                                   <  1:1004177/149‑1 (MQ=255)
                                                 acaacaAAAAGCAGCAGCGACCCGCCTCCTACCACCGGCAATGGGGTCATCGGGATA‑GCGGCCCCCATTTA‑GCGGAAGCAGGAAAGCAGAAAAGCGAAAACCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCa                                                                                                    <  1:352240/148‑1 (MQ=255)
                                                  cacCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCAc                                                                                                >  1:279/1‑151 (MQ=255)
                                                      aTAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGcc                                                                                             >  1:1239153/1‑150 (MQ=255)
                                                       aaaaGCGGAAGCAAACCCCCCCCCACCCCCGGAAAGGGTATCATCTGGATG‑GCGGCGGCCATTTC‑CCCGACCAAGGAAAGCAGGAAAGCAAAACTCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCaa                                                                                           <  2:1568094/150‑1 (MQ=255)
                                                            cagcagcGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGttt                                                                                      >  2:1047409/1‑151 (MQ=255)
                                                               cagcGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTtct                                                                                   >  2:777021/1‑151 (MQ=255)
                                                                    aaccgccgcCCACCCCCGGCAGGGGTATCTTCAGAAGG‑GGGGCGGCCATTTA‑CCCGACGCGGGAAAGAAGAATAGCAAAAGCCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTtct                                                                                   <  1:697114/143‑1 (MQ=255)
                                                                    aacccccgccCATCCCCGGAAATGGGATCATCTGGATG‑GCGGCACCCATTTT‑CCCGACGCAGAAAAGCAGGATAGAAAAAACCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTtctc                                                                                  <  2:269922/142‑1 (MQ=255)
                                                                    aacccccccccACCACCGGAAAGGGGATACTCGGAATA‑GC‑GCCGCCAATTTCCCCGACGCAGAAAAGCAGGATACAAAAACCCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTa                                                                              <  1:304154/143‑1 (MQ=255)
                                                                    aaaCGCCCCCCACCGCGGAAAATGGGATCAT‑TGGAAAGGGGGCAGCTATTTC‑CCCGACGCAGGAAGGAGGGAGAGGAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCg                                                                                <  2:1013920/149‑1 (MQ=255)
                                                                     aaCGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTa                                                                              >  2:941543/1‑150 (MQ=255)
                                                                         ccgccCACCCCCGGCAATGGGATCATCTGGATG‑GGGGCAGCCAGTTT‑CCCGACGCAGGAAAGCAGGATAGAAAAACCCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAagt                                                                           <  2:120027/149‑1 (MQ=255)
                                                                         ccgccCACCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCTTTTT‑CCCGACGCAGAAAAGCAGGATAGAAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAagta                                                                          <  2:1261971/150‑1 (MQ=255)
                                                                         ccccccATCCCCGGAAAGGGGACCATCTGGATA‑GCGGCCGCCAGTTT‑CCCGCCGCAGAAAAGCAGGATAGCAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAagta                                                                          <  2:861572/147‑1 (MQ=255)
                                                                           gccCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTTCTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAagtagt                                                                        >  2:927296/1‑150 (MQ=255)
                                                                               aTCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTg                                                                   >  1:1082699/1‑151 (MQ=255)
                                                                                 cgccggaagaggaatcagcTGGATG‑GCGCCAGCACTTTC‑CCCGACGCAGGAAAGCAGAATAGCAAAAACCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTg                                                                   <  1:1194019/136‑1 (MQ=255)
                                                                                 ccccGGAAATGGGATCATCTGAATA‑GCGGCAGCCATTTT‑CCCGCCGCAGAAAAGCAGGATAGCAAAAACCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTg                                                                   <  2:576409/147‑1 (MQ=255)
                                                                                  cccGGCAATGGTATCATTTGGATA‑GCGGCAGCCGTTTT‑CCCGAGGCAGGAAAGAAGGATAGCAAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTgg                                                                  <  2:752170/148‑1 (MQ=255)
                                                                                   ccGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGt                                                               >  2:968684/1‑151 (MQ=255)
                                                                                   ccGGAAATGGAATCATGTGGATA‑GGGGCGGCCATTTT‑CCCGACGAAGAAAAGCAGGATAGAAAAAACCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGt                                                               <  2:35248/151‑1 (MQ=255)
                                                                                    cGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGt                                                               >  2:1280849/1‑150 (MQ=255)
                                                                                    cGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGt                                                               >  1:1531248/1‑150 (MQ=255)
                                                                                        aaggggATCATCTGGATA‑GCGGCGCCCATTTT‑CCCGACGCAGGAAAGCAGGATAGCAAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGc                                                          <  2:1191776/148‑1 (MQ=255)
                                                                                          tGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGcc                                                         >  2:1272284/1‑150 (MQ=255)
                                                                                           gggATCATCTGGATG‑GCGGCAGCAAGTTT‑CCCGACGAAGGAAAGCAGGATAGAAAAACTCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCa                                                        <  1:1225153/150‑1 (MQ=255)
                                                                                           gggATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCa                                                        >  1:990190/1‑150 (MQ=255)
                                                                                           gggATCATCTGGAAA‑GCGCCAGCCTTTTT‑CCCGCCGAGGAAAAGCAGGATAGCAAAAACCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTg                                                                   <  2:1158373/139‑1 (MQ=255)
                                                                                           gggATCATCGGGAGC‑GCGCCACCCATTTT‑ACCGCCGAAGAAAAGCAAGATAGAAAAACCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCaa                                                       <  1:1168223/151‑1 (MQ=255)
                                                                                           ggaaccatcGGAATG‑GCGGCAGCCAGTTC‑CCCGACGCGGAAAAGCAGGATAGAAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGcc                                                         <  2:834366/144‑1 (MQ=255)
                                                                                            gtagaatcTGGAAG‑GCGGAAGCCAGTTA‑ACAGAGAAAGAAAAACGAGAAAGCGAAAACCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTg                                                                   <  1:60046/133‑1 (MQ=255)
                                                                                                 catctggtggggggcggcCAGTTT‑ACCGACGGAGAAAAGGAGTAGAGGGAAACCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCaaa                                                      <  2:931693/135‑1 (MQ=255)
                                                                                                    tGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCgg                                              >  2:751461/1‑151 (MQ=255)
                                                                                                    tGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCgg                                              >  2:1100933/1‑151 (MQ=255)
                                                                                                     ggATG‑GCGGCAGCCTGTTC‑CCCGACGCAGAAAAGCAGAATAGAAAAAACCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCgg                                              <  1:1156800/150‑1 (MQ=255)
                                                                                                     ggATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCgg                                              >  1:806210/1‑150 (MQ=255)
                                                                                                      gATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGa                                             >  2:934243/1‑150 (MQ=255)
                                                                                                        atgccGGCCGCCGGTTT‑CCCGACACAGAAAAGGAGAATAGCGAAATCCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGaa                                            <  1:494698/146‑1 (MQ=255)
                                                                                                        atcgCGGCACCAATTTT‑CCCGGCGAAGAAAAGCAGGATAGAAAAACCCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGaa                                            <  2:710462/147‑1 (MQ=255)
                                                                                                        agggCGGCACCCATTTC‑CCCGCCGCAGAAAAGAGGGATAGAAAAAACCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAaagaag                                                  <  2:1419268/141‑1 (MQ=255)
                                                                                                             ggcggcCATTTT‑CCCGACGCAGAAAAGCGGGATAGAAAAACCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAAt                                          <  2:285096/147‑1 (MQ=255)
                                                                                                              gccacccgtttccccgccgCAGAAAAGAGGGATAGCAAAACCCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCAc                                       <  1:437042/137‑1 (MQ=255)
                                                                                                               cAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGt                                    >  1:92280/1‑151 (MQ=255)
                                                                                                                  ccAGTTT‑CCCGCCGCAGAAAAGCAGAATAGAAAAAACCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTc                                   <  2:841585/149‑1 (MQ=255)
                                                                                                                  ccAGTTT‑CCCGACGCAGAAAAGCAGGATAGCAAAAACCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGt                                    <  2:720528/148‑1 (MQ=255)
                                                                                                                    catttccccGACGAAGAAAAGCAGAATAGAAAAACCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGAc                                <  2:92289/144‑1 (MQ=255)
                                                                                                                      ttt‑ACCGACGAAGAAAAGCAGGATAGCAAAAATCCCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACaa                              <  2:720614/150‑1 (MQ=255)
                                                                                                                             gACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCAtt                        >  1:1385687/1‑150 (MQ=255)
                                                                                                                                     aaaGCGGGATAGCAAAACCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTa                       <  1:729966/143‑1 (MQ=255)
                                                                                                                                       aGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCg             >  2:666553/1‑151 (MQ=255)
                                                                                                                                       aGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCg             >  1:394775/1‑151 (MQ=255)
                                                                                                                                       aGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACACCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCg             >  2:474681/1‑151 (MQ=255)
                                                                                                                                           ggagagaaaaaaCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGaa           <  2:578730/141‑1 (MQ=255)
                                                                                                                                           ggaaaacaaaaaCCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACg         <  2:533712/143‑1 (MQ=255)
                                                                                                                                            gaaggaaaaaaCCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCg             <  2:984765/139‑1 (MQ=255)
                                                                                                                                             atagaaaaaaCCGCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACg         <  2:1557323/143‑1 (MQ=255)
                                                                                                                                                  gAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGc     >  1:1218076/1‑148 (MQ=255)
                                                                                                                                                  gAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCAGGAAATCACGGTCGACAACCCATTAGCAAACATCGAACGGTg      >  1:462072/1‑147 (MQ=255)
                                                                                                                                                  aaaaaCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAAc          <  2:1289054/142‑1 (MQ=255)
                                                                                                                                                  aaaaaCCCCCCCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGa  <  1:381122/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
GTAATCCACTTTCGATTCGATCAAAACACGAATACCGGAAGCACCGATGACACCATAAAGCAGCAGCGAAACGCCGCCCATCACCGGCAATGGGATCATCTGGATA‑GCGGCAGCCAGTTT‑ACCGACGCAGGAAAGCAGGATAGCGAAAATCGCCGCCCCGCCGATAACCCAGGTACTGTAAACACGGGTGATCGCCATCACGCCAATGTTTTCTCCGTAAGTAGTATTTGGCGTAGAGCCAAAGAAGCCGGAAATCACGGTCGACAAGCCATTAGCAAACATCGAACGGTGCAGA  >  NC_000913/2619081‑2619375

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: ‑