Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,069,747 C→T 6.6% Q31Q (CAG→CAA)  rutG ← pyrimidine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,069,7470CT 6.6% 287.1 / 6.3 91Q31Q (CAG→CAA) rutGpyrimidine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (47/38);  new base T (4/2);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGC  >  NC_000913/1069599‑1069894
                                                                                                                                                    |                                                                                                                                                   
ggACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCt                                                                                                                                                    >  1:366011/1‑150 (MQ=255)
   cgcgCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCGATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAg                                                                                                                                                 >  1:644568/1‑150 (MQ=255)
   cgcgCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAg                                                                                                                                                 <  2:1181394/150‑1 (MQ=255)
        ccgccgGTAATAAAGAAAAACAGTAACGTGCCGATCCACGACATTAAAATGGAAATATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAAc                                                                                                                                            >  2:725358/1‑150 (MQ=255)
          gccgGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACgg                                                                                                                                          >  2:1555141/1‑150 (MQ=255)
          gccgGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACgg                                                                                                                                          <  2:97983/150‑1 (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)
                       aaaaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTc                                                                                                                                                     >  2:900178/1‑126 (MQ=255)
                       aaaaaCAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTc                                                                                                                                                     <  1:900178/126‑1 (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:139723/150‑1 (MQ=255)
                                              cGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGt                                                                                                      <  1:100498/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                                                   >  1:358980/1‑150 (MQ=255)
                                                                                                 cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                   <  2:467310/150‑1 (MQ=255)
                                                                                                 cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                   <  2:901623/150‑1 (MQ=255)
                                                                                                 cacCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTg                                                   <  2:433063/150‑1 (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  >  1:1491464/1‑150 (MQ=255)
                                                                                                                                                  gTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGc  <  2:98797/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGACGCGCCCGCCGGTAATAAAGAAAAACAGTAACGTGCCGATCCCCGACATTAAAATGGAAAGATTGGGATCCAGCCCCATCAGAATCGGCATTAACACCGTCGCGCCAAACATCGCCACCGCGTGTTGAACGCCCATTACTGCCGTCTGAGCAAACGGCAATCGTTCATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGC  >  NC_000913/1069599‑1069894

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