Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,945,504 C→A 16.5% R106L (CGT→CTT)  recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,945,5040CA16.5% 790.8 / 105.3 259R106L (CGT→CTT) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (107/109);  new base A (21/22);  total (128/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

CTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGT  >  CP009273/2945355‑2945653
                                                                                                                                                     |                                                                                                                                                     
cTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAc                                                                                                                                                       >  1:1411078/1‑150 (MQ=255)
  ttttGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGa                                                                                                                                                     <  1:3129560/150‑1 (MQ=255)
     tGCCAGTTAATTTCATCGCTTACTGGAGAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCg                                                                                                                                                  <  1:2538994/150‑1 (MQ=255)
     tGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCg                                                                                                                                                  <  1:1110436/150‑1 (MQ=255)
       ccAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGcc                                                                                                                                                <  1:3445345/150‑1 (MQ=255)
         aGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCcac                                                                                                                                              >  1:2952998/1‑150 (MQ=255)
             aaTTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGAGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTACGTTGAAAAAGCGAGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGACACagag                                                                                                                                          >  1:782168/1‑150 (MQ=255)
              aTTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAgg                                                                                                                                         >  2:1004940/1‑150 (MQ=255)
               tttCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGa                                                                                                                                        <  1:1198843/150‑1 (MQ=255)
                 tCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGatc                                                                                                                                      <  2:844610/150‑1 (MQ=255)
                 tCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGatc                                                                                                                                      <  1:280415/150‑1 (MQ=255)
                  cATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGatca                                                                                                                                     <  1:4227422/150‑1 (MQ=255)
                   aTCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGatcat                                                                                                                                    <  1:2009049/150‑1 (MQ=255)
                    tCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGatcatc                                                                                                                                   >  1:1858659/1‑150 (MQ=255)
                     cGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCg                                                                                                                                  <  1:4283874/150‑1 (MQ=255)
                     cGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCg                                                                                                                                  >  1:1457327/1‑150 (MQ=255)
                      gCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCgg                                                                                                                                 <  1:259309/150‑1 (MQ=255)
                       cTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGc                                                                                                                                >  2:2367857/1‑150 (MQ=255)
                        ttACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCg                                                                                                                               >  2:2531244/1‑150 (MQ=255)
                              gAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGtt                                                                                                                         >  1:2055926/1‑150 (MQ=255)
                              gAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGtt                                                                                                                         >  2:3694279/1‑150 (MQ=255)
                               aaaaaGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTc                                                                                                                        <  2:2513696/150‑1 (MQ=255)
                                  aaGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCAtc                                                                                                                     <  1:4061326/150‑1 (MQ=255)
                                   aGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCAtct                                                                                                                    >  2:2632736/1‑150 (MQ=255)
                                   aGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCAtct                                                                                                                    >  2:5582/1‑150 (MQ=255)
                                    gTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCAtctc                                                                                                                   >  1:1665332/1‑150 (MQ=255)
                                    gTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCAtctc                                                                                                                   >  1:4535698/1‑150 (MQ=255)
                                     tttGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTcc                                                                                                                  >  1:1734488/1‑150 (MQ=255)
                                           cAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCctgct                                                                                                            <  1:17859/150‑1 (MQ=255)
                                            aGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCctgctg                                                                                                           >  2:1444591/1‑150 (MQ=255)
                                              ggTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGAc                                                                                                         <  2:347552/150‑1 (MQ=255)
                                               gTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGAcc                                                                                                        >  2:990739/1‑150 (MQ=255)
                                                tttGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCg                                                                                                       <  2:3888515/150‑1 (MQ=255)
                                                  tGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCt                                                                                                     <  1:3667397/150‑1 (MQ=255)
                                                  tGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCt                                                                                                     <  1:3728368/150‑1 (MQ=255)
                                                  tGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCt                                                                                                     >  2:3771766/1‑150 (MQ=255)
                                                      ccAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTgag                                                                                                 >  1:3440015/1‑150 (MQ=255)
                                                       cAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTgaga                                                                                                >  2:71379/1‑150 (MQ=255)
                                                            gagCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCa                                                                                           >  1:3818886/1‑150 (MQ=255)
                                                            gagCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCa                                                                                           <  1:3221283/150‑1 (MQ=255)
                                                             agCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCag                                                                                          >  2:2911925/1‑150 (MQ=255)
                                                                 tcatcaACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCagcaag                                                                                      >  2:2837422/1‑150 (MQ=255)
                                                                 tcatcaACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCagcaag                                                                                      >  2:2236002/1‑150 (MQ=255)
                                                                  catcaACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCagcaagc                                                                                     <  1:3467603/150‑1 (MQ=255)
                                                                   atcaACCTCAATGTCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCagcaagca                                                                                    <  1:2096599/150‑1 (MQ=255)
                                                                   atcaACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCagcaagca                                                                                    <  2:4273762/150‑1 (MQ=255)
                                                                      aaCCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCAttc                                                                                 >  1:2951320/1‑150 (MQ=255)
                                                                        ccTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCAttctt                                                                               >  2:4512167/1‑150 (MQ=255)
                                                                         cTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCAttcttc                                                                              <  2:3536620/150‑1 (MQ=255)
                                                                          tCAATGGCATGATTCACTTCGTCGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTcc                                                                             >  2:1669473/1‑150 (MQ=255)
                                                                           cAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTccc                                                                            >  1:2525598/1‑150 (MQ=255)
                                                                            aaTGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCa                                                                           <  2:2405216/150‑1 (MQ=255)
                                                                               ggCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAATCCAGCAAGCATTCTTCCCAAtt                                                                        >  2:1471226/1‑150 (MQ=255)
                                                                                        cACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTc                                                               >  2:416283/1‑150 (MQ=255)
                                                                                        cACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTc                                                               <  2:1679410/150‑1 (MQ=255)
                                                                                           ttCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCAcc                                                            <  2:2037370/150‑1 (MQ=255)
                                                                                             cGTTGAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGAt                                                         >  2:1093449/1‑150 (MQ=255)
                                                                                                tGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGAttt                                                       >  2:1968176/1‑150 (MQ=255)
                                                                                                   aaaaGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCAc                                                    >  2:4120/1‑150 (MQ=255)
                                                                                                   aaaaGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCAc                                                    <  2:2659/150‑1 (MQ=255)
                                                                                                    aaaGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGAttt                                                       <  1:451130/146‑1 (MQ=255)
                                                                                                    aaaGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGAttt                                                       <  1:451292/146‑1 (MQ=255)
                                                                                                           gCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATAATCGGCTAGGGTTCATCACCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTACCCAATTATGTAGCTCACCGATTTCACAGACACAg                                            >  2:2540647/1‑150 (MQ=255)
                                                                                                             cACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGt                                          <  1:3694279/150‑1 (MQ=255)
                                                                                                              aCTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTc                                         <  2:391587/150‑1 (MQ=255)
                                                                                                                 gtgCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGcc                                      >  1:868635/1‑150 (MQ=255)
                                                                                                                    cgcTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAAc                                   <  2:4081157/150‑1 (MQ=255)
                                                                                                                        cGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCg                               >  1:949606/1‑150 (MQ=255)
                                                                                                                         gTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACAGATTTAACTGACACAGGTCGCCAACAGCgg                              >  1:2649883/1‑150 (MQ=255)
                                                                                                                          ttACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGa                             >  1:2330529/1‑150 (MQ=255)
                                                                                                                          ttACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGa                             >  1:64059/1‑150 (MQ=255)
                                                                                                                          ttACACCACATGCAATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGa                             >  1:303182/1‑150 (MQ=255)
                                                                                                                          ttACACCACAGGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGa                             <  2:4021612/150‑1 (MQ=255)
                                                                                                                           tACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGAt                            <  2:1809623/150‑1 (MQ=255)
                                                                                                                           tACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGAt                            <  2:3150547/150‑1 (MQ=255)
                                                                                                                               ccacATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGa                        <  2:1730296/150‑1 (MQ=255)
                                                                                                                                cacaTGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGAc                       <  2:1408598/150‑1 (MQ=255)
                                                                                                                                cacaTGCGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGAc                       <  1:2841172/150‑1 (MQ=255)
                                                                                                                                  caTGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGc                     >  2:2622801/1‑150 (MQ=255)
                                                                                                                                   aTGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGcc                    <  1:2701987/150‑1 (MQ=255)
                                                                                                                                      cGATTCAAATAACGCCGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCg                 <  1:2434530/150‑1 (MQ=255)
                                                                                                                                      cGATTCAAATAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCg                 >  2:2571862/1‑150 (MQ=255)
                                                                                                                                        aTTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGtt               <  2:2709115/150‑1 (MQ=255)
                                                                                                                                          tCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTAt             <  1:1004940/150‑1 (MQ=255)
                                                                                                                                          tCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTAt             <  1:5582/150‑1 (MQ=255)
                                                                                                                                            gaaTAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTTACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGCTGCGACGCCTCGTTAttt           >  1:4240640/2‑150 (MQ=255)
                                                                                                                                            aaaTAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTAttt           >  1:4112607/1‑150 (MQ=255)
                                                                                                                                            aaaTAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATATTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTAttt           >  1:1304953/1‑150 (MQ=255)
                                                                                                                                            aaaTAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTAttt           <  1:3158952/150‑1 (MQ=255)
                                                                                                                                            aaaTAAAGAAGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTAttt           <  1:892198/150‑1 (MQ=255)
                                                                                                                                              aTAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTc         <  1:2911925/150‑1 (MQ=255)
                                                                                                                                              aTAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTc         >  1:3672834/1‑150 (MQ=255)
                                                                                                                                               tAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTcc        >  2:4279583/1‑150 (MQ=255)
                                                                                                                                               tAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTcc        >  2:4360224/1‑150 (MQ=255)
                                                                                                                                                aaaGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCa       >  2:3667495/1‑150 (MQ=255)
                                                                                                                                                aaaGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCa       <  2:4535698/150‑1 (MQ=255)
                                                                                                                                                aaaGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCa       <  2:1665301/150‑1 (MQ=255)
                                                                                                                                                 aaGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAg      <  1:3916325/150‑1 (MQ=255)
                                                                                                                                                  agaagaTCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGt     <  1:195267/150‑1 (MQ=255)
                                                                                                                                                   gACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTc    >  1:1085271/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTGAAAAAGCGTGCCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGT  >  CP009273/2945355‑2945653

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: