Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,036,569 G→T 9.8% L265L (CTG→CTT srkA → stress response kinase A; MazF antagonist; Thr/Ser kinase involved in Cpx stress response

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,036,5690GT 9.8% 1009.5 / 54.3 284L265L (CTG→CTTsrkAstress response kinase A; MazF antagonist; Thr/Ser kinase involved in Cpx stress response
Reads supporting (aligned to +/- strand):  ref base G (129/127);  new base T (12/16);  total (141/143)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 284 total aligned reads displayed.

GATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAACAGGCA  >  CP009273/4036423‑4036716
                                                                                                                                                  |                                                                                                                                                   
gATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGAt                                                                                                                                                   <  1:4133799/149‑1 (MQ=255)
  tCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTg                                                                                                                                                 <  2:3079469/149‑1 (MQ=255)
   cTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGa                                                                                                                                                <  1:752673/149‑1 (MQ=255)
    tGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGaa                                                                                                                                               >  2:4081309/1‑149 (MQ=255)
     ggATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAAc                                                                                                                                              <  2:3882751/149‑1 (MQ=255)
      gatgatGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAAcc                                                                                                                                             >  2:2047337/1‑149 (MQ=255)
        tgatgCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAACCtt                                                                                                                                           >  2:2198818/1‑149 (MQ=255)
          atgCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTa                                                                                                                                         >  2:1845667/1‑149 (MQ=255)
              aCGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTAcgcg                                                                                                                                     >  2:4160530/1‑149 (MQ=255)
                 tAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCa                                                                                                                                  >  2:2700643/1‑149 (MQ=255)
                 tAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCa                                                                                                                                  >  2:1036680/1‑149 (MQ=255)
                   aTGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATg                                                                                                                                >  2:522648/1‑149 (MQ=255)
                   aTGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACGATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATg                                                                                                                                >  2:2484108/1‑149 (MQ=255)
                        ccAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCTCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGttt                                                                                                                           <  2:122449/149‑1 (MQ=255)
                          aGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTgg                                                                                                                         >  2:2974637/1‑149 (MQ=255)
                              gTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGTAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGGTTGGTtta                                                                                                                     >  1:3150422/1‑149 (MQ=255)
                                   ggATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATc                                                                                                                <  2:4694182/149‑1 (MQ=255)
                                     aTTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCtt                                                                                                              <  1:1635408/149‑1 (MQ=255)
                                     aTTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCtt                                                                                                              <  2:1685884/149‑1 (MQ=255)
                                     aTTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCtt                                                                                                              >  1:3341086/1‑149 (MQ=255)
                                     aTTAGTGGATGTTGCTCAATGGCGAAAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCtt                                                                                                              <  2:1992246/149‑1 (MQ=255)
                                       ttGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGc                                                                                                            <  1:298672/149‑1 (MQ=255)
                                        tgtgGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGcc                                                                                                           <  2:2300858/149‑1 (MQ=255)
                                        tgtgGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGcc                                                                                                           <  1:3199645/149‑1 (MQ=255)
                                         gtgGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCt                                                                                                          >  1:2356386/1‑149 (MQ=255)
                                         gtgGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCt                                                                                                          >  1:3298844/1‑149 (MQ=255)
                                         gtgGATGGTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCt                                                                                                          <  1:4262179/149‑1 (MQ=255)
                                          tgGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAACCTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTg                                                                                                         >  2:180344/1‑149 (MQ=255)
                                          tgGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTg                                                                                                         >  2:1428999/1‑149 (MQ=255)
                                           ggATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTgg                                                                                                        >  2:1002324/1‑149 (MQ=255)
                                           ggATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCACTGCGTGTGGTTTATTATCTTGCCTgg                                                                                                        >  1:2425438/1‑149 (MQ=255)
                                                ttGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAAt                                                                                                   <  2:1067008/149‑1 (MQ=255)
                                                 tGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATg                                                                                                  <  2:2655681/149‑1 (MQ=255)
                                                 tGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATg                                                                                                  <  2:2655354/149‑1 (MQ=255)
                                                  gCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTGTATTATCTTGCCTGGCTAATgc                                                                                                 >  1:3501910/1‑149 (MQ=255)
                                                   cTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATgcg                                                                                                >  1:818535/1‑149 (MQ=255)
                                                    tCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATgcgg                                                                                               >  2:3529808/1‑149 (MQ=255)
                                                    tCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATgcgg                                                                                               <  1:3071033/149‑1 (MQ=255)
                                                      aaTGGCGATACAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTGAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATgcggcg                                                                                             <  2:2594246/149‑1 (MQ=255)
                                                          gCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTgg                                                                                         >  2:2014119/1‑149 (MQ=255)
                                                            gATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCTATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGc                                                                                       >  1:4107719/1‑149 (MQ=255)
                                                            gATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGc                                                                                       <  2:1667750/149‑1 (MQ=255)
                                                            gATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGc                                                                                       <  2:1145819/149‑1 (MQ=255)
                                                               aaaGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGa                                                                                    >  1:4623486/1‑149 (MQ=255)
                                                               aaaGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGa                                                                                    <  1:1285607/149‑1 (MQ=255)
                                                                aaGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGAt                                                                                   >  1:1167321/1‑149 (MQ=255)
                                                                 aGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATc                                                                                  <  2:2835074/149‑1 (MQ=255)
                                                                   ccGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATccc                                                                                <  2:4586147/149‑1 (MQ=255)
                                                                     gAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCcgc                                                                              >  1:2087767/1‑149 (MQ=255)
                                                                      agcagcGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCcgcg                                                                             >  1:3398547/1‑149 (MQ=255)
                                                                             gATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGa                                                                      >  1:386103/1‑149 (MQ=255)
                                                                             gATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGa                                                                      <  2:683713/149‑1 (MQ=255)
                                                                              aTGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGaa                                                                     >  2:4012191/1‑149 (MQ=255)
                                                                                gCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGaaaa                                                                   >  1:2260495/1‑149 (MQ=255)
                                                                                   aCTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAAtt                                                                >  1:1894834/1‑149 (MQ=255)
                                                                                    cTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAAttt                                                               <  1:4128631/149‑1 (MQ=255)
                                                                                     tGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTc                                                              <  1:606670/149‑1 (MQ=255)
                                                                                      ggAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTcc                                                             <  1:2556408/149‑1 (MQ=255)
                                                                                       gAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTccc                                                            >  2:3660908/1‑149 (MQ=255)
                                                                                       gAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTccc                                                            <  2:4371199/149‑1 (MQ=255)
                                                                                       gAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTccc                                                            >  2:2931507/1‑149 (MQ=255)
                                                                                         aaCTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGt                                                          >  2:3349116/1‑149 (MQ=255)
                                                                                         aaCTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGt                                                          <  1:4202937/149‑1 (MQ=255)
                                                                                         aaCTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGt                                                          >  2:4133155/1‑149 (MQ=255)
                                                                                          aCTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTg                                                         <  1:4744625/149‑1 (MQ=255)
                                                                                          aCTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTg                                                         >  2:4238654/1‑149 (MQ=255)
                                                                                           cTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTgg                                                        <  1:2906676/149‑1 (MQ=255)
                                                                                              ttattGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTa                                                     <  2:884490/149‑1 (MQ=255)
                                                                                                 ttGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAAcc                                                  >  1:21878/1‑149 (MQ=255)
                                                                                                   gAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCgg                                                >  2:1453044/1‑149 (MQ=255)
                                                                                                    aaGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCggg                                               <  1:893412/149‑1 (MQ=255)
                                                                                                    aaGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCggg                                               <  1:1755936/149‑1 (MQ=255)
                                                                                                     aGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCgggg                                              >  2:3579130/1‑149 (MQ=255)
                                                                                                      gCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTCTTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGa                                             >  2:1354576/1‑149 (MQ=255)
                                                                                                      gCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGa                                             >  1:1541875/1‑149 (MQ=255)
                                                                                                          aTGAAGAGTTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGAt                                         >  2:2532371/1‑149 (MQ=255)
                                                                                                           tGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGAtt                                        <  1:454515/149‑1 (MQ=255)
                                                                                                             aagaaTTTAGCGAGTTCGACACCGCTGAAATCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTAc                                      <  2:2390833/149‑1 (MQ=255)
                                                                                                                 atgtAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGc                                  <  2:4452922/146‑1 (MQ=255)
                                                                                                                    tAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGc                               >  2:937661/1‑149 (MQ=255)
                                                                                                                     aGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCg                              >  1:368314/1‑149 (MQ=255)
                                                                                                                        gAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACa                           >  1:1752082/1‑149 (MQ=255)
                                                                                                                          gTTCGACACCGCTGAACTCGGACTGATGGAACCTGTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGa                         <  2:4645344/149‑1 (MQ=255)
                                                                                                                          gTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGa                         <  2:1644030/149‑1 (MQ=255)
                                                                                                                            tCGACACCGCTGAAATCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACg                       >  2:3804734/1‑149 (MQ=255)
                                                                                                                             cGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACgg                      <  2:2966059/149‑1 (MQ=255)
                                                                                                                               acacCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCg                    >  1:3795638/1‑149 (MQ=255)
                                                                                                                                 acCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGAc                  >  2:2359905/1‑149 (MQ=255)
                                                                                                                                 acCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGAc                  <  1:2621748/149‑1 (MQ=255)
                                                                                                                                 acCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGAc                  <  1:1042131/149‑1 (MQ=255)
                                                                                                                                  ccGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACt                 >  2:2867327/1‑149 (MQ=255)
                                                                                                                                   cGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACtt                <  1:2969659/149‑1 (MQ=255)
                                                                                                                                       gAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTtat            <  2:2788315/149‑1 (MQ=255)
                                                                                                                                       gAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGAGTACTGGCTGCGACAGACGGCGACTTTtat            >  2:3412219/1‑149 (MQ=255)
                                                                                                                                          aTCGGACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGa         >  1:4358043/1‑149 (MQ=255)
                                                                                                                                            cGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAAc       >  1:4724736/1‑149 (MQ=255)
                                                                                                                                            cGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAAc       <  1:2272085/149‑1 (MQ=255)
                                                                                                                                             ggACTTATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAACa      <  1:1965454/149‑1 (MQ=255)
                                                                                                                                             ggACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAACa      <  1:1002324/149‑1 (MQ=255)
                                                                                                                                                cTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAACAGGc   <  2:4290689/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGAATTTAGCGAGTTCGACACCGCTGAAATCGGACTGATTGAACCTTTACGCGCCATGCGTTTGGTTTATTATCTTGCCTGGCTAATGCGGCGTTGGGCTGATCCCGCGTTCCCGAAAAATTTCCCGTGGTTAACCGGGGAAGATTACTGGCTGCGACAGACGGCGACTTTTATAGAACAGGCA  >  CP009273/4036423‑4036716

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: