Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,725,747 G→T 100% R596S (CGC→AGC)  clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,725,7470GT98.5% 704.1 / 0.7 220R596S (CGC→AGC) clpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (113/103);  total (116/104)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 220 total aligned reads displayed.

GACACCGAGTGTTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCGATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGAAAC  >  CP009273/2725602‑2725891
                                                                                                                                                 |                                                                                                                                                
 acacCGAGTGTTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTAtt                                                                                                                                               >  2:206089/1‑148 (MQ=255)
     cGAGTGTTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGa                                                                                                                                           <  1:2957657/148‑1 (MQ=255)
      gAGTGTTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGAt                                                                                                                                          <  1:2517322/148‑1 (MQ=255)
          gtTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGc                                                                                                                                      <  2:3252247/148‑1 (MQ=255)
          gtTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGc                                                                                                                                      <  1:969837/148‑1 (MQ=255)
            ttCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCa                                                                                                                                    >  1:2163544/1‑148 (MQ=255)
               tcCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGcc                                                                                                                                 <  2:625313/148‑1 (MQ=255)
                 cATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGcccc                                                                                                                               >  1:1126622/1‑148 (MQ=255)
                    aaaCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCa                                                                                                                            >  1:753699/1‑148 (MQ=255)
                    aaaCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCa                                                                                                                            >  1:735534/1‑148 (MQ=255)
                    aaaCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCa                                                                                                                            >  1:3160915/1‑148 (MQ=255)
                    aaaCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCa                                                                                                                            >  1:2319965/1‑148 (MQ=255)
                    aaaCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCa                                                                                                                            <  1:172211/148‑1 (MQ=255)
                      aCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACg                                                                                                                          >  2:95009/1‑148 (MQ=255)
                      aCTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACg                                                                                                                          >  2:29416/1‑148 (MQ=255)
                        tctgACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGc                                                                                                                        <  2:2882069/145‑1 (MQ=255)
                         cGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCt                                                                                                                       >  1:1574495/1‑148 (MQ=255)
                                atCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTAcgacga                                                                                                                >  1:1073011/1‑148 (MQ=255)
                                atCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTAcgacga                                                                                                                >  2:3175912/1‑148 (MQ=255)
                                atCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTAcgacga                                                                                                                <  1:1549378/148‑1 (MQ=255)
                                 tCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGaa                                                                                                               >  2:1063778/1‑148 (MQ=255)
                                 tCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGaa                                                                                                               >  2:1883846/1‑148 (MQ=255)
                                  cGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACCACGAAt                                                                                                              >  2:1782546/1‑148 (MQ=255)
                                    aTACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAg                                                                                                            <  2:568956/148‑1 (MQ=255)
                                      aCGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCg                                                                                                          <  2:741411/148‑1 (MQ=255)
                                      aCGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCg                                                                                                          <  1:2674299/148‑1 (MQ=255)
                                       cGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGt                                                                                                         <  1:2650538/148‑1 (MQ=255)
                                        ggACCATCGCCTCGTCGCTATCAAACATACAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGtt                                                                                                        <  2:2563602/148‑1 (MQ=255)
                                          aCCATCGCCTCGTCGCTATCAAACCTAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAg                                                                                                      <  1:78735/148‑1 (MQ=255)
                                          aCCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAg                                                                                                      >  2:2568374/1‑148 (MQ=255)
                                          aCCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAg                                                                                                      <  2:2849973/148‑1 (MQ=255)
                                          aCCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAg                                                                                                      <  1:3401951/148‑1 (MQ=255)
                                          aCCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAg                                                                                                      >  2:808560/1‑148 (MQ=255)
                                            cATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCTCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGAt                                                                                                    <  2:3141954/148‑1 (MQ=255)
                                            cATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGAt                                                                                                    <  1:338757/148‑1 (MQ=255)
                                             aTCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATa                                                                                                   <  2:2463933/148‑1 (MQ=255)
                                                gCCTCGTCGCTATCAAACAGAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCg                                                                                                <  2:154268/148‑1 (MQ=255)
                                                  cTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCa                                                                                              >  1:656461/1‑148 (MQ=255)
                                                   tcgtcgCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCAt                                                                                             >  1:2594885/1‑148 (MQ=255)
                                                   tcgtcgCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCAt                                                                                             >  1:3097353/1‑148 (MQ=255)
                                                   tcgtcgCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCAt                                                                                             <  2:2311622/148‑1 (MQ=255)
                                                   tcgtcgCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCAt                                                                                             <  1:654267/148‑1 (MQ=255)
                                                      tcgCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCaa                                                                                          <  1:801703/148‑1 (MQ=255)
                                                       cgCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAAc                                                                                         >  1:1924458/1‑148 (MQ=255)
                                                        gCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAAcc                                                                                        >  1:1315303/1‑148 (MQ=255)
                                                          tATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGc                                                                                      <  1:1047822/148‑1 (MQ=255)
                                                             cAAACCTAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTc                                                                                   <  1:579442/148‑1 (MQ=255)
                                                             cAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTc                                                                                   >  2:706871/1‑148 (MQ=255)
                                                              aaaCATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCg                                                                                  <  1:616935/148‑1 (MQ=255)
                                                                    aaaGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGa                                                                            >  2:2817334/1‑148 (MQ=255)
                                                                        ttCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCaa                                                                        >  2:1572561/1‑148 (MQ=255)
                                                                           gCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTa                                                                     <  2:192366/148‑1 (MQ=255)
                                                                             cAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTAcg                                                                   >  1:1643923/1‑148 (MQ=255)
                                                                             cAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTAcg                                                                   <  2:261449/148‑1 (MQ=255)
                                                                              aGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTAcgc                                                                  >  1:372633/1‑148 (MQ=255)
                                                                                  ccTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATg                                                              <  1:690123/148‑1 (MQ=255)
                                                                                  ccTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATg                                                              >  2:1828680/1‑148 (MQ=255)
                                                                                  ccTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATg                                                              >  1:128584/1‑148 (MQ=255)
                                                                                   cTTACAAAGCTCTGTTTTCCCCAGACCAGTTGGGCCGAGGAAAAGGAATGAACCAATCGGTCTATTTGGATGCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCTACCGCTTCGTTCTGACCAATTTCGCGATgg                                                             <  1:494901/148‑1 (MQ=255)
                                                                                     tACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTg                                                           <  2:2163544/148‑1 (MQ=255)
                                                                                       cAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCa                                                         >  2:2747659/1‑148 (MQ=255)
                                                                                        aaaGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAg                                                        >  2:2871418/1‑148 (MQ=255)
                                                                                          aGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGtt                                                      >  2:2803368/1‑148 (MQ=255)
                                                                                          aGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGtt                                                      <  2:2261306/148‑1 (MQ=255)
                                                                                          aGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGtt                                                      >  1:1686503/1‑148 (MQ=255)
                                                                                           gCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTc                                                     >  1:227118/1‑148 (MQ=255)
                                                                                            ctctGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCt                                                    >  2:1135338/1‑148 (MQ=255)
                                                                                             tctGTTTTCCCCACACCAGTTGGGACGAGGAACAGGAATGAACCAATCGGGCGATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGAAAAATTACGCGATGGTGCAGTTCtt                                                   >  1:2413259/1‑148 (MQ=255)
                                                                                              ctGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTg                                                  <  2:128584/148‑1 (MQ=255)
                                                                                                    tCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCAt                                            <  2:1344353/148‑1 (MQ=255)
                                                                                                     ccccACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATa                                           >  2:3458747/1‑148 (MQ=255)
                                                                                                            ccAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACgcagca                                    <  2:872486/148‑1 (MQ=255)
                                                                                                             cAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGcagcag                                   <  1:2568374/148‑1 (MQ=255)
                                                                                                              aGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGt                                  >  1:1456435/1‑148 (MQ=255)
                                                                                                               gTTGGGCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGtt                                 <  2:2502256/148‑1 (MQ=255)
                                                                                                                   ggCCGAGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTc                             >  1:1998031/1‑148 (MQ=255)
                                                                                                                       gAGGCACAGGAATGAACCAATCGTGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAAGAGCGTTAGATACCGCAACAACCGCTTCGTTCTGACCAATTCCGCGAGGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCgcgc                         <  2:1996785/148‑1 (MQ=255)
                                                                                                                        aGGAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGcgct                        >  1:1172920/1‑148 (MQ=255)
                                                                                                                         ggAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGATCCATACGCAGCAGTTTTTCGcgctc                       >  2:749989/1‑148 (MQ=255)
                                                                                                                          gAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGcgctcg                      >  1:2001469/1‑148 (MQ=255)
                                                                                                                          gAACAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGcgctcg                      >  1:227798/1‑148 (MQ=255)
                                                                                                                           aaCAGGAATGAACCAATCGGGCGATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGcgctcgc                     >  1:1405069/1‑148 (MQ=255)
                                                                                                                             cAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGTAGCAGTTTTTCGCGCTCGCtt                   >  1:1954565/1‑148 (MQ=255)
                                                                                                                             cAGGAATGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCtt                   <  2:477670/148‑1 (MQ=255)
                                                                                                                                 aaTGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCca               <  1:1869792/148‑1 (MQ=255)
                                                                                                                                   tGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCcatc             >  1:3043596/1‑148 (MQ=255)
                                                                                                                                   tGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCcatc             >  2:2265213/1‑148 (MQ=255)
                                                                                                                                   tGAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTGTCcatc             >  2:2781659/1‑148 (MQ=255)
                                                                                                                                    gAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCcatca            <  2:590288/148‑1 (MQ=255)
                                                                                                                                    gAACCAATCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCcatca            <  2:1867987/148‑1 (MQ=255)
                                                                                                                                         aaTCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCga       >  2:2787805/1‑148 (MQ=255)
                                                                                                                                         aaTCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCga       <  2:1315303/148‑1 (MQ=255)
                                                                                                                                          aTCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCgag      <  2:1007577/148‑1 (MQ=255)
                                                                                                                                           tCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCgaga     <  1:2874929/148‑1 (MQ=255)
                                                                                                                                           tCGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCgaga     >  1:3508827/1‑148 (MQ=255)
                                                                                                                                            cGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGaa    <  2:2656890/148‑1 (MQ=255)
                                                                                                                                            cGGGCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGaa    >  1:819865/1‑148 (MQ=255)
                                                                                                                                             gggCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGaaa   <  2:1175047/148‑1 (MQ=255)
                                                                                                                                              ggCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGAAAc  <  1:1204357/148‑1 (MQ=255)
                                                                                                                                              ggCTATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGAAAc  <  1:1226512/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GACACCGAGTGTTTCTCCATAAACTCGGACATATCGATACGGACCATCGCCTCGTCGCTATCAAACATAAAGTTCGCCAGCGCCTTACAAAGCTCTGTTTTCCCCACACCAGTTGGGCCGAGGAACAGGAATGAACCAATCGGGCGATTTGGATCCGCCAGCCCCGCACGGCTACGACGAATAGCGTTAGATACCGCATCAACCGCTTCGTTCTGACCAATTACGCGATGGTGCAGTTCTTGCTCCATACGCAGCAGTTTTTCGCGCTCGCTTTCCATCATGCGAGAAAC  >  CP009273/2725602‑2725891

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: