Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,753,076 G→A 11.3% D245N (GAT→AAT)  yfjI → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,753,0760GA11.3% 783.3 / 65.1 248D245N (GAT→AAT) yfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (102/118);  new base A (16/12);  total (118/130)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 260 total aligned reads displayed.

AGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGC  >  CP009273/2752927‑2753225
                                                                                                                                                     |                                                                                                                                                     
   tGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGata                                                                                                                                                    <  1:3688619/150‑1 (MQ=255)
   tGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGata                                                                                                                                                    <  1:3587474/150‑1 (MQ=255)
      cGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATgg                                                                                                                                                 <  2:649172/150‑1 (MQ=255)
      cGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATgg                                                                                                                                                 >  2:3929193/1‑150 (MQ=255)
       gTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGc                                                                                                                                                <  1:3810178/150‑1 (MQ=255)
       gTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGc                                                                                                                                                >  1:825205/1‑150 (MQ=255)
       gTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGc                                                                                                                                                <  1:2555361/150‑1 (MQ=255)
       gTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGc                                                                                                                                                >  2:1900968/1‑150 (MQ=255)
         ttATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCaa                                                                                                                                              <  1:4337504/150‑1 (MQ=255)
               aTAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGgga                                                                                                                                        >  1:2656499/1‑150 (MQ=255)
                     tgtgGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGggattggat                                                                                                                                  >  2:1279917/1‑150 (MQ=255)
                          gATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGAtttttt                                                                                                                             <  1:4546745/150‑1 (MQ=255)
                           aTGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTg                                                                                                                            <  2:1691511/150‑1 (MQ=255)
                            tGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGc                                                                                                                           >  2:3888113/1‑150 (MQ=255)
                                   atatTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGc                                                                                                                    >  2:863961/1‑150 (MQ=255)
                                    tatTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGcc                                                                                                                   <  2:2968358/150‑1 (MQ=255)
                                      tttACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTc                                                                                                                 <  1:2921467/150‑1 (MQ=255)
                                       ttACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCa                                                                                                                <  1:1239711/150‑1 (MQ=255)
                                       ttACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCa                                                                                                                >  1:1501784/1‑150 (MQ=255)
                                        tACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCAt                                                                                                               >  1:848414/1‑150 (MQ=255)
                                         aCGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCAtg                                                                                                              >  2:1144033/1‑150 (MQ=255)
                                         aCGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCAtg                                                                                                              >  2:3174439/1‑150 (MQ=255)
                                             tGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGcc                                                                                                          <  1:229269/150‑1 (MQ=255)
                                              ggaaaggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCa                                                                                                         <  1:2772736/150‑1 (MQ=255)
                                                 aaggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGcc                                                                                                      <  1:1178567/150‑1 (MQ=255)
                                                  aggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCt                                                                                                     <  1:2727135/150‑1 (MQ=255)
                                                  aggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCt                                                                                                     <  1:744941/150‑1 (MQ=255)
                                                  aggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCt                                                                                                     <  2:3204642/150‑1 (MQ=255)
                                                  aggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCt                                                                                                     <  2:3112506/150‑1 (MQ=255)
                                                   ggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTg                                                                                                    <  1:4565055/150‑1 (MQ=255)
                                                     aaaaaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCt                                                                                                  >  2:1472667/1‑150 (MQ=255)
                                                      aaaaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCtt                                                                                                 >  1:1476323/1‑150 (MQ=255)
                                                       aaaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTc                                                                                                <  2:2205002/150‑1 (MQ=255)
                                                         aCGAGCCCGAGAAATTAATTAGAGCTGCAAGAATATCACTGTCGCGGATGGTCCAGCCTAATGTTTTTAAGGGTGATATCGACAGGAAAGGAGATATGGCATAGGGGATTGGATTTTTTGCTCGTTGCCTCATGTGCCAGCCTGCTTCaa                                                                                              <  2:4427468/150‑1 (MQ=255)
                                                          cGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAac                                                                                             >  1:168202/1‑150 (MQ=255)
                                                                cGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAgg                                                                                       >  2:1237563/1‑150 (MQ=255)
                                                                    aaaTTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAAc                                                                                   <  1:1053413/150‑1 (MQ=255)
                                                                     aattaattAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACa                                                                                  <  2:3604/150‑1 (MQ=255)
                                                                     aattaattAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACa                                                                                  >  2:275303/1‑150 (MQ=255)
                                                                     aattaattAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACa                                                                                  >  2:275212/1‑150 (MQ=255)
                                                                      attaattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAg                                                                                 >  2:868701/1‑150 (MQ=255)
                                                                         aattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaa                                                                              <  1:1611190/150‑1 (MQ=255)
                                                                         aattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaa                                                                              >  2:2442447/1‑150 (MQ=255)
                                                                         aattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaa                                                                              <  1:4176069/150‑1 (MQ=255)
                                                                          attaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaaa                                                                             >  2:4374156/1‑150 (MQ=255)
                                                                            taGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAAt                                                                           >  2:331971/1‑150 (MQ=255)
                                                                            taGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAAt                                                                           <  1:1939751/150‑1 (MQ=255)
                                                                             agagATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAAtt                                                                          <  1:2339134/150‑1 (MQ=255)
                                                                              gagaTGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAAttt                                                                         <  2:1544262/150‑1 (MQ=255)
                                                                              gagaTGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTTACAGAAAAAttt                                                                         <  2:3087479/150‑1 (MQ=255)
                                                                                  tGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCaa                                                                     <  1:3866190/150‑1 (MQ=255)
                                                                                  tGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCaa                                                                     >  2:922531/1‑150 (MQ=255)
                                                                                  tGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCaa                                                                     >  1:1042370/1‑150 (MQ=255)
                                                                                  tGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCaa                                                                     >  2:2575059/1‑150 (MQ=255)
                                                                                    cAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAcc                                                                   >  1:127361/1‑150 (MQ=255)
                                                                                    cAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAccc                                                                  <  1:566803/150‑1 (MQ=255)
                                                                                        aaTAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAt                                                               >  2:2924831/1‑150 (MQ=255)
                                                                                          tAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAttt                                                             >  1:4065485/1‑150 (MQ=255)
                                                                                           aaCACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtttt                                                            >  2:3118635/1‑150 (MQ=255)
                                                                                           aaCACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtttt                                                            >  1:778146/1‑150 (MQ=255)
                                                                                             cacTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtttttt                                                          >  2:3477228/1‑150 (MQ=255)
                                                                                              acTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTc                                                         <  1:2299049/150‑1 (MQ=255)
                                                                                              acTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAATTTTTc                                                         >  1:1756126/1‑150 (MQ=255)
                                                                                                   cGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATg                                                    >  2:1366999/1‑150 (MQ=255)
                                                                                                    gCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGa                                                   <  1:1558533/150‑1 (MQ=255)
                                                                                                    gCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGa                                                   <  1:1820155/150‑1 (MQ=255)
                                                                                                     cTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGaa                                                  >  2:4574872/1‑150 (MQ=255)
                                                                                                       gATGGTCCAGCCTAATGTTTTTAAGGGTCATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACa                                                <  2:3406047/150‑1 (MQ=255)
                                                                                                        aTGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCAAACCCAATTTTTTCAAATGAACAt                                               <  1:1894519/150‑1 (MQ=255)
                                                                                                        aTGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAt                                               <  2:2708386/150‑1 (MQ=255)
                                                                                                         tGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAtt                                              >  1:1094895/1‑150 (MQ=255)
                                                                                                          ggTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAttt                                             <  2:2335719/150‑1 (MQ=255)
                                                                                                          ggTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAttt                                             <  1:4402963/150‑1 (MQ=255)
                                                                                                                gCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGt                                       <  2:4184050/150‑1 (MQ=255)
                                                                                                                gCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGt                                       >  1:1699877/1‑150 (MQ=255)
                                                                                                                gCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGt                                       >  2:2010532/1‑150 (MQ=255)
                                                                                                                 ccTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTa                                      >  1:364289/1‑150 (MQ=255)
                                                                                                                 ccTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTa                                      <  2:1737116/150‑1 (MQ=255)
                                                                                                                    aaTGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAttt                                   >  1:2016282/1‑150 (MQ=255)
                                                                                                                           ttAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCTCGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGTTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAAc                            <  2:393789/150‑1 (MQ=255)
                                                                                                                           ttAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAAc                            <  1:3006850/150‑1 (MQ=255)
                                                                                                                            tAAGGGTTATATCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACg                           <  1:2316695/150‑1 (MQ=255)
                                                                                                                                 gTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTa                      <  1:2894172/150‑1 (MQ=255)
                                                                                                                                     tatCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGa                  >  2:3415295/1‑150 (MQ=255)
                                                                                                                                      atCGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaa                 >  1:201179/1‑150 (MQ=255)
                                                                                                                                       tCGACAGGAAAGCAGATATGGCAAAGGGGAATGGATTTTTTGCACGCTGCCTCAAGTGCCAGCCAGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTAATCAAATGAACATATGACGGTAATTAACAAACATCTAAAGGaaa                >  2:4114450/1‑150 (MQ=255)
                                                                                                                                        cGACAGGAAAGGAAATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAAt               >  2:71914/1‑150 (MQ=255)
                                                                                                                                         gACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAAtt              >  1:3548694/1‑150 (MQ=255)
                                                                                                                                          aCAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACTGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTg             <  1:3187221/150‑1 (MQ=255)
                                                                                                                                          aCAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTg             <  2:3797995/150‑1 (MQ=255)
                                                                                                                                            aGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGtt           >  1:420479/1‑150 (MQ=255)
                                                                                                                                            aGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGtt           >  1:4327749/1‑150 (MQ=255)
                                                                                                                                              gAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTCTTTCACCAACGTCTTATGGAAATTGTTaa         >  2:3665690/1‑150 (MQ=255)
                                                                                                                                               aaaGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAAt        >  1:4110287/1‑150 (MQ=255)
                                                                                                                                                aaGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATg       <  1:1152428/150‑1 (MQ=255)
                                                                                                                                                 aGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATga      >  1:3499802/1‑150 (MQ=255)
                                                                                                                                                  ggAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATgag     >  2:202945/1‑150 (MQ=255)
                                                                                                                                                   gagaTATGGCAAAGGGGATTGGATTTTTTGTACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATgaga    <  1:3916839/150‑1 (MQ=255)
                                                                                                                                                     gaTATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGc  <  1:2129139/150‑1 (MQ=255)
                                                                                                                                                     gaTATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGc  >  2:3821297/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGC  >  CP009273/2752927‑2753225

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: