Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 412,004 A→T 5.5% R135R (CGT→CGA sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273412,0040AT 5.5% 965.2 / 31.9 316R135R (CGT→CGAsbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (170/127);  new base T (8/10);  total (178/138)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

AAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGA  >  CP009273/411855‑412153
                                                                                                                                                     |                                                                                                                                                     
aaaTGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCa                                                                                                                                                       <  2:2867187/150‑1 (MQ=255)
 aaTGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATATTCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCAc                                                                                                                                                      <  1:1527871/150‑1 (MQ=255)
 aaTGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCAc                                                                                                                                                      <  1:4032779/150‑1 (MQ=255)
 aaTGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCAAGTAAATGCTGCTGTTTTTCAATACCTTTAAGCCCCGCCTGGCTGGTAATAATGTCAc                                                                                                                                                      >  2:459401/1‑150 (MQ=255)
  aTGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCAcg                                                                                                                                                     <  2:4094432/150‑1 (MQ=255)
    gTCCCGTGGCGGTGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCTCCTGGCTGGTAATAATGTCAcgcg                                                                                                                                                   >  2:4786982/1‑150 (MQ=255)
     tCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCgg                                                                                                                                                  <  2:1488612/150‑1 (MQ=255)
          tGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTa                                                                                                                                             <  2:3841260/150‑1 (MQ=255)
          tGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTa                                                                                                                                             <  1:1987007/150‑1 (MQ=255)
            gCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTaaa                                                                                                                                           <  1:1108444/150‑1 (MQ=255)
              gatgatGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTaaaaa                                                                                                                                         <  1:1004957/150‑1 (MQ=255)
              gatgatGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTaaaaa                                                                                                                                         >  1:2100606/1‑150 (MQ=255)
              gatgatGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTaaaaa                                                                                                                                         >  1:4332642/1‑150 (MQ=255)
               atgatgGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAAc                                                                                                                                        <  2:4965433/150‑1 (MQ=255)
               atgatgGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAAc                                                                                                                                        <  1:3899486/150‑1 (MQ=255)
                tgatgGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACg                                                                                                                                       <  1:4354768/150‑1 (MQ=255)
                   tgGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGaa                                                                                                                                    <  1:658337/150‑1 (MQ=255)
                   tgGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGaa                                                                                                                                    <  1:2797710/150‑1 (MQ=255)
                   tgGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATAGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATTCCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGaa                                                                                                                                    >  1:2655605/1‑150 (MQ=255)
                    gggCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACATAAAAACGGAAt                                                                                                                                   >  1:2039127/1‑150 (MQ=255)
                      gCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCAGAGTGTGGGTGGTCATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATgg                                                                                                                                 >  2:2496066/1‑150 (MQ=255)
                        agagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATgggg                                                                                                                               >  1:1466358/1‑150 (MQ=255)
                        agagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATgggg                                                                                                                               >  2:5084029/1‑150 (MQ=255)
                        agagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATgggg                                                                                                                               >  1:516117/1‑150 (MQ=255)
                        agagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATgggg                                                                                                                               >  1:3595820/1‑150 (MQ=255)
                        agagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATgggg                                                                                                                               <  2:3646711/150‑1 (MQ=255)
                         gagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGc                                                                                                                              <  2:612359/150‑1 (MQ=255)
                         gagGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGc                                                                                                                              >  2:558151/1‑150 (MQ=255)
                          agGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGca                                                                                                                             <  1:3204894/150‑1 (MQ=255)
                           ggCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCTCGCGGACGTAAAAACGGAATGGGGcac                                                                                                                            <  2:2373586/150‑1 (MQ=255)
                             cTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAg                                                                                                                          <  1:3045566/150‑1 (MQ=255)
                             cTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAg                                                                                                                          >  2:4303911/1‑150 (MQ=255)
                              tGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGc                                                                                                                         <  1:3414777/150‑1 (MQ=255)
                                  cgccgcGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTg                                                                                                                     >  2:1840212/1‑150 (MQ=255)
                                    ccgcGCAGTTTGCAGGCATCGGCATAGTTTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTgcg                                                                                                                   <  2:3404/150‑1 (MQ=255)
                                    ccgcGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTgcg                                                                                                                   <  1:1576509/150‑1 (MQ=255)
                                      gcgcAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGcc                                                                                                                 <  1:689027/150‑1 (MQ=255)
                                                cAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCg                                                                                                       >  2:3196590/1‑150 (MQ=255)
                                                 aGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGt                                                                                                      >  2:3338785/1‑150 (MQ=255)
                                                  ggCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGACTGGGGCACAGCACTGCGCCTGGCGTCCCGTc                                                                                                     >  1:4358748/1‑150 (MQ=255)
                                                  ggCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTc                                                                                                     >  1:3653684/1‑150 (MQ=255)
                                                   gCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTCTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTcg                                                                                                    >  2:1114798/1‑150 (MQ=255)
                                                             aGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCa                                                                                          >  1:4677974/1‑150 (MQ=255)
                                                               tgttgtTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAg                                                                                        >  2:3795259/1‑150 (MQ=255)
                                                                 ttgttgGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGAt                                                                                      <  1:891701/150‑1 (MQ=255)
                                                                   gttgGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGAttt                                                                                    <  1:4104639/150‑1 (MQ=255)
                                                                   gttgGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGAttt                                                                                    >  1:5224051/1‑150 (MQ=255)
                                                                      ggTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTgcg                                                                                 >  2:858278/1‑150 (MQ=255)
                                                                       gTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTgcgg                                                                                >  2:2588190/1‑150 (MQ=255)
                                                                       gTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTgcgg                                                                                >  1:3847062/1‑150 (MQ=255)
                                                                       gTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTgcgg                                                                                >  1:610034/1‑150 (MQ=255)
                                                                        taataaTCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTgcggc                                                                               >  1:1397686/1‑150 (MQ=255)
                                                                        taataaTCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAAGACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTgcggc                                                                               <  2:4875287/150‑1 (MQ=255)
                                                                           aaatCGGTAAACTCTTCCAGAAAAAGCTGCTGTTTTTCAATACCGTTAATCCCCTCCGGGCTGGTAAAAATGACACGCGGACGAAAAAACGGAATGGGTCACAGCACTTCGACTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCa                                                                            >  2:1342706/2‑150 (MQ=255)
                                                                            aatCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCAt                                                                           <  2:2901075/150‑1 (MQ=255)
                                                                             atCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATg                                                                          >  1:4291735/1‑150 (MQ=255)
                                                                                ggTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTcc                                                                       <  2:2305503/150‑1 (MQ=255)
                                                                                      cGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCt                                                                 >  1:1072257/1‑150 (MQ=255)
                                                                                       gCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTg                                                                >  1:2400628/1‑150 (MQ=255)
                                                                                       gCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTg                                                                >  2:4279013/1‑150 (MQ=255)
                                                                                        cTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTgg                                                               >  1:5127319/1‑150 (MQ=255)
                                                                                         tGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGc                                                              >  1:4190097/1‑150 (MQ=255)
                                                                                            cAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGAc                                                           >  2:5199505/1‑150 (MQ=255)
                                                                                            cAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGAc                                                           <  1:5075042/150‑1 (MQ=255)
                                                                                             aGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGAcc                                                          >  1:758829/1‑150 (MQ=255)
                                                                                             aGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCAGGTCCGGCGCTGGCGAcc                                                          >  2:1734009/1‑150 (MQ=255)
                                                                                              gTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCa                                                         <  1:1797382/150‑1 (MQ=255)
                                                                                                  aTGCTGCTGTTTTTCAATACCGTTCAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGgt                                                     <  2:4945105/150‑1 (MQ=255)
                                                                                                  aTGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCTCGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGgt                                                     >  1:4480606/1‑150 (MQ=255)
                                                                                                  aTGCTGCTGTTTTTCAAGACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGgt                                                     <  2:5275975/150‑1 (MQ=255)
                                                                                                     ctgctgTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGgtagt                                                  <  2:157337/150‑1 (MQ=255)
                                                                                                     ctgctgTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGgtagt                                                  >  1:2749923/1‑150 (MQ=255)
                                                                                                        ctgTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCTCGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTAtt                                               <  1:1204917/150‑1 (MQ=255)
                                                                                                        ctgTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTAtt                                               >  1:2219447/1‑150 (MQ=255)
                                                                                                                cAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACg                                       >  1:4414570/1‑150 (MQ=255)
                                                                                                                  aTACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGcc                                     >  1:2554978/1‑150 (MQ=255)
                                                                                                                     ccGTTAAGCCCCGCCTGGCGGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCCAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGACCGCCATg                                  >  2:3295039/1‑150 (MQ=255)
                                                                                                                      cGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGa                                 <  1:1270006/150‑1 (MQ=255)
                                                                                                                       gTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGat                                >  2:835148/1‑150 (MQ=255)
                                                                                                                       gTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCATCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGat                                >  2:693472/1‑150 (MQ=255)
                                                                                                                        ttAAGCCCCGCCTGGCTGGTAATAATGTCTCGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCTCGGTAGTATTGAGGAACGCCATGata                               <  1:3482243/150‑1 (MQ=255)
                                                                                                                         tAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGatat                              <  2:2704530/150‑1 (MQ=255)
                                                                                                                         tAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGatat                              <  2:5112441/150‑1 (MQ=255)
                                                                                                                          aaGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATc                             <  1:1610458/150‑1 (MQ=255)
                                                                                                                          aaGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATc                             >  2:2418490/1‑150 (MQ=255)
                                                                                                                           aGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTTGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATcg                            <  1:4843171/150‑1 (MQ=255)
                                                                                                                                  ccTGGCTGGTAATAATGTCTCGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGatt                     >  1:2520707/1‑150 (MQ=255)
                                                                                                                                  ccTGGCTGGTAATAATGTCTCGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGatt                     >  1:5165079/1‑150 (MQ=255)
                                                                                                                                    tggccggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGattca                   <  2:3229048/150‑1 (MQ=255)
                                                                                                                                     ggctggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGattcat                  <  1:4955515/150‑1 (MQ=255)
                                                                                                                                     ggctggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGattcat                  <  2:2100606/150‑1 (MQ=255)
                                                                                                                                     ggctggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGattcat                  <  1:558151/150‑1 (MQ=255)
                                                                                                                                      gctggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGattcatt                 >  2:5235633/1‑150 (MQ=255)
                                                                                                                                       ctggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGattcattc                <  1:1754930/150‑1 (MQ=255)
                                                                                                                                         ggTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAg              >  1:995747/1‑150 (MQ=255)
                                                                                                                                          gTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGc             <  1:4547551/150‑1 (MQ=255)
                                                                                                                                            aataatGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGGAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGt           >  2:3391482/1‑150 (MQ=255)
                                                                                                                                             ataatGTCACGCGGACGTAAAAACGGAATGGGGAACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTg          >  1:1914161/1‑150 (MQ=255)
                                                                                                                                              taatGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTgg         >  2:3067667/1‑150 (MQ=255)
                                                                                                                                               aatGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGc        >  1:3125004/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAATCGCTGCCAGTAAATGCTGCTGTTTTTCAATACCGTTAAGCCCCGCCTGGCTGGTAATAATGTCACGCGGACGTAAAAACGGAATGGGGCACAGCACTGCGCCTGGCGTCCCGTCGCGACGAGGCAAGATTTGCGGCGCATGTCCGGCGCTGGCGACCACGGTAGTATTGAGGAACGCCATGATATCGCGCGATTCATTCAGCGTGGCGACCGA  >  CP009273/411855‑412153

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: