Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,750,424 G→T 48.1% E303* (GAA→TAA)  intA → CP4‑57 prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,750,4240GT48.1% 153.8 / inf 234E303* (GAA→TAA) intACP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base G (63/58);  new base T (60/53);  total (123/111)
Fisher's exact test for biased strand distribution p-value = 8.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

CTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTT  >  CP009273/2750275‑2750573
                                                                                                                                                     |                                                                                                                                                     
 tgAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCgaga                                                                                                                                                      <  1:2299288/150‑1 (MQ=255)
    ggCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAAtt                                                                                                                                                   >  2:3540337/1‑150 (MQ=255)
    ggCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAAtt                                                                                                                                                   >  1:1907121/1‑150 (MQ=255)
      cTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTa                                                                                                                                                 <  2:3817405/150‑1 (MQ=255)
       tCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTAt                                                                                                                                                >  2:253789/1‑150 (MQ=255)
        cGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATc                                                                                                                                               >  1:92340/1‑150 (MQ=255)
         gATGGGATGAGATCGATTTCAATGCGAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCt                                                                                                                                              <  2:1187931/150‑1 (MQ=255)
           tGGGATGAGATCGATTTCAATGCTAGCGAATTGAAAATTCCTGCAGCTCGAATGAAGATGAACCGTGACCATACGGTTCCACTATATGATGGGGCTCTTGCTATTCTGTAAATGATTAATCCTCTCAGTGGTGGCCGATAATTTATCttt                                                                                                                                            <  2:1954374/150‑1 (MQ=255)
            gggATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTc                                                                                                                                           >  1:1819732/1‑150 (MQ=255)
            gggATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTc                                                                                                                                           >  2:2444112/1‑150 (MQ=255)
             cgaagagaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAAAGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTAGTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATATATCTATcc                                                                                                                                          >  2:1385722/5‑150 (MQ=255)
              gATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCt                                                                                                                                         <  2:2922192/150‑1 (MQ=255)
              gATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCt                                                                                                                                         <  1:3213230/150‑1 (MQ=255)
              gATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCt                                                                                                                                         <  2:1017370/150‑1 (MQ=255)
               aTGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTa                                                                                                                                        <  2:1247695/150‑1 (MQ=255)
                 gagaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGc                                                                                                                                      >  1:1124320/1‑150 (MQ=255)
                 gagaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGc                                                                                                                                      <  1:684343/150‑1 (MQ=255)
                  anaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGcc                                                                                                                                     >  2:4332860/1‑150 (MQ=255)
                  agaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGcc                                                                                                                                     >  2:3893096/1‑150 (MQ=255)
                   gaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCg                                                                                                                                    >  2:1989157/1‑150 (MQ=255)
                    aTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGt                                                                                                                                   >  1:4163607/1‑150 (MQ=255)
                       gATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATc                                                                                                                                >  1:1531981/1‑150 (MQ=255)
                        aTTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCa                                                                                                                               <  1:3753373/150‑1 (MQ=255)
                           tCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGc                                                                                                                            <  2:617002/150‑1 (MQ=255)
                            cAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGcc                                                                                                                           >  1:1427456/1‑150 (MQ=255)
                            cAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGcc                                                                                                                           >  2:445910/1‑150 (MQ=255)
                             aaTGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGccc                                                                                                                          <  2:2589605/150‑1 (MQ=255)
                               tGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaa                                                                                                                        >  1:1028203/1‑150 (MQ=255)
                                gCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCccaac                                                                                                                       <  2:2689396/150‑1 (MQ=255)
                                  tAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCccaacca                                                                                                                     <  1:666786/150‑1 (MQ=255)
                                  tAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCccaacca                                                                                                                     >  2:2199590/1‑150 (MQ=255)
                                     cGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaacc                                                                                                                  <  2:106798/150‑1 (MQ=255)
                                      gAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaacca                                                                                                                 >  2:4350127/1‑150 (MQ=255)
                                         tGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATg                                                                                                              >  2:2126613/1‑150 (MQ=255)
                                          ggAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGa                                                                                                             >  1:1887735/1‑150 (MQ=255)
                                               aTTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGc                                                                                                        >  2:2020862/1‑150 (MQ=255)
                                                  ccTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCaa                                                                                                     >  2:395913/1‑150 (MQ=255)
                                                  ccTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCaa                                                                                                     <  1:2616394/150‑1 (MQ=255)
                                                   cTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACCATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCaaa                                                                                                    >  2:275776/1‑150 (MQ=255)
                                                     gCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACa                                                                                                  <  1:2989099/150‑1 (MQ=255)
                                                       aGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGt                                                                                                <  2:3766070/150‑1 (MQ=255)
                                                       aGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGt                                                                                                <  1:4612051/150‑1 (MQ=255)
                                                        gCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTg                                                                                               >  1:1203013/1‑150 (MQ=255)
                                                            gAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATATCCAAACAGTGAATg                                                                                           >  1:2067217/1‑150 (MQ=255)
                                                               tGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcag                                                                                        >  2:4042656/1‑150 (MQ=255)
                                                                gAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcagc                                                                                       >  1:3849237/1‑150 (MQ=255)
                                                                gAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcagc                                                                                       <  2:812282/150‑1 (MQ=255)
                                                                gAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcagc                                                                                       >  1:2457033/1‑150 (MQ=255)
                                                                 aaGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcagca                                                                                      <  2:1539962/150‑1 (MQ=255)
                                                                        aCCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGc                                                                               <  1:2924717/150‑1 (MQ=255)
                                                                         ccGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCg                                                                              >  2:3873016/1‑150 (MQ=255)
                                                                             gACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCt                                                                          >  1:821598/1‑150 (MQ=255)
                                                                             gACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCt                                                                          >  2:579116/1‑150 (MQ=255)
                                                                             gACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCt                                                                          >  1:123833/1‑150 (MQ=255)
                                                                              aCCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTg                                                                         <  1:4440945/150‑1 (MQ=255)
                                                                               ccATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTgg                                                                        <  1:698716/150‑1 (MQ=255)
                                                                               ccATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTgg                                                                        <  1:1009836/150‑1 (MQ=255)
                                                                                cATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGc                                                                       <  2:3067279/150‑1 (MQ=255)
                                                                                cATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGc                                                                       <  1:762569/150‑1 (MQ=255)
                                                                                cATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGc                                                                       >  1:4309443/1‑150 (MQ=255)
                                                                                cATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGc                                                                       <  2:4076205/150‑1 (MQ=255)
                                                                                 aTACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCt                                                                      >  2:1666968/1‑150 (MQ=255)
                                                                                  tACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCtt                                                                     <  1:3582893/150‑1 (MQ=255)
                                                                                   aCGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTa                                                                    <  1:1112821/150‑1 (MQ=255)
                                                                                        tCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTaggagg                                                               <  1:3600766/150‑1 (MQ=255)
                                                                                         ccACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGgt                                                              >  2:3172736/1‑150 (MQ=255)
                                                                                          cACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGgtg                                                             >  1:2173635/1‑150 (MQ=255)
                                                                                          cACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGgtg                                                             >  1:153621/1‑150 (MQ=255)
                                                                                              aTCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACt                                                         <  1:2432071/150‑1 (MQ=255)
                                                                                               tCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACtt                                                        <  2:3434010/150‑1 (MQ=255)
                                                                                               tCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACtt                                                        <  1:1612618/150‑1 (MQ=255)
                                                                                                cTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTg                                                       <  1:2582573/150‑1 (MQ=255)
                                                                                                cTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTg                                                       >  2:400726/1‑150 (MQ=255)
                                                                                                   aTGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGttt                                                    >  2:3005422/1‑150 (MQ=255)
                                                                                                   aTGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGttt                                                    <  1:2621859/150‑1 (MQ=255)
                                                                                                     ggggCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTca                                                  <  1:2875955/150‑1 (MQ=255)
                                                                                                      gggCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTcac                                                 >  1:1073407/1‑150 (MQ=255)
                                                                                                       ggCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTcaca                                                >  1:3140635/1‑150 (MQ=255)
                                                                                                        gCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCacac                                               >  1:2201266/1‑150 (MQ=255)
                                                                                                          tctTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACgg                                             >  2:3966426/1‑150 (MQ=255)
                                                                                                           ctTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGc                                            >  2:4369215/1‑150 (MQ=255)
                                                                                                             tGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCtt                                          >  1:2462621/1‑150 (MQ=255)
                                                                                                             tGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCtt                                          <  1:3589006/150‑1 (MQ=255)
                                                                                                               cTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGc                                        >  1:60879/1‑150 (MQ=255)
                                                                                                                tATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCg                                       <  1:2723756/150‑1 (MQ=255)
                                                                                                                    cTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCt                                   >  2:4499110/1‑150 (MQ=255)
                                                                                                                       gAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATc                                <  1:494931/150‑1 (MQ=255)
                                                                                                                        aaaTGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCg                               <  2:1496251/150‑1 (MQ=255)
                                                                                                                          atgatgAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGcc                             <  1:3825520/150‑1 (MQ=255)
                                                                                                                            gatgaAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAg                           >  2:4563606/1‑150 (MQ=255)
                                                                                                                              tgaAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTa                         <  1:3009578/150‑1 (MQ=255)
                                                                                                                               gaAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTAc                        <  1:3422015/150‑1 (MQ=255)
                                                                                                                                 aGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACgg                      >  1:1304837/1‑150 (MQ=255)
                                                                                                                                     tctcAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACt                  <  1:23336/150‑1 (MQ=255)
                                                                                                                                             gtggCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGagg          <  1:1161017/150‑1 (MQ=255)
                                                                                                                                             gtggCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGagg          >  2:1010159/1‑150 (MQ=255)
                                                                                                                                              tggCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATTCAGCACTCAAGCGTGATGGCTTAGGAGGTGTACTTGTTTCACACGGCTAGCTATATATCGCCAGTACGCcacacactcagta         >  2:731360/1‑143 (MQ=255)
                                                                                                                                              tggCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGagga         >  1:1816979/1‑150 (MQ=255)
                                                                                                                                              tggCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGagga         >  1:3708039/1‑150 (MQ=255)
                                                                                                                                                   gagaATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGAt    >  2:1734388/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTT  >  CP009273/2750275‑2750573

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: