Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,757,927 A→T 24.0% S196T (TCC→ACC)  yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,757,9270AT24.0% 593.9 / 207.1 271S196T (TCC→ACC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base A (92/114);  new base T (31/34);  total (123/148)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 277 total aligned reads displayed.

GTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGT  >  CP009273/2757778‑2758076
                                                                                                                                                     |                                                                                                                                                     
 tCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGaa                                                                                                                                                      <  2:3026205/150‑1 (MQ=255)
  cATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAAc                                                                                                                                                     >  1:1252938/1‑150 (MQ=255)
   aTCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACt                                                                                                                                                    >  2:197184/1‑150 (MQ=255)
     cTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTga                                                                                                                                                  <  1:832414/150‑1 (MQ=255)
     cTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTga                                                                                                                                                  <  1:832481/150‑1 (MQ=255)
      ttCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTgag                                                                                                                                                 <  2:3749903/150‑1 (MQ=255)
      ttCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTgag                                                                                                                                                 >  1:2245176/1‑150 (MQ=255)
        cGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCGTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGaa                                                                                                                                               <  2:249095/150‑1 (MQ=255)
         gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAAt                                                                                                                                              <  2:669305/150‑1 (MQ=255)
         gACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAAt                                                                                                                                              <  2:331470/150‑1 (MQ=255)
          acacGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATg                                                                                                                                             <  2:1183948/150‑1 (MQ=255)
           cacGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGa                                                                                                                                            >  1:1485476/1‑150 (MQ=255)
                  ctTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATg                                                                                                                                     >  2:55374/1‑150 (MQ=255)
                    tGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTa                                                                                                                                   >  1:3756852/1‑150 (MQ=255)
                      aaCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGa                                                                                                                                 >  1:3835835/1‑150 (MQ=255)
                       aCCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGAt                                                                                                                                <  2:1214535/150‑1 (MQ=255)
                         cACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCt                                                                                                                              <  2:1088105/150‑1 (MQ=255)
                               aGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAg                                                                                                                        >  1:2006311/1‑150 (MQ=255)
                               aGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAg                                                                                                                        >  2:3724978/1‑150 (MQ=255)
                                gTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAgg                                                                                                                       <  2:2670650/150‑1 (MQ=255)
                                   ggTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCa                                                                                                                    <  1:1728246/150‑1 (MQ=255)
                                      cTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGcc                                                                                                                 <  2:2750316/150‑1 (MQ=255)
                                      cTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGcc                                                                                                                 <  2:1897602/150‑1 (MQ=255)
                                      cTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGcc                                                                                                                 >  2:4097670/1‑150 (MQ=255)
                                           cGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATg                                                                                                            <  1:155847/150‑1 (MQ=255)
                                           attAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATg                                                                                                            <  1:2699965/148‑1 (MQ=255)
                                                gTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCtt                                                                                                       <  2:3458056/150‑1 (MQ=255)
                                                 tAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTg                                                                                                      >  2:3743842/1‑150 (MQ=255)
                                                 gaGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCACGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTg                                                                                                      <  2:1381903/149‑1 (MQ=255)
                                                  aGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGa                                                                                                     <  1:3104794/150‑1 (MQ=255)
                                                  aGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGa                                                                                                     >  1:637947/1‑150 (MQ=255)
                                                   ggACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGAt                                                                                                    >  1:1213573/1‑150 (MQ=255)
                                                    gACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATc                                                                                                   >  2:4172652/1‑150 (MQ=255)
                                                     aCTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATcc                                                                                                  >  2:1291446/1‑150 (MQ=255)
                                                       ttGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCt                                                                                                >  2:3418383/1‑150 (MQ=255)
                                                        tGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTa                                                                                               <  1:4274578/150‑1 (MQ=255)
                                                        tGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTa                                                                                               <  1:3264135/150‑1 (MQ=255)
                                                        tGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTa                                                                                               <  2:3062271/150‑1 (MQ=255)
                                                        tGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTa                                                                                               <  2:3587881/150‑1 (MQ=255)
                                                         gTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAg                                                                                              >  1:904667/1‑150 (MQ=255)
                                                           catcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGaa                                                                                            <  1:3151929/150‑1 (MQ=255)
                                                                gCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTcc                                                                                       <  2:2076933/150‑1 (MQ=255)
                                                                   cTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGt                                                                                    >  2:2555274/1‑150 (MQ=255)
                                                                       gAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCGTGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTcgc                                                                                <  1:3650855/150‑1 (MQ=255)
                                                                        aaTGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTcgcg                                                                               <  1:2612772/150‑1 (MQ=255)
                                                                            cGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAgg                                                                           <  1:1330081/150‑1 (MQ=255)
                                                                                gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATg                                                                       >  1:648115/1‑150 (MQ=255)
                                                                                  gaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaa                                                                     >  2:4045348/1‑150 (MQ=255)
                                                                                    aTGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                                                                   >  1:296666/1‑150 (MQ=255)
                                                                                       cTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAacac                                                                >  1:4195300/1‑150 (MQ=255)
                                                                                        tCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACt                                                               >  2:2776394/1‑150 (MQ=255)
                                                                                           cccTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTcc                                                            <  2:2958598/150‑1 (MQ=255)
                                                                                           cccTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTcc                                                            >  2:508132/1‑150 (MQ=255)
                                                                                            ccTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTccc                                                           >  2:4268854/1‑150 (MQ=255)
                                                                                                 cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      >  1:500955/1‑150 (MQ=255)
                                                                                                 cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTa                                                      >  2:3238980/1‑150 (MQ=255)
                                                                                                 cgcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTAGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACTCTTCCCATGTa                                                      >  2:3154636/1‑150 (MQ=255)
                                                                                                  gcgTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAg                                                     <  1:3983024/150‑1 (MQ=255)
                                                                                                    gTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGCAAGCACTTCCCATGTAGcc                                                   >  2:2308510/1‑150 (MQ=255)
                                                                                                    gTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGcc                                                   >  1:1081298/1‑150 (MQ=255)
                                                                                                    gTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGcc                                                   >  2:765214/1‑150 (MQ=255)
                                                                                                    gTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGcc                                                   <  1:23807/150‑1 (MQ=255)
                                                                                                        caGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCAt                                               <  1:3474864/150‑1 (MQ=255)
                                                                                                          gTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCg                                             >  1:978370/1‑150 (MQ=255)
                                                                                                            aTCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGAt                                           <  1:659309/150‑1 (MQ=255)
                                                                                                             tCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATc                                          <  1:3000489/150‑1 (MQ=255)
                                                                                                               aaGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCa                                        >  1:7663/1‑150 (MQ=255)
                                                                                                                aGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCaa                                       <  1:1021663/150‑1 (MQ=255)
                                                                                                                   aTGCGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCt                                    <  2:1635496/150‑1 (MQ=255)
                                                                                                                    tGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTg                                   >  1:1276914/1‑150 (MQ=255)
                                                                                                                      cGGCACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                 >  1:2383050/1‑150 (MQ=255)
                                                                                                                        gCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAAc                               <  2:1671389/150‑1 (MQ=255)
                                                                                                                         cACCAAATGTCTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACg                              >  1:1769690/1‑150 (MQ=255)
                                                                                                                           ccAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCt                            <  2:3326527/150‑1 (MQ=255)
                                                                                                                           ccAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCt                            <  2:1119779/150‑1 (MQ=255)
                                                                                                                            cAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCtt                           <  2:46275/150‑1 (MQ=255)
                                                                                                                             aaaTGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTta                          <  1:4295624/150‑1 (MQ=255)
                                                                                                                             aaaTGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTta                          <  1:3446107/150‑1 (MQ=255)
                                                                                                                               aTGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTtata                        <  2:125313/150‑1 (MQ=255)
                                                                                                                                   cTTCGGCGACGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTg                    <  1:3490869/150‑1 (MQ=255)
                                                                                                                                    ttCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGa                   <  2:1896622/150‑1 (MQ=255)
                                                                                                                                     tcagCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAg                  >  1:2135357/4‑150 (MQ=255)
                                                                                                                                     tCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAg                  <  1:3593830/150‑1 (MQ=255)
                                                                                                                                          gACGGTTTCGGAACTGAGAATGACCCCCTGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGt             <  2:3115042/150‑1 (MQ=255)
                                                                                                                                          gACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGt             >  2:3944561/1‑150 (MQ=255)
                                                                                                                                           aCGGTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGtt            <  2:4289202/150‑1 (MQ=255)
                                                                                                                                             ggTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaa          >  1:3460177/1‑150 (MQ=255)
                                                                                                                                              gTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa         >  1:45000/1‑150 (MQ=255)
                                                                                                                                              gTTTCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa         >  1:2219174/1‑150 (MQ=255)
                                                                                                                                              gTTTCGGACCTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATTTTCTTGATCCCGAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa         <  2:392533/150‑1 (MQ=255)
                                                                                                                                              gTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa         <  2:392101/150‑1 (MQ=255)
                                                                                                                                               tttCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaaa        <  1:4114485/150‑1 (MQ=255)
                                                                                                                                                ttCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaaaa       <  2:3835875/150‑1 (MQ=255)
                                                                                                                                                ttCGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaaaa       >  1:2556521/1‑150 (MQ=255)
                                                                                                                                                  cGGTACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAAt     >  1:696664/1‑150 (MQ=255)
                                                                                                                                                  cGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAAt     <  1:3355010/150‑1 (MQ=255)
                                                                                                                                                  cGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACATGAGTCTGTTAAAAAAt     <  1:55374/150‑1 (MQ=255)
                                                                                                                                                   ggAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAAtt    >  1:1912971/1‑150 (MQ=255)
                                                                                                                                                    gAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTg   <  2:2110027/150‑1 (MQ=255)
                                                                                                                                                     aaCTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGt  >  2:4181280/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGT  >  CP009273/2757778‑2758076

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: