Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,759,184 C→T 18.3% intergenic (‑49/‑93) yfjM ← / → rnlA CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,759,1840CT18.3% 841.7 / 166.6 333intergenic (‑49/‑93)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base C (125/147);  new base T (32/29);  total (157/176)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 339 total aligned reads displayed.

AGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTGA  >  CP009273/2759036‑2759332
                                                                                                                                                    |                                                                                                                                                    
  aaTTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTGGACGATCCTGCCCTATAGGTGGAGTGcccc                                                                                                                                                   >  1:233780/1‑150 (MQ=255)
  aaTTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTCTAGTTGGAGTGtccc                                                                                                                                                   >  2:5651305/1‑150 (MQ=255)
    ttttttCTTGTCTGCTGGCTAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                                                                                                                 <  2:2260259/150‑1 (MQ=255)
     tttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTa                                                                                                                                                <  1:967011/150‑1 (MQ=255)
       tttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTAtt                                                                                                                                              <  1:1069770/150‑1 (MQ=255)
        ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCGGCCCTGGAGTTGGGGTGCCCCCtcgtc                                                                                                                                             >  1:1756139/1‑146 (MQ=255)
          cTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAg                                                                                                                                           >  2:3518058/1‑150 (MQ=255)
             gTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTg                                                                                                                                        <  1:2171738/150‑1 (MQ=255)
             gTCAGCAGGCGATTAATGACTCATTAACGAGTGGCTGCTGTGTGAGCTCAGGCAGAGGTGACCGAACCGCTGGACATCATCCCCCATTTTGTGGTCTGCCCCCCCGGTAGCCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTg                                                                                                                                        <  2:2645851/150‑1 (MQ=255)
              tCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                                                                                                                       <  1:2674010/150‑1 (MQ=255)
                   tgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCGCTCTAGTTGGAGTGCCCCCTATTCAGGTGGCTgaa                                                                                                                                  >  2:4972718/1‑150 (MQ=255)
                     gCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAAtt                                                                                                                                <  2:5201121/150‑1 (MQ=255)
                       gAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTa                                                                                                                              <  1:4445076/150‑1 (MQ=255)
                        aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCGTAGGGTTGGAGTGCCCCCTGGTCAGGTGGCTAAAGTTAg                                                                                                                             >  1:2456057/1‑150 (MQ=255)
                        aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                                                                             >  1:3883941/1‑150 (MQ=255)
                         ataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGTTGGACAACATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGATGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                                                                                                            <  2:4764483/150‑1 (MQ=255)
                         ataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGc                                                                                                                            <  1:1119602/150‑1 (MQ=255)
                          taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                                           <  1:1157259/150‑1 (MQ=255)
                          taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                                           <  1:1703955/150‑1 (MQ=255)
                           aatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCaa                                                                                                                          <  2:2838310/150‑1 (MQ=255)
                            atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCAGGGTGTGGACTGCCCCCCCGGGAGACGAGCCTGCCCGGGGGGGGGAGTGCCCCCTCTTCGGGGGGCTAACTTTTGCaaa                                                                                                                         >  2:3700727/1‑150 (MQ=255)
                             ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                                        <  2:1312939/150‑1 (MQ=255)
                               aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                                                                                                                      <  1:2809060/150‑1 (MQ=255)
                                 tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTAGAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                                                                                    <  2:3852190/150‑1 (MQ=255)
                                  aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTGGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                                   >  1:1735284/1‑150 (MQ=255)
                                   aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                                                                                                                  <  1:4951456/150‑1 (MQ=255)
                                   aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                                                                                                                  <  2:1759600/150‑1 (MQ=255)
                                     tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                                                <  2:5039621/150‑1 (MQ=255)
                                      atagAGTTGCTTCATGGTGTTAACAGTCAGAGTTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTTCACCCACGGTAGACGATCCTGCCCTATATTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                                                                               <  2:3429802/147‑1 (MQ=255)
                                      aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTCTAGTTGGAGTGCCCCCTATTCTGGTGGCTAAATTTAGCACATCAATACACt                                                                                                               >  2:581603/1‑150 (MQ=255)
                                      aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                                                                               <  1:5061037/150‑1 (MQ=255)
                                      aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                                                                               <  1:3914481/150‑1 (MQ=255)
                                       tCGATTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACtt                                                                                                              <  2:1627315/150‑1 (MQ=255)
                                        cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACtgc                                                                                                             >  2:4461741/1‑148 (MQ=255)
                                             tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                                        <  2:3307457/150‑1 (MQ=255)
                                             tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATGTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCATATAGTTGGAGTGCCCCCGATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                                        >  1:2762330/1‑150 (MQ=255)
                                              gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg                                                                                                       <  1:589445/150‑1 (MQ=255)
                                               ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                                      >  1:497675/1‑150 (MQ=255)
                                                  ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGGTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat                                                                                                   >  2:1099415/1‑150 (MQ=255)
                                                   ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTCAATTTAGCAAATCAATACACTTCAGGGGGGgatt                                                                                                  >  2:4777687/1‑150 (MQ=255)
                                                   ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatt                                                                                                  >  2:4638105/1‑150 (MQ=255)
                                                     ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat                                                                                                >  1:3438884/1‑150 (MQ=255)
                                                       tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg                                                                                              <  1:4370945/150‑1 (MQ=255)
                                                       tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTGAGCAAATCAATACACGTCAGGGGGGTATTATTg                                                                                              >  2:4628171/1‑150 (MQ=255)
                                                         agaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                                                            <  2:3772418/150‑1 (MQ=255)
                                                         agaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                                                            <  1:4127472/150‑1 (MQ=255)
                                                              aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTGGAGTTGGAGTGCCCCCTATTCAGGTGGCTAAGTTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGTGtt                                                                                       >  1:215702/1‑150 (MQ=255)
                                                                 aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                                                    <  1:1349553/150‑1 (MQ=255)
                                                                   gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcccc                                                                                  >  1:448183/1‑150 (MQ=255)
                                                                     ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                                                                                <  1:3114420/150‑1 (MQ=255)
                                                                        taacGAAACGATGGACATCATCCCCCATTGTGTGGACTGCCCCCCCGGTAGACGATCCTGACCTCTAGTTGGAGGTCCCCATATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCTTGGGGGTATTATTGTAGAGTTTCCCCATATg                                                                             <  2:4135816/147‑1 (MQ=255)
                                                                         aCCGATACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTTTTGTAGAGTTTCCCCATATGt                                                                            <  2:466741/150‑1 (MQ=255)
                                                                           cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                                                          <  1:3733333/150‑1 (MQ=255)
                                                                            gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTc                                                                         <  1:4762693/150‑1 (MQ=255)
                                                                            gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTc                                                                         <  1:886162/150‑1 (MQ=255)
                                                                             aaaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCt                                                                        >  1:4210687/1‑150 (MQ=255)
                                                                                cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                                                     <  1:3566202/150‑1 (MQ=255)
                                                                                  aTGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATggg                                                                   >  1:2561937/1‑150 (MQ=255)
                                                                                   tGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAGGGGa                                                                  >  1:3160304/1‑150 (MQ=255)
                                                                                     gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                                                <  1:1525140/150‑1 (MQ=255)
                                                                                     gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                                                <  1:4265817/150‑1 (MQ=255)
                                                                                      aCATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATcc                                                               <  2:32063/150‑1 (MQ=255)
                                                                                       catcatCACCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCa                                                              <  2:1468778/150‑1 (MQ=255)
                                                                                             cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt                                                        <  1:1934446/150‑1 (MQ=255)
                                                                                             cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt                                                        <  1:2521235/150‑1 (MQ=255)
                                                                                             cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt                                                        <  1:5498313/150‑1 (MQ=255)
                                                                                                      tgtgGACTGCCCCCCCGGTAGACGATCCTGCCCTAGAGTTGGAGTGCCCCCTCTTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAg                                               <  1:4808418/150‑1 (MQ=255)
                                                                                                        tgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGa                                             >  2:4851045/1‑150 (MQ=255)
                                                                                                            cggCCCCCCCGGTATACGATCCTGCCCGATATTTGGAGTGTCCCCCATTCAGGTTGCTAACTTTAGCAAATCAATACACTTACGGGGGGTATTAGTGTCGAGTTTCCCCATATGTGTCTATGGGATCCAGGAAATGACAATAAGGAGTTa                                         <  2:686421/148‑1 (MQ=255)
                                                                                                             tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTAc                                        >  2:4290687/1‑150 (MQ=255)
                                                                                                             tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTCCTATGGGATCCAGGAAATGACAATCAGGAGTTAc                                        >  2:1193193/1‑150 (MQ=255)
                                                                                                              gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACa                                       >  2:1585133/1‑150 (MQ=255)
                                                                                                                ccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa                                     >  1:3212389/1‑150 (MQ=255)
                                                                                                                ccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa                                     >  1:351986/1‑150 (MQ=255)
                                                                                                                   cccGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAAc                                  <  2:3899020/150‑1 (MQ=255)
                                                                                                                     cGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACtt                                <  1:3041168/150‑1 (MQ=255)
                                                                                                                     cGGTAGACGATCCGGCCCTATAGTGTGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACtt                                <  2:1026226/150‑1 (MQ=255)
                                                                                                                      ggTAGACGATCCTGCCCTATTGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTa                               <  2:2141441/150‑1 (MQ=255)
                                                                                                                      ggTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTa                               <  1:454596/150‑1 (MQ=255)
                                                                                                                      ggTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTa                               >  2:1969896/1‑150 (MQ=255)
                                                                                                                        tAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaaa                             >  1:3652859/1‑150 (MQ=255)
                                                                                                                        tAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaaa                             <  1:5323590/150‑1 (MQ=255)
                                                                                                                        tAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTaaa                             >  2:756292/1‑150 (MQ=255)
                                                                                                                           tcGATCCTGCACTATTTTTGGAGTGCCCCGTATTCAGGTGGATATATTTATCAAATAAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATGAGGAGTTACAAAAACTTAAATCt                          <  1:2970551/149‑1 (MQ=255)
                                                                                                                             gATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTgg                        >  1:5167197/1‑150 (MQ=255)
                                                                                                                             gATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTgg                        >  1:5430666/1‑150 (MQ=255)
                                                                                                                             gATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTgg                        >  1:5039257/1‑150 (MQ=255)
                                                                                                                              aTCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATAGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGt                       >  2:5156894/1‑150 (MQ=255)
                                                                                                                                ccTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCa                     >  2:861607/1‑150 (MQ=255)
                                                                                                                                 cTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAg                    >  2:4216518/1‑150 (MQ=255)
                                                                                                                                  tGCCCTATACTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAgg                   <  2:3023104/150‑1 (MQ=255)
                                                                                                                                   gCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACGTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAggg                  <  1:1224176/150‑1 (MQ=255)
                                                                                                                                      cTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCaa               <  1:4638240/150‑1 (MQ=255)
                                                                                                                                         taGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAata            <  1:4625329/150‑1 (MQ=255)
                                                                                                                                             tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCga        <  1:2707549/150‑1 (MQ=255)
                                                                                                                                              ggAGTGCCCCCTATTCAGGTGGCTAAATTTAGAAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCgag       <  1:775128/150‑1 (MQ=255)
                                                                                                                                               gAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCgaga      >  2:1894771/1‑150 (MQ=255)
                                                                                                                                               gAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCgaga      <  1:693714/150‑1 (MQ=255)
                                                                                                                                                 gTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACt    >  2:2747470/1‑150 (MQ=255)
                                                                                                                                                   gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTGa  >  1:181033/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTGA  >  CP009273/2759036‑2759332

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: