Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,960,036 G→A 92.7% R87H (CGC→CAC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,960,0360GA92.7% 779.2 / 36.8 234R87H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (7/10);  new base A (104/113);  total (111/123)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 241 total aligned reads displayed.

TCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACT  >  CP009273/3959889‑3960182
                                                                                                                                                   |                                                                                                                                                  
tCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccacc                                                                                                                                                  >  2:148586/1‑150 (MQ=255)
tCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccacc                                                                                                                                                  <  1:3206126/150‑1 (MQ=255)
 ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccaccg                                                                                                                                                 >  1:880548/1‑150 (MQ=255)
      agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCa                                                                                                                                            <  2:2894162/150‑1 (MQ=255)
        cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAc                                                                                                                                          <  1:250596/150‑1 (MQ=255)
         agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAcc                                                                                                                                         >  1:1580397/1‑150 (MQ=255)
            aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTcc                                                                                                                                      >  1:2890516/1‑150 (MQ=255)
              gCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGc                                                                                                                                    >  2:2577598/1‑150 (MQ=255)
              gCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGc                                                                                                                                    >  2:3301370/1‑150 (MQ=255)
               cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCa                                                                                                                                   >  2:4840041/1‑150 (MQ=255)
                aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                  >  1:1022431/1‑150 (MQ=255)
                aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                  >  2:1871168/1‑150 (MQ=255)
                    agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                              <  2:3655602/150‑1 (MQ=255)
                    agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                              <  2:1924193/150‑1 (MQ=255)
                      tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTGTCAACCTCCGCACTGGTGa                                                                                                                            >  2:2985165/1‑150 (MQ=255)
                         cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAc                                                                                                                         <  1:3994364/150‑1 (MQ=255)
                          gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAcc                                                                                                                        <  2:1786034/150‑1 (MQ=255)
                          gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAcc                                                                                                                        <  2:397771/150‑1 (MQ=255)
                           aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                       >  1:2937224/1‑150 (MQ=255)
                           aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                       <  1:2196529/150‑1 (MQ=255)
                            aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                      <  1:1887900/150‑1 (MQ=255)
                            aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                      >  2:146156/1‑150 (MQ=255)
                                 tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTg                                                                                                                 <  1:236425/150‑1 (MQ=255)
                                     tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa                                                                                                             <  1:2668248/150‑1 (MQ=255)
                                       gTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                                                                                           >  2:4241173/1‑150 (MQ=255)
                                          aTGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                        <  2:3997397/150‑1 (MQ=255)
                                                 aCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                                                                                                 >  1:3385522/1‑150 (MQ=255)
                                                   tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                                                                                               >  1:3718940/1‑150 (MQ=255)
                                                   tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                                                                                               >  1:17298/1‑150 (MQ=255)
                                                     gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                                                                                             >  2:1961888/1‑150 (MQ=255)
                                                      agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                            <  1:4642165/150‑1 (MQ=255)
                                                         tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                         >  1:4341819/1‑150 (MQ=255)
                                                         tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                         <  1:1908335/150‑1 (MQ=255)
                                                           ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                                                                       <  1:1329551/150‑1 (MQ=255)
                                                           ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                                                                       >  2:239025/1‑150 (MQ=255)
                                                            tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                      >  2:2787309/1‑150 (MQ=255)
                                                             gCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                                                                                     <  1:377157/150‑1 (MQ=255)
                                                              cAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGaa                                                                                    <  2:678940/150‑1 (MQ=255)
                                                                   tggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTa                                                                               <  1:4715193/150‑1 (MQ=255)
                                                                      atTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                            <  1:983257/150‑1 (MQ=255)
                                                                       tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                           <  1:2847047/150‑1 (MQ=255)
                                                                           gTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                       >  1:4137809/1‑150 (MQ=255)
                                                                           gTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                       >  1:4137817/1‑150 (MQ=255)
                                                                           gTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                       >  2:4512505/1‑150 (MQ=255)
                                                                           gTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                       >  1:1305372/1‑150 (MQ=255)
                                                                             ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                                     <  1:3751524/150‑1 (MQ=255)
                                                                                ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAGTTTGCGctgctg                                                                  >  2:3813056/1‑150 (MQ=255)
                                                                                  ccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                                <  1:2178725/150‑1 (MQ=255)
                                                                                   cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                               <  1:1206371/150‑1 (MQ=255)
                                                                                   cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                               <  1:1144433/150‑1 (MQ=255)
                                                                                     ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                             <  2:4115646/150‑1 (MQ=255)
                                                                                     ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                             >  1:2044375/1‑150 (MQ=255)
                                                                                      tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                            <  2:2501976/150‑1 (MQ=255)
                                                                                      tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                            <  2:1074074/150‑1 (MQ=255)
                                                                                      tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                            <  1:4238536/150‑1 (MQ=255)
                                                                                      tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                            >  2:1244867/1‑150 (MQ=255)
                                                                                        cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                          <  2:3495435/150‑1 (MQ=255)
                                                                                             aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                     <  1:2799451/150‑1 (MQ=255)
                                                                                              cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                    >  1:1370868/1‑150 (MQ=255)
                                                                                               aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                   >  1:4542056/1‑150 (MQ=255)
                                                                                                 cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                 >  1:1025116/1‑150 (MQ=255)
                                                                                                 cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                 >  1:475770/1‑150 (MQ=255)
                                                                                                  tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                <  2:4828092/150‑1 (MQ=255)
                                                                                                  tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                <  1:754193/150‑1 (MQ=255)
                                                                                                  tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                >  2:3390007/1‑150 (MQ=255)
                                                                                                   ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                               <  2:2038316/150‑1 (MQ=255)
                                                                                                   ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                               <  2:4654032/150‑1 (MQ=255)
                                                                                                   ccTACCTCACCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                               <  1:573117/150‑1 (MQ=255)
                                                                                                     tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                             <  2:1295780/150‑1 (MQ=255)
                                                                                                     tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                             <  2:1252219/150‑1 (MQ=255)
                                                                                                     tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                             <  2:4531084/150‑1 (MQ=255)
                                                                                                     tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                             <  2:2668592/150‑1 (MQ=255)
                                                                                                         tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                         <  2:2294903/150‑1 (MQ=255)
                                                                                                           gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                       >  1:1523455/1‑150 (MQ=255)
                                                                                                              ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                    <  2:3307125/150‑1 (MQ=255)
                                                                                                              ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGATCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                    >  2:3538020/1‑150 (MQ=255)
                                                                                                               gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                   >  1:4527861/1‑150 (MQ=255)
                                                                                                               gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                   >  2:2427396/1‑150 (MQ=255)
                                                                                                               gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                   <  1:2124973/150‑1 (MQ=255)
                                                                                                                   tgatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                               <  1:3566309/150‑1 (MQ=255)
                                                                                                                   tgatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                               <  1:3566275/150‑1 (MQ=255)
                                                                                                                    gatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                              >  1:2385998/1‑150 (MQ=255)
                                                                                                                          aTCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                        >  1:3791725/1‑150 (MQ=255)
                                                                                                                           tCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                       >  2:4191928/1‑150 (MQ=255)
                                                                                                                              aggTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca                    <  2:4722241/148‑1 (MQ=255)
                                                                                                                              aCGTTTCCCCTCGCCAAATCCACCTTTTCAACCTCCGCAGGGGTTATACCATCTCTGGTATGTTTCGCCCGCCTAAAGATGGTGAACGCTTTTTTGCGCTGCTGAATGTTAACGAATTTAACTTCGACAAACCTGAAAACGCCCGcaaca                    <  2:1021554/150‑1 (MQ=255)
                                                                                                                                  ttCCCCTAGCCAAATCAACCTTTACAACCTCCGCACTGGAGATACCATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTGAAAGCTATAAAGCGCTGCTCAAAGTAAACGACGTTAACTTCGACAAACCTGAAAACGCCCGCACCAAAAt                >  2:1870287/1‑150 (MQ=255)
                                                                                                                                   tCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc               >  1:109409/1‑150 (MQ=255)
                                                                                                                                    ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc              >  1:3007312/1‑150 (MQ=255)
                                                                                                                                    ccccTACCCAAATCCGCCGTTTAAACCTCCGAACGGGTGATACCATCTCTGGTAAGATTCGCCCGCCTAAATATGGTGAACGATATTTTGCGCTGCTTAAATTTAACTAATTTAACTTCGACAACCCTTAAAACGCCCGCAACAAAATcc              <  2:881233/150‑1 (MQ=255)
                                                                                                                                      ccTAGCCAAATCCACCTTTTCAACCTCCGCACTGGTGATACCTTCTCTGTTAAGATTCGCCCGCCTAAAGACGGTGATCGCTATTTTGCGCTGCTGAAATTTAACGATGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc            <  2:490253/150‑1 (MQ=255)
                                                                                                                                      ccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc            >  1:3306472/1‑150 (MQ=255)
                                                                                                                                       cTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct           >  1:3940203/1‑150 (MQ=255)
                                                                                                                                        tAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt          >  1:1779043/1‑150 (MQ=255)
                                                                                                                                        tAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt          <  2:811402/150‑1 (MQ=255)
                                                                                                                                           ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga       <  2:3616396/150‑1 (MQ=255)
                                                                                                                                               aTCCACCTTTTCAACCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc   <  2:2142532/150‑1 (MQ=255)
                                                                                                                                               aTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc   <  2:2675192/150‑1 (MQ=255)
                                                                                                                                               aTCCACCGTTTCAACCGCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc   <  2:578876/150‑1 (MQ=255)
                                                                                                                                                tCCACCGTTTAAACCTCCGCACTGGTGATACCATCTCTGGTATGTTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGTAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACt  <  2:1280522/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACT  >  CP009273/3959889‑3960182

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: