Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,167,908 T→A 25.5% W80R (TGG→AGG)  nusG → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,167,9080TA25.5% 569.9 / 232.0 282W80R (TGG→AGG) nusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (112/97);  new base A (38/34);  total (150/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 283 total aligned reads displayed.

CAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAA  >  CP009273/4167760‑4168057
                                                                                                                                                    |                                                                                                                                                     
cAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTg                                                                                                                                                      <  2:1585931/150‑1 (MQ=255)
 aaaTTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTgg                                                                                                                                                     <  1:480324/150‑1 (MQ=255)
  aaTTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGc                                                                                                                                                    >  1:4307935/1‑150 (MQ=255)
   aTTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCa                                                                                                                                                   >  2:3889839/1‑150 (MQ=255)
    ttACACAACCTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCAc                                                                                                                                                  <  2:3595858/150‑1 (MQ=255)
      acacAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCt                                                                                                                                                >  1:1263123/1‑150 (MQ=255)
       cacaCCAGGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTg                                                                                                                                               <  1:3290095/150‑1 (MQ=255)
       cacaACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTg                                                                                                                                               >  2:1062931/1‑150 (MQ=255)
       cacaACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTg                                                                                                                                               >  1:2435687/1‑150 (MQ=255)
          aacaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTg                                                                                                                                            <  2:2334443/150‑1 (MQ=255)
          aacaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTg                                                                                                                                            >  1:1834545/1‑150 (MQ=255)
           acaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTgc                                                                                                                                           <  2:3212985/150‑1 (MQ=255)
           acaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTgc                                                                                                                                           <  2:4060885/150‑1 (MQ=255)
           acaTGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTgc                                                                                                                                           <  1:558332/150‑1 (MQ=255)
              tGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCa                                                                                                                                        <  1:3438729/150‑1 (MQ=255)
               ggAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAg                                                                                                                                       <  1:774890/150‑1 (MQ=255)
                gAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGc                                                                                                                                      >  1:729929/1‑150 (MQ=255)
                  aGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGt                                                                                                                                    <  2:1743041/150‑1 (MQ=255)
                   gATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTa                                                                                                                                   >  2:2920244/1‑150 (MQ=255)
                   gATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTa                                                                                                                                   >  2:3129861/1‑150 (MQ=255)
                        gtttgGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACcgcg                                                                                                                              >  2:2108291/1‑150 (MQ=255)
                        gtttgGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACcgcg                                                                                                                              >  1:4169738/1‑150 (MQ=255)
                           tgGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCgtgt                                                                                                                           >  2:1657398/1‑150 (MQ=255)
                               gAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATg                                                                                                                       >  1:572104/1‑150 (MQ=255)
                                 aGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATggg                                                                                                                     >  2:869805/1‑150 (MQ=255)
                                  gTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGc                                                                                                                    <  2:2857409/150‑1 (MQ=255)
                                   tCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCt                                                                                                                   >  2:4299113/1‑150 (MQ=255)
                                      tGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCa                                                                                                                >  1:955436/1‑150 (MQ=255)
                                       ggTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                                                                                                               <  1:3970181/150‑1 (MQ=255)
                                       ggTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCAt                                                                                                               <  2:4359563/150‑1 (MQ=255)
                                          aCCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcgg                                                                                                            <  2:2580836/150‑1 (MQ=255)
                                           ccAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcggc                                                                                                           >  1:1960007/1‑150 (MQ=255)
                                            cAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcggcg                                                                                                          <  2:1971269/150‑1 (MQ=255)
                                             aaCCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcggcgg                                                                                                         >  2:4338641/1‑150 (MQ=255)
                                             aaCCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcggcgg                                                                                                         >  2:135996/1‑150 (MQ=255)
                                             aaCCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATcggcgg                                                                                                         <  1:1700261/150‑1 (MQ=255)
                                               ccGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTa                                                                                                       <  2:2746000/150‑1 (MQ=255)
                                                 gaagaaGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACt                                                                                                     >  1:1159594/1‑150 (MQ=255)
                                                 gaagaaGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACt                                                                                                     >  1:2052888/1‑150 (MQ=255)
                                                  aagaagTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACtt                                                                                                    >  2:268632/1‑150 (MQ=255)
                                                   agaagTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTc                                                                                                   >  1:136629/1‑150 (MQ=255)
                                                   agaagTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACCTc                                                                                                   >  2:152117/1‑150 (MQ=255)
                                                    gaagTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTcc                                                                                                  <  2:915736/150‑1 (MQ=255)
                                                      agTGGTTGAAATCCGTGGCGGTCAGCGTCGCAACAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa                                                                                                <  2:2164788/150‑1 (MQ=255)
                                                       gTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGAt                                                                                               <  1:3064484/150‑1 (MQ=255)
                                                       gTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGAt                                                                                               >  2:2602982/1‑150 (MQ=255)
                                                        tGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATc                                                                                              <  2:52608/150‑1 (MQ=255)
                                                         ggTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCg                                                                                             <  2:1283926/150‑1 (MQ=255)
                                                          gTTGAAATCCGTGGCGGTCAGCGTCGCAAACGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGt                                                                                            <  2:3720837/150‑1 (MQ=255)
                                                            tGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTcc                                                                                          <  2:847079/150‑1 (MQ=255)
                                                              aaaTCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTg                                                                                        >  2:3403039/1‑150 (MQ=255)
                                                               aaTCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTgc                                                                                       >  1:2560217/1‑150 (MQ=255)
                                                                aTCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTgcg                                                                                      >  1:3879751/1‑150 (MQ=255)
                                                                   cGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCa                                                                                   <  2:2016482/150‑1 (MQ=255)
                                                                       gCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCa                                                                               <  2:805401/150‑1 (MQ=255)
                                                                        cGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAg                                                                              >  2:3416519/1‑150 (MQ=255)
                                                                        cGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAg                                                                              <  1:3132277/150‑1 (MQ=255)
                                                                         ggTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGc                                                                             <  1:3910926/150‑1 (MQ=255)
                                                                          gTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCg                                                                            >  1:1580166/1‑150 (MQ=255)
                                                                            cAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGAt                                                                          <  1:4006805/150‑1 (MQ=255)
                                                                            cAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGAt                                                                          >  1:3120591/1‑150 (MQ=255)
                                                                              gCGTCGCAAAAGCGAACGTAAAGTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATaa                                                                        <  2:2652702/150‑1 (MQ=255)
                                                                                  cGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGagaaagaa                                                                    >  2:1851371/1‑150 (MQ=255)
                                                                                  cGCAAAAGCGAACGTAAATTCGTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                    <  2:3280839/150‑1 (MQ=255)
                                                                                   gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaag                                                                   >  1:2756165/1‑150 (MQ=255)
                                                                                    cAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGt                                                                  <  1:612783/150‑1 (MQ=255)
                                                                                       aaGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGa                                                               <  2:1254219/150‑1 (MQ=255)
                                                                                          cGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGc                                                            >  2:1090127/1‑150 (MQ=255)
                                                                                             aCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAt                                                         <  2:2282219/150‑1 (MQ=255)
                                                                                              cGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATACAGAAGTCGATGCGAtt                                                        >  2:2770577/1‑150 (MQ=255)
                                                                                              cGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAtt                                                        >  2:58866/1‑150 (MQ=255)
                                                                                              cGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAtt                                                        >  1:793389/1‑150 (MQ=255)
                                                                                               gTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTa                                                       <  1:3007329/150‑1 (MQ=255)
                                                                                                tAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCgatgat                                                      >  1:4258714/1‑150 (MQ=255)
                                                                                                 aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGTCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATg                                                     <  1:2285608/150‑1 (MQ=255)
                                                                                                 aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATg                                                     <  1:112413/150‑1 (MQ=255)
                                                                                                 aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATg                                                     >  2:1514545/1‑150 (MQ=255)
                                                                                                  aaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGa                                                    <  2:671518/150‑1 (MQ=255)
                                                                                                  aaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGa                                                    <  2:1114682/150‑1 (MQ=255)
                                                                                                     tcttcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAAcc                                                 >  2:2985066/1‑150 (MQ=255)
                                                                                                         cccTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCt                                             >  1:83858/1‑150 (MQ=255)
                                                                                                                tACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                      >  1:2847754/1‑150 (MQ=255)
                                                                                                                 aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAgg                                     <  1:3827409/150‑1 (MQ=255)
                                                                                                                  cGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGt                                    <  1:2879345/150‑1 (MQ=255)
                                                                                                                    tCCTCGTTCAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTg                                  >  2:3066876/1‑150 (MQ=255)
                                                                                                                         gTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGAt                             <  1:1327340/150‑1 (MQ=255)
                                                                                                                          ttCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATa                            >  1:3684156/1‑150 (MQ=255)
                                                                                                                          ttCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATa                            >  1:1493895/1‑150 (MQ=255)
                                                                                                                            cAGATGGTGATGAACGACGCGAGCTGGCACCTGGGGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAg                          <  2:3076652/150‑1 (MQ=255)
                                                                                                                            cAGATGGTGATGAACGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAg                          >  2:3522878/1‑150 (MQ=255)
                                                                                                                              gATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGcc                        <  1:1053661/150‑1 (MQ=255)
                                                                                                                               agggTGATGAACGACGCGAGCGGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCGTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCcg                       <  2:2435687/148‑1 (MQ=255)
                                                                                                                                        aaCGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGaaa              >  1:2827288/1‑150 (MQ=255)
                                                                                                                                        aaCGACGCGAGCAGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGaaa              >  2:408649/1‑150 (MQ=255)
                                                                                                                                             cgcgAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCt         >  2:2401658/1‑150 (MQ=255)
                                                                                                                                                 aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGttt     >  1:1554652/1‑150 (MQ=255)
                                                                                                                                                 aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGttt     <  2:3480451/150‑1 (MQ=255)
                                                                                                                                                 aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGttt     >  1:1006673/1‑150 (MQ=255)
                                                                                                                                                 aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGttt     >  2:370831/1‑150 (MQ=255)
                                                                                                                                                   cTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGa   >  1:4037473/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAATTACACAACATGGAAGATTTGTTTGGTGAAGTCATGGTACCAACCGAAGAAGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAA  >  CP009273/4167760‑4168057

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: