Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,318,517 A→T 100% T522T (ACA→ACT crfC → clamp‑binding sister replication fork colocalization protein, dynamin‑related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,318,5170AT99.6% 590.5 / ‑6.0 180T522T (ACA→ACTcrfCclamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (82/97);  total (83/97)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
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 182 total aligned reads displayed.

GCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTT  >  CP009273/4318372‑4318661
                                                                                                                                                 |                                                                                                                                                
gcgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGActct                                                                                                                                                >  2:3105866/1‑148 (MQ=255)
 cgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTg                                                                                                                                               <  1:2085397/148‑1 (MQ=255)
     aTTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTAc                                                                                                                                           <  2:2062415/148‑1 (MQ=255)
      ttAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACa                                                                                                                                          >  1:554336/1‑148 (MQ=255)
       tAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAg                                                                                                                                         >  1:859856/1‑148 (MQ=255)
           cATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACgg                                                                                                                                     <  2:160223/148‑1 (MQ=255)
              gAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGgcgc                                                                                                                                  >  1:3099088/1‑148 (MQ=255)
                    gAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAaca                                                                                                                            <  2:1786704/148‑1 (MQ=255)
                      gCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACAcc                                                                                                                          >  2:147957/1‑148 (MQ=255)
                         ggCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATc                                                                                                                       >  1:2955448/1‑148 (MQ=255)
                          gCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCa                                                                                                                      >  2:2169372/1‑148 (MQ=255)
                           cccTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAg                                                                                                                     >  2:3250133/1‑148 (MQ=255)
                            ccTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGt                                                                                                                    <  1:840790/148‑1 (MQ=255)
                                 ccTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGa                                                                                                               <  2:195984/148‑1 (MQ=255)
                                   tCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACt                                                                                                             >  1:1272866/1‑148 (MQ=255)
                                    ccGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACtt                                                                                                            >  2:3194250/1‑148 (MQ=255)
                                    ccGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACtt                                                                                                            >  1:2137827/1‑148 (MQ=255)
                                      gCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTa                                                                                                          <  2:3042506/148‑1 (MQ=255)
                                        cAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACg                                                                                                        <  2:2335257/148‑1 (MQ=255)
                                           ccACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTg                                                                                                     <  1:717865/148‑1 (MQ=255)
                                           ccACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTg                                                                                                     <  1:3228673/148‑1 (MQ=255)
                                               tttCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGtt                                                                                                 >  2:2223999/1‑148 (MQ=255)
                                                ttCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGttt                                                                                                <  2:142310/148‑1 (MQ=255)
                                                     cGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCg                                                                                           <  1:1922463/148‑1 (MQ=255)
                                                     cGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCg                                                                                           <  1:478353/148‑1 (MQ=255)
                                                       tCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAg                                                                                         <  1:2440434/148‑1 (MQ=255)
                                                        cAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGc                                                                                        <  1:416670/148‑1 (MQ=255)
                                                        cAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGc                                                                                        <  1:1427923/148‑1 (MQ=255)
                                                         aaCAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCt                                                                                       <  1:2047130/148‑1 (MQ=255)
                                                          aCAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCtt                                                                                      >  1:500889/1‑148 (MQ=255)
                                                           cAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTg                                                                                     >  1:3103927/1‑148 (MQ=255)
                                                            aGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGa                                                                                    <  2:2763873/148‑1 (MQ=255)
                                                            aGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTAAAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGa                                                                                    <  1:3198980/148‑1 (MQ=255)
                                                             gCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGaa                                                                                   >  2:1674962/1‑148 (MQ=255)
                                                              cAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAAt                                                                                  <  2:136226/148‑1 (MQ=255)
                                                              cAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAAt                                                                                  >  1:2475530/1‑148 (MQ=255)
                                                               aaGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATc                                                                                 <  1:1427554/148‑1 (MQ=255)
                                                               aaGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATc                                                                                 <  1:2609385/148‑1 (MQ=255)
                                                               aaGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATc                                                                                 <  1:2608063/148‑1 (MQ=255)
                                                               aaGATGCGCTGAAGGTGCAGTTAGCCGCCGCCCCCCACGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATc                                                                                 <  2:1629166/148‑1 (MQ=255)
                                                                aGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATcc                                                                                >  2:1550310/1‑148 (MQ=255)
                                                                aGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATcc                                                                                >  1:1865913/1‑148 (MQ=255)
                                                                aGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATcc                                                                                >  1:1869082/1‑148 (MQ=255)
                                                                aGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATcc                                                                                >  1:1869086/1‑148 (MQ=255)
                                                                  aTGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGATGTGAGACTCTGTTACAGACGGCGCAGAACACCATCCGTCGCGACTTTAAGCTGCGTTTTGCCGAGCTTGAATCCAc                                                                              >  2:2032643/1‑148 (MQ=255)
                                                                    gcgcTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCAccc                                                                            <  2:2981177/148‑1 (MQ=255)
                                                                     cgcTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCt                                                                           <  1:3248580/148‑1 (MQ=255)
                                                                      gcTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCtt                                                                          <  1:2108690/148‑1 (MQ=255)
                                                                       cTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCttt                                                                         <  1:3105866/148‑1 (MQ=255)
                                                                        tGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTg                                                                        >  2:1380917/1‑148 (MQ=255)
                                                                           aGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCg                                                                     <  1:255421/148‑1 (MQ=255)
                                                                            ggTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTAGCCGAGCTTGAATCCACACTTTGCAgg                                                                    >  1:3532145/1‑148 (MQ=255)
                                                                             gTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCagg                                                                   <  2:1100745/148‑1 (MQ=255)
                                                                               gCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGgtg                                                                 <  1:2980279/148‑1 (MQ=255)
                                                                                cAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGgtgt                                                                <  2:694960/148‑1 (MQ=255)
                                                                                    tAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAAc                                                            >  1:2526031/1‑148 (MQ=255)
                                                                                      gccgccTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCg                                                          >  2:3396135/1‑148 (MQ=255)
                                                                                      gccgccTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCg                                                          <  1:3038874/148‑1 (MQ=255)
                                                                                       ccgccTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGa                                                         >  2:1487995/1‑148 (MQ=255)
                                                                                       ccgccTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGa                                                         <  2:247324/148‑1 (MQ=255)
                                                                                        cgccTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGAt                                                        <  2:291836/148‑1 (MQ=255)
                                                                                          ccTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGt                                                      >  2:813364/1‑148 (MQ=255)
                                                                                           cTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGtt                                                     <  2:1869086/148‑1 (MQ=255)
                                                                                           cTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGtt                                                     <  2:1869074/148‑1 (MQ=255)
                                                                                             tGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTAt                                                   <  1:1054856/148‑1 (MQ=255)
                                                                                             tGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTAt                                                   <  1:979250/148‑1 (MQ=255)
                                                                                                ttCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCg                                                <  2:1278614/148‑1 (MQ=255)
                                                                                                 tCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGc                                               <  2:3504440/148‑1 (MQ=255)
                                                                                                 tCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGc                                               <  1:2024670/148‑1 (MQ=255)
                                                                                                  cAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcc                                              >  2:100152/1‑148 (MQ=255)
                                                                                                    ggATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcccc                                            >  2:605922/1‑148 (MQ=255)
                                                                                                     gatgatTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACTCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCa                                           >  1:2361182/1‑148 (MQ=255)
                                                                                                     gatgatTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCa                                           >  1:3348641/1‑148 (MQ=255)
                                                                                                        gatTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                        >  1:1803093/1‑148 (MQ=255)
                                                                                                        gatTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                        <  1:3223509/148‑1 (MQ=255)
                                                                                                        gatTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                        <  2:2384992/148‑1 (MQ=255)
                                                                                                               aGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGTCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGcaaca                                 >  2:2776052/1‑148 (MQ=255)
                                                                                                                gCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGcaacaa                                >  2:1772205/1‑148 (MQ=255)
                                                                                                                gCCGTTAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGcaacaa                                <  1:3170932/148‑1 (MQ=255)
                                                                                                                     tAAGCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCaa                           >  1:1568530/1‑148 (MQ=255)
                                                                                                                        gCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAAt                        >  2:2006332/1‑148 (MQ=255)
                                                                                                                        gCGAGATTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAAt                        <  2:520332/148‑1 (MQ=255)
                                                                                                                           agaTTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGa                     <  2:2663889/148‑1 (MQ=255)
                                                                                                                           agaTTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGa                     >  1:3256324/1‑148 (MQ=255)
                                                                                                                           agaTTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGa                     <  2:173655/148‑1 (MQ=255)
                                                                                                                           agaTTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGa                     <  1:1647908/148‑1 (MQ=255)
                                                                                                                            gaTTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACTCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGaa                    <  2:2361182/148‑1 (MQ=255)
                                                                                                                            gaTTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGaa                    >  2:1729758/1‑148 (MQ=255)
                                                                                                                                cGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGa                <  1:1380917/148‑1 (MQ=255)
                                                                                                                                cGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGa                <  2:500889/148‑1 (MQ=255)
                                                                                                                                  tACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGc              >  1:2629230/1‑148 (MQ=255)
                                                                                                                                   aCCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCg             >  2:914754/1‑148 (MQ=255)
                                                                                                                                     ccGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAg           <  2:2899271/148‑1 (MQ=255)
                                                                                                                                         tgtgAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAg       >  2:961473/1‑148 (MQ=255)
                                                                                                                                           tgAGACTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAggg     >  2:330097/1‑148 (MQ=255)
                                                                                                                                             agaCTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt   <  1:828264/148‑1 (MQ=255)
                                                                                                                                             agaCTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt   <  1:1144759/148‑1 (MQ=255)
                                                                                                                                             agaCTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt   <  1:100152/148‑1 (MQ=255)
                                                                                                                                             agaCTCTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt   >  1:2937517/1‑148 (MQ=255)
                                                                                                                                             agaCTCTGTTACAGAAGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATACAAACTTTGCCGGGTGTTAACAGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt   >  1:2937540/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTT  >  CP009273/4318372‑4318661

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: