Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,759,184 C→T 6.0% 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,1840CT 6.0% 997.8 / 34.2 282intergenic (‑49/‑93)yfjM/rnlACP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base C (114/151);  new base T (12/5);  total (126/156)
Fisher's exact test for biased strand distribution p-value = 4.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 292 total aligned reads displayed.

AGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTG  >  CP009273/2759036‑2759331
                                                                                                                                                     |                                                                                                                                                   
    ttttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGGTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccccc                                                                                                                                                 <  2:2274983/149‑1 (MQ=255)
     tttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                                                                                                                <  2:1095645/149‑1 (MQ=255)
     tttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCt                                                                                                                                                <  2:2601572/149‑1 (MQ=255)
       tttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATGGGTGGTGTGCCCCCTAt                                                                                                                                              >  1:612744/1‑149 (MQ=255)
       tttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTAt                                                                                                                                              <  1:4716037/149‑1 (MQ=255)
        ttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCGCTATAGTTGGAGTGCCCCCTCtt                                                                                                                                             >  2:100489/1‑149 (MQ=255)
         tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTCGAGTTGGAGTGCCCCCTATTc                                                                                                                                            >  1:1857244/1‑149 (MQ=255)
         tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTc                                                                                                                                            <  1:3157502/149‑1 (MQ=255)
           ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTGTGTGGACTGCCCCCCCGGGCGACGATCCTGCCCTAGAGGTGGCGTGCCCCCtcggccg                                                                                                                                          >  2:2854654/1‑143 (MQ=255)
             gTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGt                                                                                                                                        <  1:2994042/149‑1 (MQ=255)
                tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATAATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCGATTCAGGTGGc                                                                                                                                     >  1:1281354/1‑149 (MQ=255)
                  ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTCGACGATCCTGCCCGCTCGTTGGAGTGCCCCCTAGTCAGGTGGCTa                                                                                                                                   >  2:2349378/1‑149 (MQ=255)
                  ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTa                                                                                                                                   <  1:1021756/149‑1 (MQ=255)
                  ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCGATAGTTGGAGTGCCCCCTATTCAGGTGGCTa                                                                                                                                   >  2:2585585/1‑149 (MQ=255)
                   tgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaa                                                                                                                                  >  1:4263341/1‑149 (MQ=255)
                    ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTaaa                                                                                                                                 <  1:2529206/149‑1 (MQ=255)
                     gCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAt                                                                                                                                <  1:4386039/149‑1 (MQ=255)
                      cGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAtt                                                                                                                               <  1:2497502/149‑1 (MQ=255)
                        aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATCGTTGGAGTGCCCCCTATTCAGGTGGCTCAATTTa                                                                                                                             >  1:3638847/1‑149 (MQ=255)
                           aatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATGGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                                          >  1:1062467/1‑149 (MQ=255)
                           aatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                                          <  2:5170448/149‑1 (MQ=255)
                            atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATGTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGGTGGAGTGCCCCCTCTTCAGGTGGCTAAATTTAGCaa                                                                                                                         >  1:2980670/1‑149 (MQ=255)
                               aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCGCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                                                                                                                      >  1:646440/1‑149 (MQ=255)
                               aCTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                                                                                                                      >  1:5237017/1‑149 (MQ=255)
                                  aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                                                                                   >  1:1446522/1‑149 (MQ=255)
                                    gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                                                                                                                 <  2:3508834/149‑1 (MQ=255)
                                     tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                                                                                <  1:3660159/149‑1 (MQ=255)
                                      aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                                               <  1:4560984/149‑1 (MQ=255)
                                      aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                                               >  2:391479/1‑149 (MQ=255)
                                       tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                                                                              <  1:2700292/149‑1 (MQ=255)
                                          aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATCGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa                                                                                                           >  2:2029670/1‑149 (MQ=255)
                                           gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                                          >  1:958597/1‑149 (MQ=255)
                                           gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                                          <  2:206147/149‑1 (MQ=255)
                                             tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                                        <  1:2022523/149‑1 (MQ=255)
                                             tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                                        <  2:4157002/149‑1 (MQ=255)
                                              gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                                       <  1:1760626/149‑1 (MQ=255)
                                              gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                                       <  1:3253855/149‑1 (MQ=255)
                                               ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                                       <  2:4469142/149‑1 (MQ=255)
                                                tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                                     <  2:1684253/149‑1 (MQ=255)
                                                tgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                                     >  1:3329620/1‑149 (MQ=255)
                                                 gctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGt                                                                                                    <  2:1130866/149‑1 (MQ=255)
                                                  ctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                                   <  2:3314547/149‑1 (MQ=255)
                                                   ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCCTCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat                                                                                                  <  2:777125/149‑1 (MQ=255)
                                                    tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCGATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatt                                                                                                 >  1:2651381/1‑149 (MQ=255)
                                                     ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta                                                                                                <  2:5087/149‑1 (MQ=255)
                                                     ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtatta                                                                                                >  2:3524062/1‑149 (MQ=255)
                                                      gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat                                                                                               <  1:2574195/149‑1 (MQ=255)
                                                      gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat                                                                                               >  2:267469/1‑149 (MQ=255)
                                                       tGAGATCAGGAAGAGGTTACCGACACGATGGACATCATCCCCCAGTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                                              <  2:606277/149‑1 (MQ=255)
                                                       tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                                              <  2:5216584/149‑1 (MQ=255)
                                                       tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                                              <  2:3331432/149‑1 (MQ=255)
                                                        gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg                                                                                             <  2:67879/149‑1 (MQ=255)
                                                          gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGGGCCCCCTATTCAGGTGGCTAAATTTCGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                                                           <  1:4817733/149‑1 (MQ=255)
                                                           aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTCTTGTag                                                                                          >  1:1486168/1‑149 (MQ=255)
                                                           aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag                                                                                          <  1:4066168/149‑1 (MQ=255)
                                                            tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCGATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTCTGATTGTaga                                                                                         >  1:4135722/1‑149 (MQ=255)
                                                                gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGGAGACGATCCTGCCCGATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAAGTTAGCAAATCAATACACTGCAGGGGGGGATTATTGTAGAGttt                                                                                     >  2:632252/1‑149 (MQ=255)
                                                                 aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTc                                                                                    <  1:501438/149‑1 (MQ=255)
                                                                     ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                                                                                <  2:1910405/149‑1 (MQ=255)
                                                                       ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCata                                                                              <  1:1589999/149‑1 (MQ=255)
                                                                         aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATg                                                                            >  2:1611676/1‑149 (MQ=255)
                                                                           cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGtt                                                                          <  2:3877760/149‑1 (MQ=255)
                                                                               aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTa                                                                      <  1:826229/149‑1 (MQ=255)
                                                                                    ggACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTCGAGTTTCCCCATATGTTTCTATGGGa                                                                 >  2:329788/1‑149 (MQ=255)
                                                                                     gaaaacatcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGATAAATTTAGCAAATAAATACTCTTCAGGTGTGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAt                                                                <  1:1836870/144‑1 (MQ=255)
                                                                                     gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACCCGTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAt                                                                >  2:4533098/1‑149 (MQ=255)
                                                                                     gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAt                                                                <  2:5034858/149‑1 (MQ=255)
                                                                                     gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAt                                                                <  1:549997/149‑1 (MQ=255)
                                                                                      aCATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                                               <  2:3970207/149‑1 (MQ=255)
                                                                                      aCATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                                               <  2:4025202/149‑1 (MQ=255)
                                                                                       catcatCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATcc                                                              <  1:2372935/149‑1 (MQ=255)
                                                                                         tcatcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAg                                                            <  1:721870/149‑1 (MQ=255)
                                                                                          catcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAgg                                                           <  1:967696/149‑1 (MQ=255)
                                                                                            tcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaa                                                         >  1:1655696/1‑149 (MQ=255)
                                                                                             cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTAGAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACATTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaaa                                                        <  2:4263341/149‑1 (MQ=255)
                                                                                                  aTTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACa                                                   >  1:3332339/1‑149 (MQ=255)
                                                                                                          gACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAg                                           <  1:3307023/149‑1 (MQ=255)
                                                                                                           aCTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATACAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGt                                          <  2:2713236/149‑1 (MQ=255)
                                                                                                            cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGtt                                         <  2:1357861/149‑1 (MQ=255)
                                                                                                                    cccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaaaa                                  >  1:4677079/1‑149 (MQ=255)
                                                                                                                      cGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt                                >  2:795549/1‑149 (MQ=255)
                                                                                                                      cGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt                                <  1:3552344/149‑1 (MQ=255)
                                                                                                                       ggTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACtt                               >  1:4456705/1‑149 (MQ=255)
                                                                                                                        gTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTa                              >  2:3631562/1‑149 (MQ=255)
                                                                                                                           gACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAAt                           >  2:3205526/1‑149 (MQ=255)
                                                                                                                           gACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAAt                           >  1:3179506/1‑149 (MQ=255)
                                                                                                                              gATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTg                        <  2:2393519/149‑1 (MQ=255)
                                                                                                                              gATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTg                        >  2:2781995/1‑149 (MQ=255)
                                                                                                                               aTCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTgg                       >  1:3086853/1‑149 (MQ=255)
                                                                                                                               aTCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTgg                       >  1:2796179/1‑149 (MQ=255)
                                                                                                                                 ccTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc                     >  2:951781/1‑149 (MQ=255)
                                                                                                                                         ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCaaa             >  2:2977505/1‑149 (MQ=255)
                                                                                                                                          taGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAat            >  2:2207707/1‑149 (MQ=255)
                                                                                                                                           aGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATGATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAata           <  2:2891348/149‑1 (MQ=255)
                                                                                                                                             ttGGAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATc         >  2:1568567/1‑149 (MQ=255)
                                                                                                                                               ggAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCga       <  1:5165182/149‑1 (MQ=255)
                                                                                                                                                gAGTGTCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCgag      >  1:3824100/1‑149 (MQ=255)
                                                                                                                                                gAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCGGGCAAATATCgaga     >  2:1242272/1‑149 (MQ=255)
                                                                                                                                                  gTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGAc    >  1:3416181/1‑149 (MQ=255)
                                                                                                                                                    gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTg  <  2:4357845/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCAAATATCGAGACTG  >  CP009273/2759036‑2759331

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: