Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 412,258 C→A 6.4% D51Y (GAT→TAT)  sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273412,2580CA 6.4% 807.4 / 20.3 217D51Y (GAT→TAT) sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (97/106);  new base A (5/9);  total (102/115)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 217 total aligned reads displayed.

ACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAGGATG  >  CP009273/412114‑412403
                                                                                                                                                |                                                                                                                                                 
     aTGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACAt                                                                                                                                           >  2:1867483/1‑148 (MQ=255)
     aTGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATAGAAAACAt                                                                                                                                           >  2:3309579/1‑148 (MQ=255)
      tGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATc                                                                                                                                          <  1:1801677/148‑1 (MQ=255)
      tGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATAGAAAACATc                                                                                                                                          >  1:2076505/1‑148 (MQ=255)
       gATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCa                                                                                                                                         <  2:170086/148‑1 (MQ=255)
       gATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCa                                                                                                                                         <  1:1267512/148‑1 (MQ=255)
       gATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCa                                                                                                                                         <  2:2818395/148‑1 (MQ=255)
        atatCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCAc                                                                                                                                        >  2:2672664/1‑148 (MQ=255)
        atatCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCAc                                                                                                                                        <  1:3306881/148‑1 (MQ=255)
          atCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCg                                                                                                                                      >  2:385834/1‑148 (MQ=255)
              cgcgATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGcaa                                                                                                                                  <  2:141493/148‑1 (MQ=255)
                 gcctcattcaGCGTGGCGACCGAGTCATGGTTTCCTGCCAGGACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGcaacaa                                                                                                                               <  1:1994280/145‑1 (MQ=255)
                  attcattcAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaat                                                                                                                              >  1:2456351/1‑148 (MQ=255)
                  attcattcAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACaat                                                                                                                              >  1:1143384/1‑148 (MQ=255)
                       ttcaGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCg                                                                                                                         <  2:1158482/148‑1 (MQ=255)
                        tcaGCGTGGCGACCGAGTCATGGTTTCCTGACAGTACCACCAGATGACAGACAGTTTGCTTTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGc                                                                                                                        <  1:43303/148‑1 (MQ=255)
                         caGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCa                                                                                                                       >  2:2134982/1‑148 (MQ=255)
                          aGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCTGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCAt                                                                                                                      <  2:3171592/148‑1 (MQ=255)
                              tGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCAc                                                                                                                  <  1:824092/148‑1 (MQ=255)
                                gCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCt                                                                                                                >  2:3447541/1‑148 (MQ=255)
                                 cGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTg                                                                                                               <  2:3403674/148‑1 (MQ=255)
                                    ccGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATAGAAAACATCACCGGCAACAATAATCGCATCCACCTGATg                                                                                                            >  1:184531/1‑148 (MQ=255)
                                       aGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGt                                                                                                         >  2:1100029/1‑148 (MQ=255)
                                       aGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGt                                                                                                         >  2:209475/1‑148 (MQ=255)
                                       aGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGt                                                                                                         >  2:679206/1‑148 (MQ=255)
                                         tCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGttt                                                                                                       >  1:2479080/1‑148 (MQ=255)
                                          cATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTtg                                                                                                      >  2:1719821/1‑148 (MQ=255)
                                          cATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTtg                                                                                                      <  1:123658/148‑1 (MQ=255)
                                          cATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTtg                                                                                                      >  2:2612175/1‑148 (MQ=255)
                                          cATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTtg                                                                                                      <  1:3501770/148‑1 (MQ=255)
                                             ggTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGc                                                                                                   >  1:66482/1‑148 (MQ=255)
                                               tttCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTg                                                                                                 <  2:3545312/148‑1 (MQ=255)
                                                    tGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTcc                                                                                            >  1:75244/1‑148 (MQ=255)
                                                     gCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCca                                                                                           >  1:1553696/1‑148 (MQ=255)
                                                         gTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGGGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCcagcag                                                                                       >  1:2101798/1‑148 (MQ=255)
                                                          tACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCagcagc                                                                                      <  2:135404/148‑1 (MQ=255)
                                                             caccaGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAg                                                                                   >  1:1373386/1‑148 (MQ=255)
                                                              accaGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGt                                                                                  >  2:3026570/1‑148 (MQ=255)
                                                              accaGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGt                                                                                  <  2:2716128/148‑1 (MQ=255)
                                                                  gATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAg                                                                              <  1:727755/148‑1 (MQ=255)
                                                                  gATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAg                                                                              <  1:244293/148‑1 (MQ=255)
                                                                   aTGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATGCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGa                                                                             >  2:1579777/1‑148 (MQ=255)
                                                                   aTGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGa                                                                             >  2:2169735/1‑148 (MQ=255)
                                                                      acccccAGTTTGCTGCAACTTGACCACACAACGGTTGTATCACGAGCGGGCGTAACTGGGCGGCGAGCCGGTAACGCACACATCACCGGCAACACTAAACGCATCCACCCGAACGGCTTGTGCTGTCTCCAGCAGCCAGTCAAGaaaa                                                                          <  2:1462255/144‑1 (MQ=255)
                                                                      aCAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGaaaa                                                                          >  2:1650618/1‑148 (MQ=255)
                                                                        aGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGc                                                                        <  1:3133694/148‑1 (MQ=255)
                                                                         gCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGcc                                                                       <  1:1121589/148‑1 (MQ=255)
                                                                         gCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGcc                                                                       >  1:3153019/1‑148 (MQ=255)
                                                                         gCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGcc                                                                       <  1:1243351/148‑1 (MQ=255)
                                                                         gCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGcc                                                                       <  2:3185303/148‑1 (MQ=255)
                                                                         gCCAGTGTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGcc                                                                       <  2:3276737/148‑1 (MQ=255)
                                                                          ccAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCt                                                                      <  2:892146/148‑1 (MQ=255)
                                                                          ccAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCt                                                                      <  2:2743407/148‑1 (MQ=255)
                                                                           cAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTg                                                                     <  1:3158001/148‑1 (MQ=255)
                                                                           cAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATAGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTg                                                                     <  2:3104098/148‑1 (MQ=255)
                                                                             gTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGAt                                                                   >  2:1907715/1‑148 (MQ=255)
                                                                             gTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGAt                                                                   <  1:454442/148‑1 (MQ=255)
                                                                              tttGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGTCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATg                                                                  <  2:2903945/148‑1 (MQ=255)
                                                                              tttGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACACTAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCGGATg                                                                  >  2:3146253/1‑148 (MQ=255)
                                                                                   tGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAg                                                             <  2:1599867/148‑1 (MQ=255)
                                                                                     tAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCt                                                           <  2:2479080/148‑1 (MQ=255)
                                                                                       aaTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTc                                                         >  1:9322/1‑148 (MQ=255)
                                                                                       aaTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTc                                                         >  2:1946740/1‑148 (MQ=255)
                                                                                        aTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTcg                                                        <  1:1000513/148‑1 (MQ=255)
                                                                                             caacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCt                                                   >  2:712508/1‑148 (MQ=255)
                                                                                             caacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCt                                                   <  1:2167718/148‑1 (MQ=255)
                                                                                             caacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCt                                                   <  1:1696048/148‑1 (MQ=255)
                                                                                              aacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCtt                                                  >  2:1793799/1‑148 (MQ=255)
                                                                                               acaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCttt                                                 >  2:3289453/1‑148 (MQ=255)
                                                                                                  aaaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTAc                                              >  2:2565226/1‑148 (MQ=255)
                                                                                                   aaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt                                             >  2:1756566/1‑148 (MQ=255)
                                                                                                   aaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt                                             >  2:781097/1‑148 (MQ=255)
                                                                                                   aaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt                                             >  2:781052/1‑148 (MQ=255)
                                                                                                    aCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTg                                            >  1:2577419/1‑148 (MQ=255)
                                                                                                    aCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTg                                            >  1:3474955/1‑148 (MQ=255)
                                                                                                     cGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGt                                           >  1:2160503/1‑148 (MQ=255)
                                                                                                      ggTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTa                                          >  1:1189114/1‑148 (MQ=255)
                                                                                                        ttGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGa                                        <  1:649401/148‑1 (MQ=255)
                                                                                                           tataACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGt                                     >  1:267487/1‑148 (MQ=255)
                                                                                                            ataACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGtt                                    <  1:1776199/148‑1 (MQ=255)
                                                                                                            ataACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGtt                                    <  2:2031462/148‑1 (MQ=255)
                                                                                                             taACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTc                                   <  2:2302681/148‑1 (MQ=255)
                                                                                                                    cGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCga                            >  1:3399332/1‑148 (MQ=255)
                                                                                                                     gggCGTAACTGGGCGGCGAGCCGGTATCGAAAACCTCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCgag                           <  1:1488162/148‑1 (MQ=255)
                                                                                                                      ggCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCgaga                          >  2:475183/1‑148 (MQ=255)
                                                                                                                         gTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGc                       <  1:2975016/148‑1 (MQ=255)
                                                                                                                         gTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGc                       <  1:3414729/148‑1 (MQ=255)
                                                                                                                          tAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGTTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGcc                      <  1:2645213/148‑1 (MQ=255)
                                                                                                                            aCTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAg                    >  1:1020540/1‑148 (MQ=255)
                                                                                                                             cTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGt                   >  1:2228825/1‑148 (MQ=255)
                                                                                                                                      gAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGgtgt          <  1:3026570/148‑1 (MQ=255)
                                                                                                                                       aGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGtgtg         <  1:874868/148‑1 (MQ=255)
                                                                                                                                         ccGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGaa       >  1:2702949/1‑148 (MQ=255)
                                                                                                                                         ccGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGaa       >  1:3353686/1‑148 (MQ=255)
                                                                                                                                           ggTATCGACAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAgg     <  1:888254/148‑1 (MQ=255)
                                                                                                                                           ggTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAgg     <  1:696384/148‑1 (MQ=255)
                                                                                                                                             tATCGAAAACATCACCGGCAACAATAATCGCATCCCCCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAGGAt   <  2:1476563/148‑1 (MQ=255)
                                                                                                                                              aTCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAGGATg  <  1:2766303/148‑1 (MQ=255)
                                                                                                                                              aTCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAGGATg  <  2:1789976/148‑1 (MQ=255)
                                                                                                                                              aTCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAGGATg  >  2:2366966/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGAGTCATGGTTTCCTGCCAGTACCACCAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTGGCCGAGATGCCAGTCTGAGGTGTGAAGGATG  >  CP009273/412114‑412403

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: