Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,376,178 C→A 100% E1004* (GAA→TAA)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,376,1780CA99.0% 1041.1 / ‑1.1 319E1004* (GAA→TAA) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (181/134);  total (184/135)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 323 total aligned reads displayed.

TACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGCG  >  CP009273/4376031‑4376323
                                                                                                                                                   |                                                                                                                                                 
tACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGACGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGc                                                                                                                                                 >  2:2377574/1‑150 (MQ=255)
tACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGc                                                                                                                                                 <  1:2481798/150‑1 (MQ=255)
 aCGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCt                                                                                                                                                <  1:6988/150‑1 (MQ=255)
 aCGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCt                                                                                                                                                <  1:3506193/150‑1 (MQ=255)
  cGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTa                                                                                                                                               <  1:308620/150‑1 (MQ=255)
      gACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTgg                                                                                                                                           <  1:4109401/150‑1 (MQ=255)
      gACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTgg                                                                                                                                           >  1:3271627/1‑150 (MQ=255)
       aCCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGc                                                                                                                                          >  2:2684237/1‑150 (MQ=255)
            tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGTAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCgg                                                                                                                                     <  2:4288844/150‑1 (MQ=255)
            tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCgg                                                                                                                                     <  2:2531629/150‑1 (MQ=255)
            tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCgg                                                                                                                                     <  2:2961192/150‑1 (MQ=255)
            tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCgg                                                                                                                                     >  1:2844837/1‑150 (MQ=255)
             ctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGc                                                                                                                                    >  1:1114499/1‑150 (MQ=255)
             ctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGc                                                                                                                                    <  1:4233293/150‑1 (MQ=255)
             ctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGc                                                                                                                                    <  2:2976843/150‑1 (MQ=255)
                  cggcgTAAATACGCAGCTCGAAAATCTGGATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAgggg                                                                                                                               >  2:4177814/1‑150 (MQ=255)
                  cggcgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAgggg                                                                                                                               >  1:2159957/1‑150 (MQ=255)
                     cgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCg                                                                                                                            <  2:2199919/150‑1 (MQ=255)
                       tAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGt                                                                                                                          >  1:1857779/1‑150 (MQ=255)
                         aaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                        <  2:4387184/150‑1 (MQ=255)
                          aTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATc                                                                                                                       >  1:2360397/1‑150 (MQ=255)
                          aTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATGTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATc                                                                                                                       >  2:308181/1‑150 (MQ=255)
                             cGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTc                                                                                                                    <  1:906389/150‑1 (MQ=255)
                              gCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCa                                                                                                                   <  1:2923071/150‑1 (MQ=255)
                                aGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAac                                                                                                                 >  2:4772791/1‑150 (MQ=255)
                                  cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                               <  2:5128430/150‑1 (MQ=255)
                                  cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                               >  2:2803464/1‑150 (MQ=255)
                                  cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                               >  1:3210914/1‑150 (MQ=255)
                                  cTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                               >  2:1802442/1‑150 (MQ=255)
                                   tCGAAAATCTGACTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcacc                                                                                                              <  1:4151552/150‑1 (MQ=255)
                                    cGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcacca                                                                                                             <  2:4314470/150‑1 (MQ=255)
                                    cGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcacca                                                                                                             >  2:2307260/1‑150 (MQ=255)
                                       aaaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgc                                                                                                          >  1:4797613/1‑150 (MQ=255)
                                            tGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                     <  1:2340210/150‑1 (MQ=255)
                                             gAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGa                                                                                                    >  1:2733466/1‑150 (MQ=255)
                                               aTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                  >  2:2471666/1‑150 (MQ=255)
                                               aTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                  >  2:3645932/1‑150 (MQ=255)
                                               aTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGGTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                  <  2:4664129/150‑1 (MQ=255)
                                               aTCCCCTGTGGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                                  <  2:1943742/150‑1 (MQ=255)
                                                tCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCg                                                                                                 <  2:1527566/150‑1 (MQ=255)
                                                 ccccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGa                                                                                                >  2:170810/1‑150 (MQ=255)
                                                  cccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAg                                                                                               >  2:3606753/1‑150 (MQ=255)
                                                     tGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCa                                                                                            >  2:1386327/1‑150 (MQ=255)
                                                      gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                           <  2:2950066/150‑1 (MQ=255)
                                                      gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                           <  2:808749/150‑1 (MQ=255)
                                                       ttGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCaga                                                                                          >  1:2706189/1‑150 (MQ=255)
                                                         gCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGCGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCCCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTCGCGGCACCACCCCGCGCGTGACCGAGTCagagc                                                                                        >  2:2277582/1‑149 (MQ=255)
                                                           aGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGc                                                                                      >  2:3166587/1‑150 (MQ=255)
                                                            gATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCg                                                                                     >  1:4527362/1‑150 (MQ=255)
                                                              tCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCTGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                   <  2:11767/150‑1 (MQ=255)
                                                               ccaccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATATCCGTCACTTCTTAGCTATTCGCATCGGCAGGGGCCGGTCTCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAcc                                                                                  >  1:5054235/1‑150 (MQ=255)
                                                                caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCa                                                                                 <  1:224687/150‑1 (MQ=255)
                                                                caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGGGACCGAGTCAGAGAGCGACCa                                                                                 <  2:1108929/150‑1 (MQ=255)
                                                                 accaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAg                                                                                >  1:1046655/1‑150 (MQ=255)
                                                                 accaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAg                                                                                >  1:3627178/1‑150 (MQ=255)
                                                                   caGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                              >  2:3851093/1‑150 (MQ=255)
                                                                   caGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                              >  2:1309989/1‑150 (MQ=255)
                                                                   caGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                              <  1:3703771/150‑1 (MQ=255)
                                                                      gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTggt                                                                           <  1:3016113/150‑3 (MQ=255)
                                                                      gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                           <  2:4573303/150‑1 (MQ=255)
                                                                       aaGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATa                                                                          >  2:3359288/1‑150 (MQ=255)
                                                                       aaGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATa                                                                          >  1:3184379/1‑150 (MQ=255)
                                                                        aGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaa                                                                         <  1:2089307/150‑1 (MQ=255)
                                                                             tCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGc                                                                    <  2:5060591/150‑1 (MQ=255)
                                                                             tCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGc                                                                    >  1:2050985/1‑150 (MQ=255)
                                                                              ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCt                                                                   <  1:3584570/150‑1 (MQ=255)
                                                                              ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCt                                                                   <  1:4890318/150‑1 (MQ=255)
                                                                                ggTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCg                                                                 >  1:718647/1‑150 (MQ=255)
                                                                                 gTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCgg                                                                >  1:2436043/1‑150 (MQ=255)
                                                                                 gTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTACCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCgg                                                                <  2:3010247/150‑1 (MQ=255)
                                                                                      gggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaata                                                           <  2:4356619/150‑1 (MQ=255)
                                                                                         ttGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAAc                                                        <  2:2046013/150‑1 (MQ=255)
                                                                                             cGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCt                                                    >  2:263914/1‑150 (MQ=255)
                                                                                             cGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCt                                                    <  1:2564799/150‑1 (MQ=255)
                                                                                              gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCtt                                                   >  2:48140/1‑150 (MQ=255)
                                                                                               aTCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTg                                                  <  1:2961220/150‑1 (MQ=255)
                                                                                                  aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTc                                               <  1:1480138/150‑1 (MQ=255)
                                                                                                  aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTc                                               <  1:3983863/150‑1 (MQ=255)
                                                                                                   ccAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCg                                              <  1:263751/150‑1 (MQ=255)
                                                                                                     aGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                            >  1:3396054/1‑150 (MQ=255)
                                                                                                     aGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                            <  2:3906296/150‑1 (MQ=255)
                                                                                                        gAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACg                                         >  2:4622207/1‑150 (MQ=255)
                                                                                                        gAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACg                                         <  1:3767278/150‑1 (MQ=255)
                                                                                                          gcagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAt                                       <  1:1685119/150‑1 (MQ=255)
                                                                                                              cGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtc                                   <  2:1577070/150‑1 (MQ=255)
                                                                                                              cGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtc                                   >  1:579426/1‑150 (MQ=255)
                                                                                                                   cgcTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                              <  2:3544626/150‑1 (MQ=255)
                                                                                                                   cgcTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                              >  2:4333035/1‑150 (MQ=255)
                                                                                                                    gcTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTa                             <  1:3000826/150‑1 (MQ=255)
                                                                                                                     cTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTAc                            <  1:1410295/150‑1 (MQ=255)
                                                                                                                      tGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACg                           >  2:2607879/1‑150 (MQ=255)
                                                                                                                        cGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGt                         >  1:3137980/1‑150 (MQ=255)
                                                                                                                           tGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTccc                      >  1:4631345/1‑150 (MQ=255)
                                                                                                                             aGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAg                    <  2:2165991/150‑1 (MQ=255)
                                                                                                                                   aTTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTg              >  1:4127904/1‑150 (MQ=255)
                                                                                                                                        cGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggt         >  1:964597/1‑150 (MQ=255)
                                                                                                                                        cGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggt         >  2:1054349/1‑150 (MQ=255)
                                                                                                                                            aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg     <  2:604579/150‑1 (MQ=255)
                                                                                                                                               tcttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGgcg  <  1:3166587/150‑1 (MQ=255)
                                                                                                                                               tcttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGgcg  <  2:2570790/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGGCG  >  CP009273/4376031‑4376323

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