Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,069,751 G→T 7.0% A30D (GCT→GAT) ‡ rutG ← pyrimidine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,069,7510GT 7.0% 303.1 / 8.0 92A30D (GCT→GAT) ‡rutGpyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (48/37);  new base T (4/3);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATC  >  NC_000913/1069602‑1069898
                                                                                                                                                     |                                                                                                                                                   
cgcgCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCGATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAg                                                                                                                                                     >  1:644568/1‑150 (MQ=255)
cgcgCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAg                                                                                                                                                     <  2:1181394/150‑1 (MQ=255)
     ccgccgGTAATAAAGAAAAACAGTAACGTGCCGATCCACGACATTAAAATGGAAATATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAAc                                                                                                                                                >  2:725358/1‑150 (MQ=255)
       gccgGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACgg                                                                                                                                              <  2:97983/150‑1 (MQ=255)
       gccgGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACgg                                                                                                                                              >  2:1555141/1‑150 (MQ=255)
        ccgGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTTTGCTCAGACGGc                                                                                                                                             <  1:1360220/150‑1 (MQ=255)
          ggTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCaa                                                                                                                                           >  1:544995/1‑150 (MQ=255)
            taataaAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATc                                                                                                                                         >  1:802273/1‑150 (MQ=255)
              ataaAGAAAAACAGTGACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGt                                                                                                                                       >  2:53586/1‑150 (MQ=255)
               taaAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTTTGCTCAGACGGCagtaatg                                                                                                                                      <  2:1130939/150‑7 (MQ=255)
               taaAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGtt                                                                                                                                      <  2:412967/150‑1 (MQ=255)
                 aaGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCa                                                                                                                                    >  1:98613/1‑150 (MQ=255)
                 aaGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCa                                                                                                                                    >  2:352361/1‑150 (MQ=255)
                  aGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATCGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGAGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCAt                                                                                                                                   >  1:276709/1‑150 (MQ=255)
                   gAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATc                                                                                                                                  >  1:467310/1‑150 (MQ=255)
                     aaaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCg                                                                                                                                >  1:1145813/1‑150 (MQ=255)
                       aaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGgc                                                                                                                              >  1:592018/1‑150 (MQ=255)
                                gTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCAcg                                                                                                                     <  2:161194/150‑1 (MQ=255)
                                gTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCACCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCAcg                                                                                                                     >  2:578641/1‑150 (MQ=255)
                                 tGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCAcgc                                                                                                                    >  1:433063/1‑150 (MQ=255)
                                   ccGATCCCCTACATTAAAATGGAAAGATTTGTATCCATCCCCATCAGAATCTGCTTTAACACCGTCGCGCCAAACATCGGCACCGCGTGTTGAACGCCCATTACTTCCGTCTGATGAAACTGGAATCGTTCATCCTGCTCGACCAcgccg                                                                                                                  <  2:752758/150‑1 (MQ=255)
                                    cGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCAcgccgc                                                                                                                 >  1:1461554/1‑150 (MQ=255)
                                         cccGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGctctct                                                                                                            >  1:323522/1‑150 (MQ=255)
                                         cccGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGctctct                                                                                                            <  2:1439048/150‑1 (MQ=255)
                                         cccGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGctctct                                                                                                            >  2:684048/1‑150 (MQ=255)
                                           cGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGt                                                                                                          <  1:100498/150‑1 (MQ=255)
                                           cGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGt                                                                                                          <  1:139723/150‑1 (MQ=255)
                                              cATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTaga                                                                                                       <  2:592018/150‑1 (MQ=255)
                                               aTTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTagag                                                                                                      >  1:664412/1‑150 (MQ=255)
                                                ttAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAgg                                                                                                     <  2:1459462/150‑1 (MQ=255)
                                                ttAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAgg                                                                                                     >  2:1632163/1‑150 (MQ=255)
                                                    aaTGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGa                                                                                                 >  1:108813/1‑150 (MQ=255)
                                                          aaGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTaa                                                                                           <  2:544995/150‑1 (MQ=255)
                                                           aGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAAc                                                                                          <  1:282040/150‑1 (MQ=255)
                                                            gATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACt                                                                                         >  1:5404/1‑150 (MQ=255)
                                                            gATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACt                                                                                         >  1:98797/1‑150 (MQ=255)
                                                              ttGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGc                                                                                       >  1:996844/1‑150 (MQ=255)
                                                                   aTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTg                                                                                  <  2:515292/150‑1 (MQ=255)
                                                                   aTCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTg                                                                                  <  2:1342360/150‑1 (MQ=255)
                                                                    tCCAGCGCCATCAGAATCGGCATTAACACCGTCGCTCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGa                                                                                 <  2:1376466/150‑1 (MQ=255)
                                                                      cAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAgg                                                                               <  2:323522/150‑1 (MQ=255)
                                                                       aGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGa                                                                              <  2:1587816/150‑1 (MQ=255)
                                                                         ccccATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGaaa                                                                            >  1:967834/1‑150 (MQ=255)
                                                                          cccATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGaaaa                                                                           >  2:688520/1‑150 (MQ=255)
                                                                           ccATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAAc                                                                          >  2:392329/1‑150 (MQ=255)
                                                                            cATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAAcc                                                                         >  1:1571327/1‑150 (MQ=255)
                                                                              tCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGa                                                                       >  2:344607/1‑150 (MQ=255)
                                                                                aGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCTCGTGTTGAACTCCCATTACTGCCGTCTGAGCAAAAGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                     <  2:559652/150‑1 (MQ=255)
                                                                                aGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                     >  2:701070/1‑150 (MQ=255)
                                                                                aGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                     >  1:1562994/1‑150 (MQ=255)
                                                                                aGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAAc                                                                     <  1:293571/150‑1 (MQ=255)
                                                                                 gAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACa                                                                    >  2:1087955/1‑150 (MQ=255)
                                                                                 gAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACa                                                                    >  1:230155/1‑150 (MQ=255)
                                                                                     cGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGc                                                                <  1:1469115/150‑1 (MQ=255)
                                                                                      ggCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCATCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGcc                                                               >  2:1456854/1‑150 (MQ=255)
                                                                                           tAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAg                                                          >  2:426403/1‑150 (MQ=255)
                                                                                              cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                       <  2:433063/150‑1 (MQ=255)
                                                                                              cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                       <  2:467310/150‑1 (MQ=255)
                                                                                              cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                       <  2:901623/150‑1 (MQ=255)
                                                                                              cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                       >  1:358980/1‑150 (MQ=255)
                                                                                                ccGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtc                                                     >  1:1535265/1‑150 (MQ=255)
                                                                                                ccGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtc                                                     <  2:116299/150‑1 (MQ=255)
                                                                                                 cGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGCCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGtct                                                    >  1:1624896/1‑150 (MQ=255)
                                                                                                   tCGCGCCATACATCGCCACCGCGTGTTGAACGCCCTTTACTGCCGTCTGAGCAAACGTCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCTGCTGTCTcc                                                  <  2:276709/150‑1 (MQ=255)
                                                                                                   tCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTcc                                                  <  2:1461554/150‑1 (MQ=255)
                                                                                                    cgcgCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCt                                                 >  2:1154685/1‑150 (MQ=255)
                                                                                                      cgcCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTa                                               <  2:366011/150‑1 (MQ=255)
                                                                                                      cgcCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTa                                               >  1:1495044/1‑150 (MQ=255)
                                                                                                       gcCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTaa                                              <  1:392329/150‑1 (MQ=255)
                                                                                                          aaaCATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAagga                                           <  2:1186439/150‑1 (MQ=255)
                                                                                                            aCATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAaggagg                                         <  1:725358/150‑1 (MQ=255)
                                                                                                               tCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTa                                      <  1:492070/150‑1 (MQ=255)
                                                                                                               tCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTa                                      <  1:53586/150‑1 (MQ=255)
                                                                                                                cGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAACACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTaa                                     >  2:1061413/1‑150 (MQ=255)
                                                                                                                  ccACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACa                                   <  1:1555141/150‑1 (MQ=255)
                                                                                                                          tgtTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGgc                           >  1:375415/1‑150 (MQ=255)
                                                                                                                           gtTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGgcg                          <  1:352361/150‑1 (MQ=255)
                                                                                                                             tGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCa                        <  2:483223/150‑1 (MQ=255)
                                                                                                                             tGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCa                        <  2:664412/150‑1 (MQ=255)
                                                                                                                             tGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCa                        >  2:675169/1‑150 (MQ=255)
                                                                                                                                  gCCCATTACTGCCGTCTGAGCAAACGGCAATCCTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAgc                   >  2:755610/1‑150 (MQ=255)
                                                                                                                                      acgattGCCGTTTGCTCAGACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGt               >  2:1195141/6‑150 (MQ=255)
                                                                                                                                       ttACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTg              <  2:1145813/150‑1 (MQ=255)
                                                                                                                                       cgattGCCGTTTGCTCAGACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTg              >  2:778498/5‑150 (MQ=255)
                                                                                                                                         attGCCGTTTGCTCAGACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGAt            >  1:973970/3‑150 (MQ=255)
                                                                                                                                         aCTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGAt            >  1:632464/1‑150 (MQ=255)
                                                                                                                                            gCCGTTTGCTCAGACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAAc         >  1:863971/1‑150 (MQ=255)
                                                                                                                                            gCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAAc         >  1:1331640/1‑150 (MQ=255)
                                                                                                                                               gTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGc      <  2:98797/150‑1 (MQ=255)
                                                                                                                                               gTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGc      >  1:1491464/1‑150 (MQ=255)
                                                                                                                                                  tGAGCAAACGGCATTCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGAt   <  2:644568/150‑1 (MQ=255)
                                                                                                                                                   gAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATc  >  2:1515054/1‑150 (MQ=255)
                                                                                                                                                   gAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATc  >  2:402679/1‑150 (MQ=255)
                                                                                                                                                   gAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTAGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATc  >  1:111352/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATC  >  NC_000913/1069602‑1069898

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