Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 849,632 C→T 15.0% A156T (GCC→ACC)  rhtA ← threonine and homoserine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913849,6320CT15.0% 1951.9 / 264.8 743A156T (GCC→ACC) rhtAthreonine and homoserine efflux system
Reads supporting (aligned to +/- strand):  ref base C (316/315);  new base T (57/55);  total (373/370)
Fisher's exact test for biased strand distribution p-value = 9.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 759 total aligned reads displayed.

CAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCG  >  NC_000913/849484‑849780
                                                                                                                                                    |                                                                                                                                                    
     ccAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCa                                                                                                                                                 <  3:3381043/149‑1 (MQ=255)
     ccAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCa                                                                                                                                                 <  3:5817549/149‑1 (MQ=255)
       agagTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGt                                                                                                                                               <  2:4735065/149‑1 (MQ=255)
       agagTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGt                                                                                                                                               <  1:6962948/149‑1 (MQ=255)
          gTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGcc                                                                                                                                            <  1:4089486/149‑1 (MQ=255)
              ttCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAgcg                                                                                                                                        >  4:665179/1‑149 (MQ=255)
              ttCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAgcg                                                                                                                                        >  2:4043551/1‑149 (MQ=255)
              ttCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAgcg                                                                                                                                        >  1:2519478/1‑149 (MQ=255)
                cACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGca                                                                                                                                      <  1:4791317/149‑1 (MQ=255)
                 aCCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGcac                                                                                                                                     >  3:4974070/1‑149 (MQ=255)
                  ccAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAgc                                                                                                                                         <  4:2639052/144‑1 (MQ=255)
                      ccTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCg                                                                                                                                <  2:2768471/149‑1 (MQ=255)
                       cTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCgg                                                                                                                               <  3:5200091/149‑1 (MQ=255)
                        tGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGt                                                                                                                              <  4:5566807/149‑1 (MQ=255)
                         gAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGtt                                                                                                                             <  3:6594275/149‑1 (MQ=255)
                           aGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCGACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTaa                                                                                                                           >  2:3815687/1‑149 (MQ=255)
                              gcTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATc                                                                                                                        >  4:6105655/1‑149 (MQ=255)
                               cTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGTCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCg                                                                                                                       <  2:5343676/149‑1 (MQ=255)
                                tCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGa                                                                                                                      >  3:6064534/1‑149 (MQ=255)
                                tCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCAAGTGCCAGCGCACAGCCGGTTAAATCGa                                                                                                                      >  4:633620/1‑149 (MQ=255)
                                tCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGtt                                                                                                                             <  1:3602567/142‑1 (MQ=255)
                                tCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGa                                                                                                                      >  2:157822/1‑149 (MQ=255)
                                        gCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGGAGGGCCATGTTCCGCTCCTCCGCGTTGCCCACTAAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCACCGCACAGCCGGTTAAATCGACATGGGaa                                                                                                              >  2:663860/1‑149 (MQ=255)
                                         cACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGaaa                                                                                                             <  4:58183/149‑1 (MQ=255)
                                          aCGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAAc                                                                                                            <  3:2401842/149‑1 (MQ=255)
                                          aCGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAAc                                                                                                            <  2:1011405/149‑1 (MQ=255)
                                            gAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGcccc                                                                                                                                                   <  3:1272145/108‑1 (MQ=255)
                                            gAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCACGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGt                                                                                                          >  4:4880390/1‑149 (MQ=255)
                                             aaaaTTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTc                                                                                                         <  3:7043037/149‑1 (MQ=255)
                                               aaTTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCtt                                                                                                       <  1:6670934/149‑1 (MQ=255)
                                               aaTTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCtt                                                                                                       <  1:3069612/149‑1 (MQ=255)
                                               aaTTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCtt                                                                                                       <  1:1627724/149‑1 (MQ=255)
                                               aaTTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCtt                                                                                                       <  4:3129002/149‑1 (MQ=255)
                                                 ttAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGc                                                                                                     >  3:4324802/1‑149 (MQ=255)
                                                 ttAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGc                                                                                                     <  1:536707/149‑1 (MQ=255)
                                                      gCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAg                                                                                                >  2:3103261/1‑149 (MQ=255)
                                                        tGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCg                                                                                              >  1:4499693/1‑149 (MQ=255)
                                                          cAATCAACGACCCAATTGCCACCGTCGCAGGGCCATGTTCCGCTGCTGCGCGTTGCCCAGTTAAAATGTAAATAGCCCAAGAAGCCCCGGCCCGCAGTGCCAGCGCACAGCCGGTTAAATCGCCATGGGAAACGTCTTGCCCCAGCGGt                                                                                            <  3:2569733/149‑1 (MQ=255)
                                                          cAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGt                                                                                            >  1:1245645/1‑149 (MQ=255)
                                                           aaTCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGTGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTa                                                                                           <  1:1002010/149‑1 (MQ=255)
                                                           aaTCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTa                                                                                           >  3:6895592/1‑149 (MQ=255)
                                                               aaCGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAg                                                                                       >  4:4593300/1‑149 (MQ=255)
                                                                aCGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAATATGTAAATAGCCAAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACAt                                                                                                                   <  4:5237924/120‑1 (MQ=255)
                                                                  gAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGaa                                                                                    >  1:3903896/1‑149 (MQ=255)
                                                                   atcccgttGCCACGGTCGCGGGGCCATGATCCGGTCCTGCGCGTTGCCCGCGTACAATGTAAATAGGCCACCAGGCCCCGGCCCCCATTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGATAGCAGGAAc                                                                                   <  2:6805004/143‑1 (MQ=255)
                                                                   aaCCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAAc                                                                                   <  2:1854417/149‑1 (MQ=255)
                                                                   aaCCAATTGCCACCGTCGCAGAGCCATGTCCCGCTCCTGCGCGTTGCCAACTTAAAATGTAAATAGCCCAACAACCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAAc                                                                                   <  2:4219298/149‑1 (MQ=255)
                                                                      cAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGa                                                                                                                      >  1:5080951/1‑111 (MQ=255)
                                                                      cAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACcac                                                                                >  4:3640979/1‑149 (MQ=255)
                                                                       aaTTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCgg                                                                                                                               <  2:6202762/101‑1 (MQ=255)
                                                                       aaTTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACcaca                                                                               <  2:1992636/149‑1 (MQ=255)
                                                                        aTTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACTAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAg                                                                              >  4:500317/1‑149 (MQ=255)
                                                                        aTTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGccc                                                                                                   >  3:4014743/1‑128 (MQ=255)
                                                                         ttGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGa                                                                             <  1:7285571/149‑1 (MQ=255)
                                                                         ttGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGa                                                                             <  4:2481133/149‑1 (MQ=255)
                                                                             cACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCaa                                                                         <  3:2753273/149‑1 (MQ=255)
                                                                             cACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCaa                                                                         >  2:3284400/1‑149 (MQ=255)
                                                                             cACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCaa                                                                         >  3:36887/1‑149 (MQ=255)
                                                                                     cAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGcc                                                                 <  2:450672/149‑1 (MQ=255)
                                                                                      aGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGAc                                                                                                                     >  3:228872/1‑96 (MQ=255)
                                                                                       gggCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAg                                                               <  2:133640/149‑1 (MQ=255)
                                                                                         gCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGca                                                             <  1:1984458/149‑1 (MQ=255)
                                                                                          ccATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGcac                                                            >  3:946178/1‑149 (MQ=255)
                                                                                            aTGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACaa                                                          >  3:1504612/1‑149 (MQ=255)
                                                                                                 ccGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGc                                                                                                     <  1:4656406/101‑1 (MQ=255)
                                                                                                   gCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGAc                                                   <  3:4593300/149‑1 (MQ=255)
                                                                                                   gCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGAc                                                   >  4:747262/1‑149 (MQ=255)
                                                                                                   gCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGAc                                                   <  1:1377993/149‑1 (MQ=255)
                                                                                                       cTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGaaa                                               <  3:1708458/149‑1 (MQ=255)
                                                                                                        tGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAAt                                              <  3:2120859/149‑1 (MQ=255)
                                                                                                         gcgcGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATc                                             >  2:923993/1‑149 (MQ=255)
                                                                                                            cgTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTAc                                          <  3:6372100/149‑1 (MQ=255)
                                                                                                            cgTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTAc                                          <  1:6123103/149‑1 (MQ=255)
                                                                                                                 cccACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGc                                     >  2:2306016/1‑149 (MQ=255)
                                                                                                                  ccACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATg                                                                                                                  >  4:421903/1‑71 (MQ=255)
                                                                                                                  ccACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCtt                                                                                                       >  1:5975806/1‑82 (MQ=255)
                                                                                                                  ccACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCtt                                                                                                       <  2:5975806/82‑1 (MQ=255)
                                                                                                                   cACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGcga                                   >  4:1713849/1‑149 (MQ=255)
                                                                                                                    aCTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGcgac                                  >  4:2925338/1‑149 (MQ=255)
                                                                                                                       tAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgag                               <  1:2177753/149‑1 (MQ=255)
                                                                                                                         aaaTGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaa                             <  2:357074/149‑1 (MQ=255)
                                                                                                                          aaTGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaag                            <  4:1204536/149‑1 (MQ=255)
                                                                                                                           aTGTAAATAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaaga                           >  4:3409385/1‑149 (MQ=255)
                                                                                                                           aTGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaaga                           >  1:3294478/1‑149 (MQ=255)
                                                                                                                              tAAAAAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaagagaa                        <  4:1511378/149‑1 (MQ=255)
                                                                                                                                aaTAGCCCAACAAGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGaaa                                               >  3:3048796/1‑124 (MQ=255)
                                                                                                                                aaTAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACa                      <  3:5413497/149‑1 (MQ=255)
                                                                                                                                     cccAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGcc                 >  4:2064206/1‑149 (MQ=255)
                                                                                                                                      ccAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCa                >  2:7017877/1‑149 (MQ=255)
                                                                                                                                      ccAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGCAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaagagaa                        >  3:622071/1‑141 (MQ=255)
                                                                                                                                        aacaaGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCgg                                       <  2:6427274/124‑1 (MQ=255)
                                                                                                                                        aacaaGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGATACCGGCGCCAcc              <  4:6999429/149‑1 (MQ=255)
                                                                                                                                         acaaGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgagaagag                          >  1:4382407/1‑136 (MQ=255)
                                                                                                                                          caaGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGc            <  1:6424205/149‑1 (MQ=255)
                                                                                                                                           aaGCCCCGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGcc           >  1:3366209/1‑149 (MQ=255)
                                                                                                                                               cccGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTAc                                          >  3:5646178/1‑114 (MQ=255)
                                                                                                                                               cccGGTCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTAc                                          <  4:2238966/114‑1 (MQ=255)
                                                                                                                                               cccGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTg       <  2:5012689/149‑1 (MQ=255)
                                                                                                                                               cccGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGAAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTg       <  3:4684990/149‑1 (MQ=255)
                                                                                                                                                  ggCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAAACAAAGACCAAGAACCGACAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGCTc    <  4:3939562/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCG  >  NC_000913/849484‑849780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: