Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,960,036 G→A 100% 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,0360GA97.1% 706.0 / 5.4 204R87H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (93/103);  total (97/107)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

CCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACT  >  CP009273/3959890‑3960182
                                                                                                                                                  |                                                                                                                                                  
ccTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccaccg                                                                                                                                                 >  2:2712482/1‑150 (MQ=255)
   gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGttt                                                                                                                                              >  2:986764/1‑150 (MQ=255)
    aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                             >  2:4887280/1‑150 (MQ=255)
    aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTc                                                                                                                                             <  1:3716032/150‑1 (MQ=255)
     agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCa                                                                                                                                            <  1:339946/150‑1 (MQ=255)
       cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAc                                                                                                                                          >  1:3809861/1‑150 (MQ=255)
        agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                                                                                                                         >  1:2746233/1‑150 (MQ=255)
        agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAcc                                                                                                                                         <  1:1718754/150‑1 (MQ=255)
               aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                  <  2:1000335/150‑1 (MQ=255)
                   agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                              >  2:629485/1‑150 (MQ=255)
                    gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTg                                                                                                                             <  2:1171072/150‑1 (MQ=255)
                      ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGAt                                                                                                                           <  1:2465956/150‑1 (MQ=255)
                      ggCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGAt                                                                                                                           >  2:2064062/1‑150 (MQ=255)
                        cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATAc                                                                                                                         >  1:543321/1‑150 (MQ=255)
                           aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                      >  2:4238066/1‑150 (MQ=255)
                            gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                     <  2:4793134/150‑1 (MQ=255)
                                 cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                                >  1:4689516/1‑150 (MQ=255)
                                   ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                              >  1:281289/1‑150 (MQ=255)
                                        gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt                                                                                                         <  1:5136732/150‑1 (MQ=255)
                                         aTGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                        >  1:2188589/1‑150 (MQ=255)
                                         aTGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                        <  1:2725824/150‑1 (MQ=255)
                                         aTGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                        <  1:2805335/150‑1 (MQ=255)
                                          tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                                                                       >  1:3298746/1‑150 (MQ=255)
                                          tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                                                                       <  2:386322/150‑1 (MQ=255)
                                            gCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                     <  2:3771253/150‑1 (MQ=255)
                                             cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                                                                                    >  2:3344301/1‑150 (MQ=255)
                                             cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                                                                                    <  2:578173/150‑1 (MQ=255)
                                              gTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccg                                                                                                   <  1:1557524/150‑1 (MQ=255)
                                                 cTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                                >  1:4426064/1‑150 (MQ=255)
                                                  tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                                                                                               >  2:3624112/1‑150 (MQ=255)
                                                     agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                            <  2:4538468/150‑1 (MQ=255)
                                                     agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                            <  2:1769348/150‑1 (MQ=255)
                                                       atatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaa                                                                                          >  1:4138423/1‑150 (MQ=255)
                                                        tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                         <  2:4202615/150‑1 (MQ=255)
                                                        tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                         <  2:1409572/150‑1 (MQ=255)
                                                         atTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAgg                                                                                        >  1:4549954/1‑150 (MQ=255)
                                                          ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                                                                       >  2:4783067/1‑150 (MQ=255)
                                                           tGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                      <  2:118428/150‑1 (MQ=255)
                                                               ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                  >  2:4678317/1‑150 (MQ=255)
                                                                 atggatTTGGTGTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCt                                                                                <  1:3575315/150‑1 (MQ=255)
                                                                    gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                             >  2:2385699/1‑150 (MQ=255)
                                                                      tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                           <  1:4826494/150‑1 (MQ=255)
                                                                      tttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                           <  1:2132311/150‑1 (MQ=255)
                                                                       ttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                          >  2:944011/1‑150 (MQ=255)
                                                                         ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                        <  2:4145483/150‑1 (MQ=255)
                                                                            ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                                     <  1:3511065/150‑1 (MQ=255)
                                                                             tcctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                    >  2:1259098/1‑150 (MQ=255)
                                                                              cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                   <  1:4520423/150‑1 (MQ=255)
                                                                               ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                                  >  1:72372/1‑150 (MQ=255)
                                                                                tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                                 <  2:4348622/150‑1 (MQ=255)
                                                                                tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                                 <  1:3428124/150‑1 (MQ=255)
                                                                                tccGGTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTGTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                                 <  1:824635/150‑1 (MQ=255)
                                                                                  cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                               >  1:3033295/1‑150 (MQ=255)
                                                                                  cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                               >  2:4097210/1‑150 (MQ=255)
                                                                                   gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                              >  1:1135613/1‑150 (MQ=255)
                                                                                    ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                             <  2:1079219/150‑1 (MQ=255)
                                                                                    ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                             <  1:3545093/150‑1 (MQ=255)
                                                                                     tCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                            <  1:4971485/150‑1 (MQ=255)
                                                                                       cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                          <  2:3242588/150‑1 (MQ=255)
                                                                                        gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                         >  2:1816139/1‑150 (MQ=255)
                                                                                        gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                         <  2:3543923/150‑1 (MQ=255)
                                                                                           gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGACGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                      >  2:4024569/1‑150 (MQ=255)
                                                                                            aCAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                     <  1:2949114/150‑1 (MQ=255)
                                                                                             cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                    >  1:2019262/1‑150 (MQ=255)
                                                                                                 tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTCTTTTGCGCTGCTGAAAGTTAACGACGTTAAc                                                >  2:2932528/1‑150 (MQ=255)
                                                                                                 tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                >  2:1609759/1‑150 (MQ=255)
                                                                                                 tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                >  1:3775836/1‑150 (MQ=255)
                                                                                                  ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                               >  1:802525/1‑150 (MQ=255)
                                                                                                  ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                               >  2:504153/1‑150 (MQ=255)
                                                                                                  ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                               >  2:3702673/1‑150 (MQ=255)
                                                                                                   cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                              <  1:4163294/150‑1 (MQ=255)
                                                                                                   cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                              >  1:1928021/1‑150 (MQ=255)
                                                                                                       cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                          <  2:2529267/150‑1 (MQ=255)
                                                                                                       cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                          <  1:1084705/150‑1 (MQ=255)
                                                                                                             ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcca                                    <  2:2202042/150‑2 (MQ=255)
                                                                                                              gTCCTGATGACATCTACGTTTCCCCTCGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                   <  2:827805/150‑1 (MQ=255)
                                                                                                               tCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGa                                  <  2:282510/150‑1 (MQ=255)
                                                                                                                  tgatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                               >  2:4757338/1‑150 (MQ=255)
                                                                                                                   gatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                              >  1:3835076/1‑150 (MQ=255)
                                                                                                                   gatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                              <  1:1564713/150‑1 (MQ=255)
                                                                                                                    atgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                             >  2:247473/1‑150 (MQ=255)
                                                                                                                       aCATCTACTTTTCCCCTAGCAAAATCCACCGTTTAAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGATGGTTAACGCTATTTTGCGCTGCTGAATGTTAACGATGTTAACTTCGACAAACCTGAAAACGccc                          <  1:4356077/150‑1 (MQ=255)
                                                                                                                        cATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                         <  2:3463557/150‑1 (MQ=255)
                                                                                                                         aTCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                        >  2:644913/1‑150 (MQ=255)
                                                                                                                         aTCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                        >  1:4841638/1‑150 (MQ=255)
                                                                                                                               gTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa                  >  1:223114/1‑150 (MQ=255)
                                                                                                                               gTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa                  >  1:1177657/1‑150 (MQ=255)
                                                                                                                                tttCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa                 >  1:2400050/1‑150 (MQ=255)
                                                                                                                                 ttCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt                <  2:4642673/150‑1 (MQ=255)
                                                                                                                                 ttCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt                >  1:838765/1‑150 (MQ=255)
                                                                                                                                    cccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct             >  2:1577864/1‑150 (MQ=255)
                                                                                                                                    cccTAGCAAAATCCACCGTTTCACCCTCCGCAGTGGTGATACCATCTCTGGTAAGATTCGCCCGCCTAAAGAAGGTGAACGCTATTTTGCGCTGCTGAATGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct             <  1:2087304/150‑1 (MQ=255)
                                                                                                                                       tcgCCAAATCCACCGTTTCACCCTCCGCTCTGTTGATACCATCTCTGGTAAGATTCGCCCGCCGAAATAAGGTGAACGCTATTTTGCGGTGCTGAAATTTCACGATGTTAACTTCGACAAACATGAAAACGCCCGCAACAAAATCCTCtt          <  1:2676014/148‑1 (MQ=255)
                                                                                                                                        aGCCAAATCCCCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGTTAAGATTCGCCCGCCTAAATAAGGTTAACGCTATTTTGCGCTGCTGAAAGTTAACGAATTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt         <  1:776496/150‑1 (MQ=255)
                                                                                                                                        aGCCAAATCCACCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCGTAATGTTACCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt         <  2:2143850/150‑1 (MQ=255)
                                                                                                                                           cAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgag      >  2:156668/1‑150 (MQ=255)
                                                                                                                                            aaaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga     >  2:4637498/1‑150 (MQ=255)
                                                                                                                                            aaaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga     >  2:4637694/1‑150 (MQ=255)
                                                                                                                                             aaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa    >  1:4419206/1‑150 (MQ=255)
                                                                                                                                               tCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACt  >  1:3998596/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACT  >  CP009273/3959890‑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: ‑