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% 927.4 / 9.7 269R87H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (134/127);  total (137/132)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 273 total aligned reads displayed.

ATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAC  >  CP009273/3959888‑3960181
                                                                                                                                                    |                                                                                                                                                 
 tCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccacc                                                                                                                                                 >  2:879513/1‑150 (MQ=255)
     gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGttt                                                                                                                                             >  1:4345228/1‑150 (MQ=255)
     gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGttt                                                                                                                                             >  2:670969/1‑150 (MQ=255)
       agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCa                                                                                                                                           >  2:1165818/1‑150 (MQ=255)
       agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCa                                                                                                                                           >  1:911203/1‑150 (MQ=255)
          agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAcc                                                                                                                                        >  2:2022598/1‑150 (MQ=255)
          agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAAcc                                                                                                                                        >  2:1470394/1‑150 (MQ=255)
           gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCt                                                                                                                                       <  1:3342761/150‑1 (MQ=255)
            caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTc                                                                                                                                      <  1:4405988/150‑1 (MQ=255)
            caCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTc                                                                                                                                      >  2:3296540/1‑150 (MQ=255)
               gCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGc                                                                                                                                   >  2:2745471/1‑150 (MQ=255)
                 aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGATATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                 >  2:3030481/1‑150 (MQ=255)
                 aaaGAGTGGCGAAGATATCTTTGGTGATGGCCTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCAc                                                                                                                                 >  2:3670922/1‑150 (MQ=255)
                  aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACt                                                                                                                                >  2:4499294/1‑150 (MQ=255)
                  aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGATTCAACCTCCGCACt                                                                                                                                >  2:2260402/1‑150 (MQ=255)
                     agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCTAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                             <  1:2007135/150‑1 (MQ=255)
                     agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGt                                                                                                                             >  2:4307491/1‑150 (MQ=255)
                      gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTg                                                                                                                            >  2:1362675/1‑150 (MQ=255)
                      gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTg                                                                                                                            >  1:2364766/1‑150 (MQ=255)
                       tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGa                                                                                                                           <  1:613950/150‑1 (MQ=255)
                       tGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGa                                                                                                                           <  1:2894484/150‑1 (MQ=255)
                         gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATa                                                                                                                         <  1:1023342/150‑1 (MQ=255)
                         gCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTACTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTATCAACCTCCGCACTGGTGATa                                                                                                                         >  2:4130782/1‑150 (MQ=255)
                            aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                      <  1:3190886/150‑1 (MQ=255)
                              gATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                    <  1:2530822/150‑1 (MQ=255)
                                tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCAtctc                                                                                                                  <  2:3212458/150‑1 (MQ=255)
                                  tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTg                                                                                                                <  1:1642814/150‑1 (MQ=255)
                                   cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                               <  2:2207605/150‑1 (MQ=255)
                                   cTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                               >  2:2169976/1‑150 (MQ=255)
                                     ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                                                                                             <  1:3749908/150‑1 (MQ=255)
                                     ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCACACTGGTGATACCATCTCTGGTa                                                                                                             >  1:1793639/1‑150 (MQ=255)
                                        gTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                                                                                          >  1:1530926/1‑150 (MQ=255)
                                            tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                                                                      >  2:2852101/1‑150 (MQ=255)
                                            tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                                                                      <  1:2100488/150‑1 (MQ=255)
                                            tGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                                                                      >  2:2851988/1‑150 (MQ=255)
                                               cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                                                                                   <  2:1992295/150‑1 (MQ=255)
                                                 tACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                                                                                                 >  1:212968/1‑150 (MQ=255)
                                                   cTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                               <  1:2846301/150‑1 (MQ=255)
                                                     ggAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                             <  2:4345228/150‑1 (MQ=255)
                                                       agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                           >  2:4148030/1‑150 (MQ=255)
                                                        gaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaaga                                                                                          <  2:1908498/150‑1 (MQ=255)
                                                        gaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaaga                                                                                          <  1:1333537/150‑1 (MQ=255)
                                                          tatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                        >  1:2428050/1‑150 (MQ=255)
                                                               cAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGaa                                                                                   >  1:1119323/1‑150 (MQ=255)
                                                                 ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                 <  2:502541/150‑1 (MQ=255)
                                                                 ggatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                                 >  2:3211347/1‑150 (MQ=255)
                                                                    tggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTa                                                                              >  2:384214/1‑150 (MQ=255)
                                                                     ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt                                                                             >  2:773901/1‑150 (MQ=255)
                                                                     ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt                                                                             <  1:1899429/150‑1 (MQ=255)
                                                                     ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt                                                                             >  1:3220264/1‑150 (MQ=255)
                                                                     ggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATATACGTTTCCCCTAGCCAAATCCAACTTTTCAACCTCCGCACTTGTGATAACATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt                                                                             <  2:4091852/150‑1 (MQ=255)
                                                                      gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                            <  1:270945/150‑1 (MQ=255)
                                                                                 ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                                 >  2:1652536/1‑150 (MQ=255)
                                                                                  tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                                >  1:2287754/1‑150 (MQ=255)
                                                                                  tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                                >  1:3538272/1‑150 (MQ=255)
                                                                                     gTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                             <  1:2632583/150‑1 (MQ=255)
                                                                                      ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            <  1:1555169/150‑1 (MQ=255)
                                                                                      ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                            <  1:2776703/150‑1 (MQ=255)
                                                                                         cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                         <  1:4053162/150‑1 (MQ=255)
                                                                                         cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                         >  2:2347308/1‑150 (MQ=255)
                                                                                          gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                        <  2:1550904/150‑1 (MQ=255)
                                                                                          gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                        <  2:1551049/150‑1 (MQ=255)
                                                                                           cAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTGTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                       <  2:1013683/150‑1 (MQ=255)
                                                                                            aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                      >  2:4191700/1‑150 (MQ=255)
                                                                                            aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                      >  2:4191736/1‑150 (MQ=255)
                                                                                            aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                      >  1:1060362/1‑150 (MQ=255)
                                                                                             gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                                                     <  1:199/150‑1 (MQ=255)
                                                                                                   tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                               >  2:2560255/1‑150 (MQ=255)
                                                                                                   tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                               >  2:265088/1‑150 (MQ=255)
                                                                                                    ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                              <  2:956594/150‑1 (MQ=255)
                                                                                                     cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                             >  1:3201369/1‑150 (MQ=255)
                                                                                                     cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                             <  2:485942/150‑1 (MQ=255)
                                                                                                      tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                            <  1:3842617/150‑1 (MQ=255)
                                                                                                        ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                          <  2:3998577/150‑1 (MQ=255)
                                                                                                          tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                        >  2:4286966/1‑150 (MQ=255)
                                                                                                           cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                       <  1:3573206/150‑1 (MQ=255)
                                                                                                             ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                     <  1:3606879/150‑1 (MQ=255)
                                                                                                              cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                    <  1:3475131/150‑1 (MQ=255)
                                                                                                                gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                  <  2:3457411/150‑1 (MQ=255)
                                                                                                                gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                  >  2:2057964/1‑150 (MQ=255)
                                                                                                                gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGAGACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                  <  2:2437180/150‑1 (MQ=255)
                                                                                                                   cTGATGACATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                               >  1:1913402/1‑150 (MQ=255)
                                                                                                                     gatgaCATCTTCGTTTCCCCTAGCCAATTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                             <  2:3423686/150‑1 (MQ=255)
                                                                                                                     gatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                             <  2:3387773/150‑1 (MQ=255)
                                                                                                                     gatgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                             >  2:1037428/1‑150 (MQ=255)
                                                                                                                       tgaCATCTACGTTTCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                           >  1:8065/1‑150 (MQ=255)
                                                                                                                                  tttCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa                <  2:1684314/150‑1 (MQ=255)
                                                                                                                                    tCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAACCCTGAAAACGCCCGCAACAAAATc              >  1:2328460/1‑150 (MQ=255)
                                                                                                                                    tCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc              >  1:2328610/1‑150 (MQ=255)
                                                                                                                                    tCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc              <  2:1690526/150‑1 (MQ=255)
                                                                                                                                    tCCCCTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc              >  2:2933603/1‑150 (MQ=255)
                                                                                                                                     ccccTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc             >  1:446827/1‑150 (MQ=255)
                                                                                                                                        cTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          <  2:3577687/150‑1 (MQ=255)
                                                                                                                                        cTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          <  2:728069/150‑1 (MQ=255)
                                                                                                                                        cTAGCCAAATCCACCGTTTCAACCTCCGCACTGGTGATACAATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct          <  1:1891840/150‑1 (MQ=255)
                                                                                                                                            ccAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga      >  1:3941563/1‑150 (MQ=255)
                                                                                                                                            ccAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga      >  1:2026285/1‑150 (MQ=255)
                                                                                                                                            ccAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga      >  2:2490625/1‑150 (MQ=255)
                                                                                                                                             cAAATCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgag     >  1:3160101/1‑150 (MQ=255)
                                                                                                                                               aaTCCACCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa   <  1:3030481/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAC  >  CP009273/3959888‑3960181

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: ‑