Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,318,898 A→C 9.1% E649D (GAA→GAC 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,8980AC 9.1% 577.0 / 50.9 217E649D (GAA→GACcrfCclamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (87/110);  new base C (16/4);  total (103/114)
Fisher's exact test for biased strand distribution p-value = 3.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

AGCAACGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTCTGGCTCTGCGT  >  CP009273/4318751‑4319045
                                                                                                                                                   |                                                                                                                                                   
aGCAACGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGaa                                                                                                                                                     <  1:2550792/148‑1 (MQ=255)
 gCAACGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAg                                                                                                                                                    >  1:128654/1‑148 (MQ=255)
    aCGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCg                                                                                                                                                 >  2:1968656/1‑148 (MQ=255)
     cGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGa                                                                                                                                                >  1:2744299/1‑148 (MQ=255)
     cGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGa                                                                                                                                                >  2:1593489/1‑148 (MQ=255)
          gagcagagcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATctct                                                                                                                                           <  2:956441/148‑1 (MQ=255)
          gagcagagcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATctct                                                                                                                                           <  2:497370/148‑1 (MQ=255)
          gagcagagcaCGCGTTTAGGAGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATctct                                                                                                                                           <  2:2610748/148‑1 (MQ=255)
           agcagagcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTa                                                                                                                                          <  1:2184692/148‑1 (MQ=255)
           agcagagcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTa                                                                                                                                          <  1:2092736/148‑1 (MQ=255)
             cagagcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTAtt                                                                                                                                        >  2:150243/1‑148 (MQ=255)
                agcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGcc                                                                                                                                     >  1:1781847/1‑148 (MQ=255)
                 gcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCca                                                                                                                                    <  1:875803/148‑1 (MQ=255)
                 gcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCca                                                                                                                                    <  1:2681656/148‑1 (MQ=255)
                 gcaCGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCca                                                                                                                                    <  1:2531613/148‑1 (MQ=255)
                    cgcgTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCcaaca                                                                                                                                 >  1:323976/1‑148 (MQ=255)
                     gcgTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCcaacaa                                                                                                                                <  1:121476/148‑1 (MQ=255)
                     gcgTTTAGGCGTTGTGCGTGAGACTTTGTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCcaacaa                                                                                                                                <  2:3239211/148‑1 (MQ=255)
                      cgTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACaaa                                                                                                                               >  2:1107945/1‑148 (MQ=255)
                         ttAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATgg                                                                                                                            <  1:69243/148‑1 (MQ=255)
                         ttAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATgg                                                                                                                            <  1:3085390/148‑1 (MQ=255)
                           aGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCt                                                                                                                          >  2:551628/1‑148 (MQ=255)
                             gCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTaa                                                                                                                        <  2:1130468/148‑1 (MQ=255)
                             gCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTaa                                                                                                                        <  1:2638728/148‑1 (MQ=255)
                               gTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAg                                                                                                                      <  2:818437/148‑1 (MQ=255)
                                ttGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGt                                                                                                                     >  1:1079419/1‑148 (MQ=255)
                                ttGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGt                                                                                                                     <  2:1718245/148‑1 (MQ=255)
                                ttGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGt                                                                                                                     <  1:2249016/148‑1 (MQ=255)
                                ttGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGt                                                                                                                     <  2:42854/148‑1 (MQ=255)
                                  gtgCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTc                                                                                                                   <  1:792416/148‑1 (MQ=255)
                                    gCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTg                                                                                                                 >  1:1411287/1‑148 (MQ=255)
                                       tGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCa                                                                                                              >  2:2524061/1‑148 (MQ=255)
                                          gaCTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCagag                                                                                                           >  2:791367/1‑148 (MQ=255)
                                           aCTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGc                                                                                                          <  2:2703483/148‑1 (MQ=255)
                                           aCTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGc                                                                                                          >  2:793723/1‑148 (MQ=255)
                                           aCTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGc                                                                                                          >  1:454178/1‑148 (MQ=255)
                                           aCTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGc                                                                                                          >  2:766946/1‑148 (MQ=255)
                                              ttttCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAgg                                                                                                       <  1:3306407/148‑1 (MQ=255)
                                                 tCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCg                                                                                                    <  2:481078/148‑1 (MQ=255)
                                                     gcTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGa                                                                                                <  2:905771/148‑1 (MQ=255)
                                                     gcTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGa                                                                                                <  2:2353585/148‑1 (MQ=255)
                                                      cTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGaa                                                                                               <  1:985103/148‑1 (MQ=255)
                                                       tGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAAt                                                                                              <  2:954813/148‑1 (MQ=255)
                                                         gTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCt                                                                                            >  1:1467196/1‑148 (MQ=255)
                                                          ttGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTg                                                                                           >  1:1345260/1‑148 (MQ=255)
                                                           tGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGt                                                                                          >  1:2961829/1‑148 (MQ=255)
                                                           tGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGt                                                                                          <  2:1430800/148‑1 (MQ=255)
                                                           tGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGt                                                                                          <  2:257444/148‑1 (MQ=255)
                                                            gAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGtt                                                                                         <  2:3287560/148‑1 (MQ=255)
                                                             aaTCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTa                                                                                        >  1:66051/1‑148 (MQ=255)
                                                             aaTCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTa                                                                                        <  2:128783/148‑1 (MQ=255)
                                                              aTCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTAc                                                                                       <  2:1253776/148‑1 (MQ=255)
                                                                cAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACgg                                                                                     >  1:2199621/1‑148 (MQ=255)
                                                                   cccGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGTCAGAGCAGGTCGATGAATCTGTTACggcag                                                                                  >  2:3289905/1‑148 (MQ=255)
                                                                     cGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACggcaggc                                                                                <  2:472573/148‑1 (MQ=255)
                                                                      gACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACggcaggca                                                                               <  1:778888/148‑1 (MQ=255)
                                                                      gACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACggcaggca                                                                               <  2:1070017/148‑1 (MQ=255)
                                                                       aCTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCAt                                                                              <  1:1100799/148‑1 (MQ=255)
                                                                         tGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGa                                                                            <  1:1630392/148‑1 (MQ=255)
                                                                          ggggACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGaa                                                                           >  2:2467056/1‑148 (MQ=255)
                                                                           gggACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAAt                                                                          >  2:1790522/1‑148 (MQ=255)
                                                                            ggACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATa                                                                         <  1:1697814/148‑1 (MQ=255)
                                                                            ggACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATa                                                                         >  2:199296/1‑148 (MQ=255)
                                                                            ggACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATa                                                                         >  2:774449/1‑148 (MQ=255)
                                                                              aCGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACt                                                                       <  1:44042/148‑1 (MQ=255)
                                                                                ggggAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACttt                                                                     <  1:1499130/148‑1 (MQ=255)
                                                                                 gggAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACtttt                                                                    <  1:3060870/148‑1 (MQ=255)
                                                                                 gggAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACtttt                                                                    <  1:551628/148‑1 (MQ=255)
                                                                                        gagaAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTg                                                             <  2:857482/148‑1 (MQ=255)
                                                                                           aaaTCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGt                                                          <  2:380645/148‑1 (MQ=255)
                                                                                              tCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCg                                                       <  2:2396890/148‑1 (MQ=255)
                                                                                              tCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCg                                                       <  1:3233632/148‑1 (MQ=255)
                                                                                                 ccGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCtt                                                    >  1:1892488/1‑148 (MQ=255)
                                                                                                 ccGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCtt                                                    <  2:285314/148‑1 (MQ=255)
                                                                                                 ccGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCtt                                                    <  2:288258/148‑1 (MQ=255)
                                                                                                     cAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATg                                                <  2:2340295/148‑1 (MQ=255)
                                                                                                      aGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGt                                               >  2:1380352/1‑148 (MQ=255)
                                                                                                      aGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGt                                               <  2:735349/148‑1 (MQ=255)
                                                                                                       gAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGtt                                              >  1:1720572/1‑148 (MQ=255)
                                                                                                           cGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGAc                                          >  2:1209790/1‑148 (MQ=255)
                                                                                                              ttGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGa                                       >  2:1852559/1‑148 (MQ=255)
                                                                                                              ttGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGa                                       <  2:2308866/148‑1 (MQ=255)
                                                                                                              ttGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGa                                       >  2:1872969/1‑148 (MQ=255)
                                                                                                                 aTTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAAtt                                    >  2:806835/1‑148 (MQ=255)
                                                                                                                 aTTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAAtt                                    >  2:1761811/1‑148 (MQ=255)
                                                                                                                   tACAGTGTGTTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTAc                                  >  1:1849782/1‑148 (MQ=255)
                                                                                                                    aCAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACa                                 >  1:3159323/1‑148 (MQ=255)
                                                                                                                        tgtgtTGCAACGAGCATTAAGCGCAGACGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGAct                             <  1:3289905/148‑2 (MQ=255)
                                                                                                                        tgtgtTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGAcg                             <  1:638584/148‑1 (MQ=255)
                                                                                                                           gtTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGt                          <  2:2024687/148‑1 (MQ=255)
                                                                                                                           gtTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGt                          <  2:1153971/148‑1 (MQ=255)
                                                                                                                                aaCGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTAcg                     <  2:2186355/148‑1 (MQ=255)
                                                                                                                                 aCGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTAcgc                    >  1:2361284/1‑148 (MQ=255)
                                                                                                                                       aTTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGt              >  1:2144723/1‑148 (MQ=255)
                                                                                                                                        ttAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTc             >  2:3095859/1‑148 (MQ=255)
                                                                                                                                         tAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTCt            >  1:770592/1‑148 (MQ=255)
                                                                                                                                         tAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTCt            >  1:2338654/1‑148 (MQ=255)
                                                                                                                                          aaGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTCTg           >  1:2846084/1‑148 (MQ=255)
                                                                                                                                           aGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCTCTGAGTTCGCTTCATGTTTGACTGAATTACAGACGCGTTTACGCGAAAGTCTgg          >  1:1153176/1‑148 (MQ=255)
                                                                                                                                                   aGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTCTGGCTCTGCGt  <  1:2467056/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCAACGGACGAGCAGAGCACGCGTTTAGGCGTTGTGCGTGAGACTTTTTCGCGCTGGTTGAATCAGCCCGACTGGGGACGGGGAAATGAGAAATCACCGACAGAAACGGTTGATTACAGTGTGTTGCAACGAGCATTAAGCGCAGAAGTCGATCTCTATTGCCAACAAATGGCTAAAGTTCTGGCAGAGCAGGTCGATGAATCTGTTACGGCAGGCATGAATACTTTTTTCGCTGAGTTCGCTTCATGTTTGACGGAATTACAGACGCGTTTACGCGAAAGTCTGGCTCTGCGT  >  CP009273/4318751‑4319045

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: