Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,251,120 T→A 24.0% T97S (ACC→TCC)  ycgV ← putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,251,1200TA24.0% 493.1 / 176.5 233T97S (ACC→TCC) ycgVputative adhesin
Reads supporting (aligned to +/- strand):  ref base T (85/92);  new base A (27/29);  total (112/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 240 total aligned reads displayed.

AGCATGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGTATCAG  >  CP009273/1250971‑1251269
                                                                                                                                                     |                                                                                                                                                     
aGCATGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGt                                                                                                                                                       >  1:3614281/1‑150 (MQ=255)
aGCATGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGt                                                                                                                                                       >  1:2848565/1‑150 (MQ=255)
   aTGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTTTGAGACGGTTAGCAGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTACTATTTCAAATCCAGAGGCGTCGTTGACATTGACGGTgtc                                                                                                                                                    >  1:144909/1‑148 (MQ=255)
   aTGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGAc                                                                                                                                                    >  1:2321787/1‑150 (MQ=255)
    tGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGAca                                                                                                                                                   <  2:2861775/150‑1 (MQ=255)
      ggTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGAcatc                                                                                                                                                 >  2:1375651/1‑150 (MQ=255)
         ataGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCAtt                                                                                                                                              <  1:1864756/150‑1 (MQ=255)
          taGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGAGTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTa                                                                                                                                             <  2:2699931/150‑1 (MQ=255)
           aGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTAc                                                                                                                                            >  1:600851/1‑150 (MQ=255)
           aGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTAc                                                                                                                                            >  1:600969/1‑150 (MQ=255)
            gTCACCAATTAAGTTAATGCCGATGGCAGATGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTACTTTTGAATATCCAGAGGCGTCGTTGACATAGAAGGTGACAACATTAcc                                                                                                                                           >  2:383716/1‑150 (MQ=255)
               aCCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAg                                                                                                                                        <  2:24518/150‑1 (MQ=255)
                ccAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAgg                                                                                                                                       >  2:2514208/1‑150 (MQ=255)
                    ttAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGa                                                                                                                                   >  2:3323781/1‑150 (MQ=255)
                    ttAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGa                                                                                                                                   >  2:238857/1‑150 (MQ=255)
                     tAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGat                                                                                                                                  >  2:1688928/1‑150 (MQ=255)
                        gTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATAtt                                                                                                                               <  2:2157292/150‑1 (MQ=255)
                            aTGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCa                                                                                                                           <  2:3194075/150‑1 (MQ=255)
                            aTGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCa                                                                                                                           <  2:3980710/150‑1 (MQ=255)
                              gCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAgg                                                                                                                         <  2:3911044/150‑1 (MQ=255)
                                cGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCg                                                                                                                       >  2:387215/1‑150 (MQ=255)
                                 gATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCgt                                                                                                                      >  1:1215243/1‑150 (MQ=255)
                                   tGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCgtgt                                                                                                                    >  2:53169/1‑150 (MQ=255)
                                   tGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCgtgt                                                                                                                    <  1:3726242/150‑1 (MQ=255)
                                    ggCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGtt                                                                                                                   <  2:3561033/150‑1 (MQ=255)
                                      cAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTc                                                                                                                 >  1:2512085/1‑150 (MQ=255)
                                      cAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTc                                                                                                                 >  1:3963370/1‑150 (MQ=255)
                                      cAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATACAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCATGCGTGTTTc                                                                                                                 >  1:3395017/1‑150 (MQ=255)
                                       agagGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTcc                                                                                                                >  2:1122663/1‑150 (MQ=255)
                                        gagGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCct                                                                                                               <  2:2986618/150‑1 (MQ=255)
                                        gagGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCct                                                                                                               <  1:2260370/150‑1 (MQ=255)
                                          ggTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCctgc                                                                                                             <  2:2249126/150‑1 (MQ=255)
                                           gTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCctgct                                                                                                            >  1:715449/1‑150 (MQ=255)
                                           gTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCctgct                                                                                                            >  1:4132191/1‑150 (MQ=255)
                                             ttGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTgg                                                                                                          <  2:1735005/150‑1 (MQ=255)
                                              tGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTggg                                                                                                         <  2:118418/150‑1 (MQ=255)
                                               gACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTgggg                                                                                                        <  2:641323/150‑1 (MQ=255)
                                                  caacaaCATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCa                                                                                                     >  2:3232034/1‑150 (MQ=255)
                                                     caacaTCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTTCCAGGCGTGTTTCCTGCTGGGGTCATTa                                                                                                  >  1:219948/1‑150 (MQ=255)
                                                     caacaTCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGGTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTa                                                                                                  >  2:2818834/1‑150 (MQ=255)
                                                      aacaTCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTACTGCTGGGGTCATTaa                                                                                                 >  1:554870/1‑149 (MQ=255)
                                                        caTCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCAGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGAAGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTacac                                                                                               >  1:245851/1‑150 (MQ=255)
                                                         aTCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGGTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACAcc                                                                                              >  1:1963733/1‑150 (MQ=255)
                                                          tCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCa                                                                                             >  2:1623111/1‑150 (MQ=255)
                                                           cTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCAt                                                                                            >  2:2156812/1‑150 (MQ=255)
                                                            tACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATa                                                                                           <  2:2273606/150‑1 (MQ=255)
                                                               tgtgAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATagag                                                                                        <  1:1140558/150‑1 (MQ=255)
                                                                gtgAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAgg                                                                                       <  2:716312/150‑1 (MQ=255)
                                                                     tCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTcc                                                                                  <  1:427313/150‑1 (MQ=255)
                                                                         ttAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGAt                                                                              >  1:2893007/1‑150 (MQ=255)
                                                                          tAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATg                                                                             >  2:2531698/1‑150 (MQ=255)
                                                                             cTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGAt                                                                          <  1:3917401/150‑1 (MQ=255)
                                                                                  aGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATcggcg                                                                     <  2:1215243/150‑1 (MQ=255)
                                                                                  aGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATcggcg                                                                     <  2:4359933/150‑1 (MQ=255)
                                                                                    ggAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGt                                                                   <  1:899007/150‑1 (MQ=255)
                                                                                     gAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTa                                                                  <  2:1037486/150‑1 (MQ=255)
                                                                                     gAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTa                                                                  >  2:694524/1‑150 (MQ=255)
                                                                                       gCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTAtc                                                                <  2:3029959/150‑1 (MQ=255)
                                                                                          gTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTcc                                                             <  2:4192506/150‑1 (MQ=255)
                                                                                           tttttGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCt                                                            <  2:1521208/150‑1 (MQ=255)
                                                                                            ttttGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTg                                                           <  2:1820539/150‑1 (MQ=255)
                                                                                               tGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTaa                                                        >  2:2022382/1‑150 (MQ=255)
                                                                                                 cccTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAAtt                                                      >  2:3373376/1‑150 (MQ=255)
                                                                                                   cTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCt                                                    <  1:3512683/150‑1 (MQ=255)
                                                                                                   cTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCt                                                    >  2:2572768/1‑150 (MQ=255)
                                                                                                    tGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTg                                                   <  2:3606918/150‑1 (MQ=255)
                                                                                                    tGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTg                                                   >  2:1123657/1‑150 (MQ=255)
                                                                                                      aataatGTTTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTc                                                 <  2:3591031/150‑1 (MQ=255)
                                                                                                       ataatGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCa                                                >  1:319281/1‑150 (MQ=255)
                                                                                                       ataatGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCa                                                <  2:740668/150‑1 (MQ=255)
                                                                                                       ataatGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCa                                                <  1:811588/150‑1 (MQ=255)
                                                                                                           tGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACcat                                            <  2:1881061/150‑1 (MQ=255)
                                                                                                           tGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACcat                                            >  1:4224453/1‑150 (MQ=255)
                                                                                                           tGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACcat                                            <  1:2631712/150‑1 (MQ=255)
                                                                                                            gATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGTTATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACcatc                                           >  1:210179/1‑150 (MQ=255)
                                                                                                                 cTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTaa                                      <  1:1341124/150‑1 (MQ=255)
                                                                                                                 cTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTACTGATGGATCGGCGGTATATCCTGTAATTCTGTAACCATCATTaa                                      >  1:2960158/1‑150 (MQ=255)
                                                                                                                    ttGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAg                                   <  1:1688928/150‑1 (MQ=255)
                                                                                                                     tGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGt                                  <  2:1034823/150‑1 (MQ=255)
                                                                                                                      gCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTa                                 >  2:3131160/1‑150 (MQ=255)
                                                                                                                           tCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATAtt                            >  1:1757860/1‑150 (MQ=255)
                                                                                                                            ccAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTg                           <  2:3312490/150‑1 (MQ=255)
                                                                                                                            ccAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTg                           >  2:1266375/1‑150 (MQ=255)
                                                                                                                            ccAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTg                           <  1:1015865/150‑1 (MQ=255)
                                                                                                                             cAGAGGCGTCGTTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTgg                          <  2:2495893/150‑1 (MQ=255)
                                                                                                                              agagGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGt                         >  1:357824/1‑150 (MQ=255)
                                                                                                                                    gtcgtTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCata                   >  1:2734960/1‑150 (MQ=255)
                                                                                                                                    gtcgtTGACATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCata                   <  2:1664534/150‑1 (MQ=255)
                                                                                                                                      cgtTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATAtt                 <  2:903718/150‑1 (MQ=255)
                                                                                                                                      cgtTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATAtt                 <  1:3282777/150‑1 (MQ=255)
                                                                                                                                      cgtTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATAtt                 >  2:2681573/1‑150 (MQ=255)
                                                                                                                                          gACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACt             >  2:1913830/1‑150 (MQ=255)
                                                                                                                                          gACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACt             >  2:1420222/1‑150 (MQ=255)
                                                                                                                                           aCATTGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACtt            >  1:2082585/1‑150 (MQ=255)
                                                                                                                                               tGACGGAGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTg        <  2:4179718/150‑1 (MQ=255)
                                                                                                                                                gACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGt       >  1:1632011/1‑150 (MQ=255)
                                                                                                                                                gACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGt       >  1:3819986/1‑150 (MQ=255)
                                                                                                                                                    gTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGTATca   <  1:3682512/150‑1 (MQ=255)
                                                                                                                                                    gTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGTATca   >  1:3999046/1‑150 (MQ=255)
                                                                                                                                                     tGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGTATcag  >  2:938756/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCATGGGTATAGTCACCAATTAAGTTAATGCCGATGGCAGAGGTTTGACCAACAACATCTACTGTGAGTCGGTTAGCTGTCAGGGAGCTGTTTTTGCCCTGAATAATGATTCCTTTTGCATATCCAGAGGCGTCGTTGACATTGACGGTGACATCATTACCCAGGTTGATATTGCCAGGCGTGTTTCCTGCTGGGGTCATTACACCATAGAGGTTTCCTGATGGATCGGCGGTATCTCCTGTAATTCTGTCACCATCATTAATAGTAATATTGGTGCCATATTGACTTCCTGTATCAG  >  CP009273/1250971‑1251269

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: ‑