Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 569,073 A→T 8.3% T83S (ACT→TCT)  rusA → DLP12 prophage; endonuclease RUS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273569,0730AT 8.3% 994.5 / 63.3 314T83S (ACT→TCT) rusADLP12 prophage; endonuclease RUS
Reads supporting (aligned to +/- strand):  ref base A (138/150);  new base T (14/12);  total (152/162)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 325 total aligned reads displayed.

TACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGC  >  CP009273/568926‑569221
                                                                                                                                                    |                                                                                                                                                    
 aCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGcactcact                                                                                                                                                    >  2:3178939/1‑149 (MQ=255)
    gcgATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTaaa                                                                                                                                                 <  2:4178012/149‑1 (MQ=255)
      gATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGc                                                                                                                                               >  2:4495148/1‑149 (MQ=255)
          aCGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGt                                                                                                                                           <  2:1644792/149‑1 (MQ=255)
           cGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGtt                                                                                                                                          <  1:3261046/149‑1 (MQ=255)
             tCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTc                                                                                                                                        <  1:4297496/149‑1 (MQ=255)
              cGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTct                                                                                                                                       <  1:4676079/149‑1 (MQ=255)
              cGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTct                                                                                                                                       <  1:2842924/149‑1 (MQ=255)
              cGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTct                                                                                                                                       <  1:2278706/149‑1 (MQ=255)
                 ccGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTctggc                                                                                                                                    <  1:1545733/149‑1 (MQ=255)
                     aTCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGga                                                                                                                                >  2:3017513/1‑149 (MQ=255)
                     aTCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGga                                                                                                                                >  1:3029698/1‑149 (MQ=255)
                        aTTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGgatga                                                                                                                             <  2:1928379/149‑1 (MQ=255)
                        aTTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGgatga                                                                                                                             <  2:1928656/149‑1 (MQ=255)
                         ttAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGgatgat                                                                                                                            <  1:4463451/149‑1 (MQ=255)
                          tAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGatgatg                                                                                                                           >  2:4429933/1‑149 (MQ=255)
                            aaaaCGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCt                                                                                                                         >  1:3405659/1‑149 (MQ=255)
                               aCGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAg                                                                                                                      >  2:257349/1‑149 (MQ=255)
                                 gCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGt                                                                                                                    <  1:61042/149‑1 (MQ=255)
                                    aTGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGt                                                                                                                 >  1:4286428/1‑149 (MQ=255)
                                     ttcTGGATATCGTCCTGGCTATGCCTGTGAAAATCCTCATTGAGTGCCTCATGCCGGATCGCCGTCGCCGTGTCCTGGATAATCTGCAAAAAGGCTCTTTTGACTCACTCTCTAAATCAGGTTTCTGGCTGGATGATGCTCAGGTCGtt                                                                                                                <  1:1768108/147‑1 (MQ=255)
                                      gCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTg                                                                                                               >  1:4320281/1‑149 (MQ=255)
                                        tGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGAt                                                                                                             >  1:1841318/1‑149 (MQ=255)
                                              tCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcg                                                                                                        <  1:3708819/149‑1 (MQ=255)
                                               cGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcgcg                                                                                                      <  2:2859177/149‑1 (MQ=255)
                                               cGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcgcg                                                                                                      >  2:1872868/1‑149 (MQ=255)
                                               cGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACcgcg                                                                                                      <  2:466996/149‑1 (MQ=255)
                                                ggCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGt                                                                                                     >  1:4495886/1‑149 (MQ=255)
                                                   cTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTtgt                                                                                                  <  1:2039982/149‑1 (MQ=255)
                                                    tGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTtgtg                                                                                                 >  1:198461/1‑149 (MQ=255)
                                                     ggCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGa                                                                                                >  2:1413118/1‑149 (MQ=255)
                                                      gCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGaa                                                                                               >  2:2378685/1‑149 (MQ=255)
                                                           gCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGc                                                                                          >  1:948403/1‑149 (MQ=255)
                                                            ccTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGcc                                                                                         >  2:649467/1‑149 (MQ=255)
                                                             cTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCt                                                                                        <  1:1661817/149‑1 (MQ=255)
                                                              tgtgAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTg                                                                                       >  2:1707951/1‑149 (MQ=255)
                                                                tgAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGtt                                                                                     >  1:4168823/1‑149 (MQ=255)
                                                                 gAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTa                                                                                    >  2:248509/1‑149 (MQ=255)
                                                                 gAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTa                                                                                    >  2:375676/1‑149 (MQ=255)
                                                                  aaaaTCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTAc                                                                                   >  1:412987/1‑149 (MQ=255)
                                                                   aaaTCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTAcc                                                                                  >  1:1734208/1‑149 (MQ=255)
                                                                        cGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAgg                                                                             >  1:102238/1‑149 (MQ=255)
                                                                         gCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGt                                                                            >  1:1642276/1‑149 (MQ=255)
                                                                         gCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGt                                                                            <  2:1109615/149‑1 (MQ=255)
                                                                          cATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTg                                                                           <  1:4056795/149‑1 (MQ=255)
                                                                           aTTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTgg                                                                          <  1:3201205/149‑1 (MQ=255)
                                                                           aTTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTgg                                                                          <  1:3493692/149‑1 (MQ=255)
                                                                              gAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAg                                                                       <  1:289330/149‑1 (MQ=255)
                                                                                gTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGc                                                                     <  2:463145/149‑1 (MQ=255)
                                                                                gTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGc                                                                     <  2:3476740/149‑1 (MQ=255)
                                                                                 tGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCt                                                                    <  1:4638032/149‑1 (MQ=255)
                                                                                    cacaTGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGa                                                                 >  2:756876/1‑149 (MQ=255)
                                                                                     acaTGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGaa                                                                <  1:3417546/149‑1 (MQ=255)
                                                                                     acaTGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGaa                                                                <  2:3626702/149‑1 (MQ=255)
                                                                                     acaTGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGaa                                                                <  1:120464/149‑1 (MQ=255)
                                                                                        tGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTg                                                             <  1:88888/149‑1 (MQ=255)
                                                                                            ggATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAACGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCa                                                         <  2:4578784/149‑1 (MQ=255)
                                                                                             gATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCAt                                                        >  2:1477743/1‑149 (MQ=255)
                                                                                             gATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCAt                                                        >  2:3799540/1‑149 (MQ=255)
                                                                                             gATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCAt                                                        >  2:3867201/1‑149 (MQ=255)
                                                                                              aTCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATc                                                       >  2:212250/1‑149 (MQ=255)
                                                                                              aTCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATc                                                       <  1:1098117/149‑1 (MQ=255)
                                                                                                 gccgtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCAcc                                                    <  2:1743514/149‑1 (MQ=255)
                                                                                                    gtcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGaa                                                 >  2:1143387/1‑149 (MQ=255)
                                                                                                     tcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGaaa                                                >  1:3809122/1‑149 (MQ=255)
                                                                                                     tcgccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGaaa                                                <  1:2637195/149‑1 (MQ=255)
                                                                                                       gccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATg                                              >  1:1691388/1‑149 (MQ=255)
                                                                                                       gccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATg                                              <  1:3216962/149‑1 (MQ=255)
                                                                                                        ccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATgg                                             <  2:1805415/149‑1 (MQ=255)
                                                                                                        ccgtGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATgg                                             <  2:4031689/149‑1 (MQ=255)
                                                                                                           tGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGga                                          <  1:3832002/149‑1 (MQ=255)
                                                                                                            gACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaa                                         <  1:681376/149‑1 (MQ=255)
                                                                                                             aCCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaat                                        <  1:777244/149‑1 (MQ=255)
                                                                                                              ccTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatg                                       >  1:1149234/1‑149 (MQ=255)
                                                                                                               cTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatga                                      >  2:527842/1‑149 (MQ=255)
                                                                                                               cTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatga                                      <  2:3628263/149‑1 (MQ=255)
                                                                                                                tGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCTCTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGgaatgaa                                     >  1:431429/1‑149 (MQ=255)
                                                                                                                   aTACTCTGCAAAAAGCCGCTTTTGACGCTCTCACTTAAGCAGGTTTCTTGCTGGATGTTGCTCAGGTCGTTGATTACCGCGTTGTGAAGTTGCCTGTTACCATTGGTGGGAGGCTGGTACTGAACATCACCGAAATGGGGAatgaatga                                  <  2:4714505/149‑1 (MQ=255)
                                                                                                                    tAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAatgat                                 >  1:1073037/1‑149 (MQ=255)
                                                                                                                      aTCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGt                               <  2:3444567/149‑1 (MQ=255)
                                                                                                                      aTCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGt                               >  1:3772706/1‑149 (MQ=255)
                                                                                                                      aTCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGt                               >  1:4064374/1‑149 (MQ=255)
                                                                                                                       tCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGtt                              >  2:1317027/1‑149 (MQ=255)
                                                                                                                       tCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGtt                              >  2:3288965/1‑149 (MQ=255)
                                                                                                                          gCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGa                           >  2:2129412/1‑149 (MQ=255)
                                                                                                                           cAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAg                          <  1:3651790/149‑1 (MQ=255)
                                                                                                                             aaaaaGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGTTCGTTGATTACAGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATCGGGAATGAATGATGTTTGAGt                         >  1:2581555/1‑149 (MQ=255)
                                                                                                                               aaaGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGttt                       <  1:1502579/149‑1 (MQ=255)
                                                                                                                               aaaGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGttt                       >  1:3934094/1‑149 (MQ=255)
                                                                                                                                    cGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAatat                  <  1:4528527/149‑1 (MQ=255)
                                                                                                                                     gCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATg                 >  2:2884743/1‑149 (MQ=255)
                                                                                                                                     gCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATg                 <  2:202711/149‑1 (MQ=255)
                                                                                                                                      cTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATgg                <  2:4287471/149‑1 (MQ=255)
                                                                                                                                       ttttGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGc               <  1:3349820/149‑1 (MQ=255)
                                                                                                                                           gACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGaa           <  2:2635560/149‑1 (MQ=255)
                                                                                                                                           gACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGaa           >  2:2944722/1‑149 (MQ=255)
                                                                                                                                           gACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGaa           <  2:322948/149‑1 (MQ=255)
                                                                                                                                           gACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGaa           >  2:941128/1‑149 (MQ=255)
                                                                                                                                            aCGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGACTGATGTTTGAGTTTAATATGGCAGAAc          >  1:1886266/1‑149 (MQ=255)
                                                                                                                                                 ctcactAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAActtctt     <  1:4156399/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGCCGTCGCCGTGACCTGGATAATCTGCAAAAAGCCGCTTTTGACGCACTCACTAAAGCAGGTTTCTGGCTGGATGATGCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGC  >  CP009273/568926‑569221

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: