Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,078,761 C→A 10.6% R141L (CGC→CTC)  fdoG ← formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,078,7610CA10.6% 1004.4 / 56.2 298R141L (CGC→CTC) fdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base C (135/130);  new base A (17/15);  total (152/146)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GTCGACCGCGAGCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGAC  >  CP009273/4078613‑4078908
                                                                                                                                                    |                                                                                                                                                   
   gACCGCGAGCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGt                                                                                                                                                  <  2:3330937/149‑1 (MQ=255)
    aCCGCGAGCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGtc                                                                                                                                                 <  2:4574433/149‑1 (MQ=255)
          aGCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCttt                                                                                                                                           <  1:4153497/149‑1 (MQ=255)
           gCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTtc                                                                                                                                          <  2:3452485/149‑1 (MQ=255)
           gCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTtc                                                                                                                                          <  1:3048126/149‑1 (MQ=255)
             aTACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTtcat                                                                                                                                        <  2:2334170/149‑1 (MQ=255)
              tACCCGGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTtcatc                                                                                                                                       <  1:3308322/149‑1 (MQ=255)
              tACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTtcatc                                                                                                                                       <  1:2455922/149‑1 (MQ=255)
               aCCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTtcatca                                                                                                                                      <  2:4077572/149‑1 (MQ=255)
                 ccAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGt                                                                                                                                    <  1:3148237/149‑1 (MQ=255)
                  cAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGtt                                                                                                                                   <  2:725953/149‑1 (MQ=255)
                   aGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCTTCAGttt                                                                                                                                  <  1:3760616/149‑1 (MQ=255)
                   aGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGttt                                                                                                                                  <  1:2928568/149‑1 (MQ=255)
                   aGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGttt                                                                                                                                  <  1:4656610/149‑1 (MQ=255)
                    gcgcgcGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTg                                                                                                                                 >  1:2811646/1‑149 (MQ=255)
                        gcgGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcga                                                                                                                             >  1:855996/1‑149 (MQ=255)
                         cgGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcgat                                                                                                                            <  2:1679961/149‑1 (MQ=255)
                          gggAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcgatg                                                                                                                           >  1:4496400/1‑149 (MQ=255)
                          gggAGAATTATTTCTTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATACCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTATCATCAGTATGgagatg                                                                                                                           >  2:175371/1‑149 (MQ=255)
                           ggAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcgatgc                                                                                                                          <  1:899847/149‑1 (MQ=255)
                           ggAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGgcgatgc                                                                                                                          <  1:943258/149‑1 (MQ=255)
                            gagaATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcgatgcg                                                                                                                         >  1:1569555/1‑149 (MQ=255)
                            gagaATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcgatgcg                                                                                                                         >  2:3891402/1‑149 (MQ=255)
                            gagaATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGgcgatgcg                                                                                                                         >  1:2297373/1‑149 (MQ=255)
                                  ttttGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCaaa                                                                                                                   <  2:1074408/149‑1 (MQ=255)
                                  ttttGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCaaa                                                                                                                   <  1:219602/149‑1 (MQ=255)
                                  ttttGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCaaa                                                                                                                   <  2:3050228/149‑1 (MQ=255)
                                  ttttGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCaaa                                                                                                                   <  2:3877516/149‑1 (MQ=255)
                                        gTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCctc                                                                                                             >  1:4163770/1‑149 (MQ=255)
                                        gTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCctc                                                                                                             >  2:2389481/1‑149 (MQ=255)
                                         ttAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCctct                                                                                                            >  2:2196828/1‑149 (MQ=255)
                                         ttAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCctct                                                                                                            >  1:3591659/1‑149 (MQ=255)
                                         ttAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCctct                                                                                                            <  2:2911870/149‑1 (MQ=255)
                                          tAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAATGCCTCtt                                                                                                           <  2:4737748/149‑1 (MQ=255)
                                            aaTAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTcc                                                                                                         >  1:2580512/1‑149 (MQ=255)
                                            aaTAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTcc                                                                                                         >  2:1159086/1‑149 (MQ=255)
                                            aaTAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTcc                                                                                                         >  1:2633313/1‑149 (MQ=255)
                                              tAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCa                                                                                                       >  2:1289442/1‑149 (MQ=255)
                                                gCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAAc                                                                                                     <  2:4263156/149‑1 (MQ=255)
                                                  cGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGGGGTTAACAGTCACGCCTTGGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTg                                                                                                   <  1:2131895/149‑1 (MQ=255)
                                                   ggTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGa                                                                                                  <  2:4778240/149‑1 (MQ=255)
                                                          ttACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTtgctg                                                                                           >  1:4366531/1‑149 (MQ=255)
                                                            aCTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGcc                                                                                         >  2:2736453/1‑149 (MQ=255)
                                                             cTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCa                                                                                        <  2:4499159/149‑1 (MQ=255)
                                                              tCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCTATCAAACGCCTCTTCCCAACTGATTTGCTGCCAt                                                                                       <  1:4210530/149‑1 (MQ=255)
                                                               cGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCAtt                                                                                      >  1:1861150/1‑149 (MQ=255)
                                                                 aCGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTa                                                                                    >  1:2071303/1‑149 (MQ=255)
                                                                    cgGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCa                                                                                 >  1:670351/1‑149 (MQ=255)
                                                                     ggAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAg                                                                                >  2:4526406/1‑149 (MQ=255)
                                                                      gAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGa                                                                               >  1:3982571/1‑149 (MQ=255)
                                                                      gAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGa                                                                               >  2:4470097/1‑149 (MQ=255)
                                                                        aGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAAc                                                                             <  2:3566324/149‑1 (MQ=255)
                                                                         gCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAAcc                                                                            <  1:1012520/149‑1 (MQ=255)
                                                                         gCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAAcc                                                                            <  2:2506128/149‑1 (MQ=255)
                                                                          cacacaGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCt                                                                           <  1:4591042/149‑1 (MQ=255)
                                                                          cacacaGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTATTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCt                                                                           <  1:1750992/149‑1 (MQ=255)
                                                                           acacaGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTg                                                                          <  2:2154269/149‑1 (MQ=255)
                                                                           acacaGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGAAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTg                                                                          <  1:3249292/149‑1 (MQ=255)
                                                                            cacaGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTgg                                                                         >  2:2101946/1‑149 (MQ=255)
                                                                                   tCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACgg                                                                  <  1:75433/149‑1 (MQ=255)
                                                                                   tCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACgg                                                                  >  1:4063993/1‑149 (MQ=255)
                                                                                        gTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTc                                                             <  2:2848114/149‑1 (MQ=255)
                                                                                         tGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTcc                                                            <  2:1483732/149‑1 (MQ=255)
                                                                                          ggAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCg                                                           >  2:3881988/1‑149 (MQ=255)
                                                                                           gagaGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCgg                                                          >  2:1021128/1‑149 (MQ=255)
                                                                                           gagaGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATAAGAACCTGGCGCACGGTATTCCgg                                                          >  1:3620942/1‑149 (MQ=255)
                                                                                            agagCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGa                                                         >  2:1603194/1‑149 (MQ=255)
                                                                                             gagCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGaa                                                        >  2:1538146/1‑149 (MQ=255)
                                                                                              agCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGaaa                                                       <  1:1197870/149‑1 (MQ=255)
                                                                                              agCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGaaa                                                       <  1:354097/149‑1 (MQ=255)
                                                                                               gCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAAc                                                      <  1:2795881/149‑1 (MQ=255)
                                                                                                ccAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACt                                                     >  1:2532018/1‑149 (MQ=255)
                                                                                                   gCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCa                                                  >  2:3645983/1‑149 (MQ=255)
                                                                                                     ggTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGa                                                <  2:1106645/149‑1 (MQ=255)
                                                                                                      gTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGACACTTCAGAc                                               >  2:3611399/1‑149 (MQ=255)
                                                                                                          aCAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCt                                           <  2:3023252/149‑1 (MQ=255)
                                                                                                          aCAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCt                                           <  2:631756/149‑1 (MQ=255)
                                                                                                               cACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCgg                                      >  1:1315633/1‑149 (MQ=255)
                                                                                                               cACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCgg                                      >  1:3373722/1‑149 (MQ=255)
                                                                                                               cACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCgg                                      >  1:984500/1‑149 (MQ=255)
                                                                                                                aCGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGa                                     <  2:280671/149‑1 (MQ=255)
                                                                                                                   ccTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTg                                  <  2:1995819/149‑1 (MQ=255)
                                                                                                                   ccTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTg                                  >  2:1552675/1‑149 (MQ=255)
                                                                                                                     ttCGGCGTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGa                                <  2:3118475/149‑1 (MQ=255)
                                                                                                                       cGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATg                              >  2:3009102/1‑149 (MQ=255)
                                                                                                                           gTTTTGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAAt                          >  2:3185396/1‑149 (MQ=255)
                                                                                                                            ttttGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATc                         >  1:1492725/1‑149 (MQ=255)
                                                                                                                              ttGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATcca                       <  1:700901/149‑1 (MQ=255)
                                                                                                                               tGCGCAATGTAGTTAGCATCGAGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATccac                      <  2:3236267/149‑1 (MQ=255)
                                                                                                                                    aaTGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGc                 <  2:3467220/149‑1 (MQ=255)
                                                                                                                                     aTGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGcc                >  2:4646492/1‑149 (MQ=255)
                                                                                                                                      tGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGccc               >  2:1350415/1‑148 (MQ=255)
                                                                                                                                         aGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAgcg            >  2:680667/1‑149 (MQ=255)
                                                                                                                                         aGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAgcg            >  2:3071940/1‑149 (MQ=255)
                                                                                                                                            tAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCt         >  2:925680/1‑149 (MQ=255)
                                                                                                                                             aGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCtt        >  1:494334/1‑149 (MQ=255)
                                                                                                                                             aGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCtt        <  2:1861150/149‑1 (MQ=255)
                                                                                                                                              gCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCttt       >  2:3860915/1‑149 (MQ=255)
                                                                                                                                                aTCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCg     >  1:3678793/1‑149 (MQ=255)
                                                                                                                                                  cgcgGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGa   <  1:3891048/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCGACCGCGAGCATACCCAGCGCGCGGGAGAATTTTTGCGTTAAATAGCCGGTTTCGTTACTCGACGCGGAAGCACACAGCATCCCGGTGGAGAGCCAGCGGTTAACAGTCACGCCTTCGGCGTTTTGCGCAATGTAGTTAGCATCGCGGTCTTCTTTCATCAGTTTGGCGATGCGATCAAACGCCTCTTCCCAACTGATTTGCTGCCATTTATCAGAACCTGGCGCACGGTATTCCGGAAACTTCAGACGGCTTTCGGAGTGGATGAAATCCACCAGGCCAGCGCCTTTCGGAC  >  CP009273/4078613‑4078908

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: