Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,500,138 G→T 37.7% I186I (ATC→ATA yhfT ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,500,1380GT37.7% 391.0 / 280.6 277I186I (ATC→ATAyhfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (81/91);  new base T (45/60);  total (126/151)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 278 total aligned reads displayed.

CACGATGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCTGAAACCA  >  CP009273/3499991‑3500286
                                                                                                                                                   |                                                                                                                                                    
cACGATGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGAt                                                                                                                                                    <  2:132996/150‑1 (MQ=255)
  cGATGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGAttt                                                                                                                                                  <  2:1787714/150‑1 (MQ=255)
   gATGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTc                                                                                                                                                 <  1:1803460/150‑1 (MQ=255)
    aTGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCg                                                                                                                                                <  2:3188946/150‑1 (MQ=255)
    aTGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCg                                                                                                                                                >  1:2953208/1‑150 (MQ=255)
    aTGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCg                                                                                                                                                <  1:4580102/150‑1 (MQ=255)
    aTGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCg                                                                                                                                                <  1:2265475/150‑1 (MQ=255)
    aTGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCg                                                                                                                                                >  1:2075346/1‑150 (MQ=255)
    aTGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCg                                                                                                                                                >  1:1935333/1‑150 (MQ=255)
       gCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGACAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTTTCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATg                                                                                                                                             <  1:3969929/150‑1 (MQ=255)
        cGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATgg                                                                                                                                            <  1:2127296/150‑1 (MQ=255)
        cGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATgg                                                                                                                                            <  1:2145306/150‑1 (MQ=255)
          atatAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGAt                                                                                                                                          <  2:3084693/150‑1 (MQ=255)
           tataGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGAtt                                                                                                                                         <  2:4882/150‑1 (MQ=255)
           tataGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGAtt                                                                                                                                         <  2:5017/150‑1 (MQ=255)
            ataGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTc                                                                                                                                        >  1:1527954/1‑150 (MQ=255)
            ataGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTc                                                                                                                                        >  2:4213417/1‑150 (MQ=255)
            ataGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTc                                                                                                                                        <  2:1305968/150‑1 (MQ=255)
             taGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCa                                                                                                                                       <  2:115198/150‑1 (MQ=255)
              aGGGTAAGTTTTTGATAATCCGTGAGGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAATTATTTCGATGGATTCAg                                                                                                                                      <  2:367200/150‑1 (MQ=255)
               gggTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAgg                                                                                                                                     >  1:2496089/1‑150 (MQ=255)
                     gTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAg                                                                                                                               <  1:4476736/150‑1 (MQ=255)
                     gTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAg                                                                                                                               <  1:2781746/150‑1 (MQ=255)
                     gTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAg                                                                                                                               <  2:2163314/150‑1 (MQ=255)
                     gTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAg                                                                                                                               >  1:3642668/1‑150 (MQ=255)
                      tttttGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGa                                                                                                                              <  2:2901866/150‑1 (MQ=255)
                      tttttGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGa                                                                                                                              <  2:3230291/150‑1 (MQ=255)
                       ttttGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGAt                                                                                                                             <  1:1071230/150‑1 (MQ=255)
                          tGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGAtgtg                                                                                                                          >  1:2426695/1‑150 (MQ=255)
                          tGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGAtgtg                                                                                                                          >  1:1961495/1‑150 (MQ=255)
                          tGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGAGTTCGATGGATTCAGGGTTAAGAtgtg                                                                                                                          >  2:2071003/1‑150 (MQ=255)
                           gATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTgg                                                                                                                         >  1:3128158/1‑150 (MQ=255)
                           gATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTgg                                                                                                                         <  2:3011261/150‑1 (MQ=255)
                            aTAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTggg                                                                                                                        >  2:788931/1‑150 (MQ=255)
                               aTCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGaaa                                                                                                                     >  1:1549985/1‑150 (MQ=255)
                                 ccGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATa                                                                                                                   <  1:164971/150‑1 (MQ=255)
                                   gTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAgc                                                                                                                 >  1:4075957/1‑150 (MQ=255)
                                    tGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAgcg                                                                                                                <  1:4638651/150‑1 (MQ=255)
                                    tGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAgcg                                                                                                                <  1:90779/150‑1 (MQ=255)
                                     gACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAgcgc                                                                                                               >  1:4196904/1‑150 (MQ=255)
                                      aCGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCa                                                                                                              >  2:4148843/1‑150 (MQ=255)
                                       cGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCAc                                                                                                             <  2:2120292/150‑1 (MQ=255)
                                       cGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCAc                                                                                                             <  2:664940/150‑1 (MQ=255)
                                            gTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTa                                                                                                        <  2:4260524/150‑1 (MQ=255)
                                            gTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTa                                                                                                        >  2:1651576/1‑150 (MQ=255)
                                            gTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTa                                                                                                        >  2:1160864/1‑150 (MQ=255)
                                               cttcAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTTATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACca                                                                                                     <  1:505791/150‑1 (MQ=255)
                                               cttcAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACca                                                                                                     <  2:2037404/150‑1 (MQ=255)
                                                  cAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCacac                                                                                                  >  1:3281906/1‑150 (MQ=255)
                                                   aaaCACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACg                                                                                                 <  1:3027056/150‑1 (MQ=255)
                                                   aaaCACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACg                                                                                                 <  2:2400461/150‑1 (MQ=255)
                                                   aaaCACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACg                                                                                                 >  1:403375/1‑150 (MQ=255)
                                                   aaaCACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACg                                                                                                 <  2:87169/150‑1 (MQ=255)
                                                     acacCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACggg                                                                                               >  1:3714107/1‑150 (MQ=255)
                                                       acCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGtc                                                                                             >  1:3499970/1‑150 (MQ=255)
                                                        ccGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGtca                                                                                            <  2:2729035/150‑1 (MQ=255)
                                                         cGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGtcat                                                                                           <  2:1197601/150‑1 (MQ=255)
                                                          gAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGtcatc                                                                                          >  2:4417575/1‑150 (MQ=255)
                                                            aaGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATAACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAg                                                                                        <  2:1105026/150‑1 (MQ=255)
                                                             aGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGt                                                                                       <  2:3924089/150‑1 (MQ=255)
                                                               cccGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGACGTGAGATATCGCGATCCCCAGCAGCATCACCATGCCAAGAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTac                                                                                     >  2:355591/1‑150 (MQ=255)
                                                                     ggCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACcaccac                                                                               >  2:1951617/1‑150 (MQ=255)
                                                                     ggCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACcaccac                                                                               >  2:2794726/1‑150 (MQ=255)
                                                                      gCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCAcg                                                                              <  2:1861761/150‑1 (MQ=255)
                                                                        aTCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCAcggc                                                                            <  1:77661/150‑1 (MQ=255)
                                                                        aTCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCAcggc                                                                            <  1:4321532/150‑1 (MQ=255)
                                                                        aTCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCAcggc                                                                            >  1:872827/1‑150 (MQ=255)
                                                                         tCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATACATATTTCGTTGGTTTCAGGGTTATGATGTGGGAATTTGCGCACGACTACCACACGGGTCATCAGTACCACCAcggcg                                                                           <  1:2694554/150‑1 (MQ=255)
                                                                             tGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGa                                                                       >  1:621288/1‑150 (MQ=255)
                                                                                cATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCa                                                                    <  2:2524806/150‑1 (MQ=255)
                                                                                   tttCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGAc                                                                 >  2:1697824/1‑150 (MQ=255)
                                                                                   tttCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGAc                                                                 >  2:3479350/1‑150 (MQ=255)
                                                                                    ttCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACt                                                                <  2:808983/150‑1 (MQ=255)
                                                                                     tcatcaCGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACtt                                                               <  1:2034934/150‑1 (MQ=255)
                                                                                      catcaCGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACttt                                                              <  1:2819579/150‑1 (MQ=255)
                                                                                       atcaCGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTg                                                             >  2:4257719/1‑150 (MQ=255)
                                                                                           cGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTc                                                         >  1:1221164/1‑150 (MQ=255)
                                                                                                aCGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGcc                                                    >  1:3821901/1‑150 (MQ=255)
                                                                                                    aGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAAc                                                >  2:3381922/1‑150 (MQ=255)
                                                                                                        cGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGt                                            <  1:699975/150‑1 (MQ=255)
                                                                                                         gTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTa                                           >  1:1382176/1‑150 (MQ=255)
                                                                                                             gTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGaaa                                       >  2:2854536/1‑150 (MQ=255)
                                                                                                                 tCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCg                                   <  1:2580995/150‑1 (MQ=255)
                                                                                                                    cgATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGcca                                >  1:2724458/1‑150 (MQ=255)
                                                                                                                     gATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGccac                               >  1:2668710/1‑150 (MQ=255)
                                                                                                                      aTCCCCAGCAGCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGccacc                              >  1:1168681/1‑150 (MQ=255)
                                                                                                                       taaccAGCAGCATCACCATGCCAATAAAGATTTCGAAGGATTCAGGGTTAAGATTTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGccacca                             >  1:529977/4‑150 (MQ=255)
                                                                                                                           cagcagCATCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCgg                         >  2:4087560/1‑150 (MQ=255)
                                                                                                                               agcaTCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAAc                     >  2:1534139/1‑150 (MQ=255)
                                                                                                                                gcaTCACCATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACa                    <  1:1575654/150‑1 (MQ=255)
                                                                                                                                  aTCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGt                  <  2:692439/150‑1 (MQ=255)
                                                                                                                                       cATGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCaaa             >  2:1851290/1‑150 (MQ=255)
                                                                                                                                       cATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCaaa             <  1:788931/150‑1 (MQ=255)
                                                                                                                                        aTGCCAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCaaaa            <  2:3732197/150‑1 (MQ=255)
                                                                                                                                           ccAATAAATATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCt         <  1:3160973/150‑1 (MQ=255)
                                                                                                                                           ccAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGTAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCt         >  1:3071476/1‑150 (MQ=255)
                                                                                                                                           ccAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCt         <  2:3592694/150‑1 (MQ=255)
                                                                                                                                             aaTAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCTGa       <  2:2871372/150‑1 (MQ=255)
                                                                                                                                                aaaGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCTGAAac    <  2:4471903/150‑1 (MQ=255)
                                                                                                                                                  aGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCTGAAacca  >  2:3628386/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACGATGGCGATATAGGGTAAGTTTTTGATAATCCGTGACGTGCGTTCTTCAAACACCGAAAGCCCGCTGGCATCAATGTCATTTTCATCACGATGACGCAGGTCGTGAGTTATCGCGATCCCCAGCAGCATCACCATGCCAATAAAGATTTCGATGGATTCAGGGTTAAGATGTGGGAAATAGCGCACGACTACCACACGGGTCATCAGTACCACCACGGCGGCGATCAGACTTTGCTTCCAGCCAAACTGGTAGAAAATCGCCACCAGCGGGAACAGTGCAAAAGCTGAAACCA  >  CP009273/3499991‑3500286

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: