Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,735,398 C→A 5.3% P735P (CCG→CCT sucA ← 2‑oxoglutarate dehydrogenase subunit E1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,735,3980CA 5.3% 602.4 / 15.1 180P735P (CCG→CCTsucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base C (79/91);  new base A (2/8);  total (81/99)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

GGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCG  >  NC_002947/4735253‑4735546
                                                                                                                                                 |                                                                                                                                                    
ggACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGcc                                                                                                                                                  <  2:4447553/150‑1 (MQ=255)
 gACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                                                                                 <  1:1033303/150‑1 (MQ=255)
 gACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                                                                                 >  1:5398210/1‑150 (MQ=255)
   cTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCt                                                                                                                                               <  1:4821931/150‑1 (MQ=255)
    tACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTg                                                                                                                                              >  1:1376279/1‑150 (MQ=255)
    tACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTg                                                                                                                                              >  1:3884467/1‑150 (MQ=255)
      ccAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACTGCTGCAGGTAACGCTCCTGACGCGCTGAGGAGTGCTCAGGCCCCTgac                                                                                                                                            <  2:1872519/150‑3 (MQ=255)
      ccAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGcc                                                                                                                                            >  1:1262850/1‑150 (MQ=255)
      ccAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGcc                                                                                                                                            >  2:2767751/1‑150 (MQ=255)
       cAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGccc                                                                                                                                           >  1:3712869/1‑150 (MQ=255)
         gCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCtt                                                                                                                                         >  2:4929000/1‑150 (MQ=255)
          cGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTc                                                                                                                                        <  2:3157815/150‑1 (MQ=255)
           ggCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCa                                                                                                                                       <  1:3260750/150‑1 (MQ=255)
             cTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATa                                                                                                                                     <  2:4103187/150‑1 (MQ=255)
               tGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGc                                                                                                                                   <  1:102581/150‑1 (MQ=255)
                gcgcAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGcc                                                                                                                                  >  2:2554163/1‑150 (MQ=255)
                gcgcAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGcc                                                                                                                                  >  2:769504/1‑150 (MQ=255)
                    aGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCgtgt                                                                                                                              >  2:108668/1‑150 (MQ=255)
                    aGCGGGCGAATGACCTGCCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCACCGGCCCCTGCCCTTCATAGCCgtgt                                                                                                                              <  2:2380951/150‑1 (MQ=255)
                           gAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGGGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagca                                                                                                                       <  2:1791181/150‑1 (MQ=255)
                           gAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagca                                                                                                                       >  2:20010/1‑150 (MQ=255)
                           gAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagca                                                                                                                       >  1:3977075/1‑150 (MQ=255)
                             aTGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagcagc                                                                                                                     <  2:1273338/150‑1 (MQ=255)
                             aTGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagcagc                                                                                                                     <  2:5103531/150‑1 (MQ=255)
                                aCCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTc                                                                                                                                        <  2:2688531/128‑1 (MQ=255)
                                aCCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATg                                                                                                                  <  1:2686976/150‑1 (MQ=255)
                                 ccTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATgg                                                                                                                 >  2:177680/1‑150 (MQ=255)
                                   tGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTc                                                                                                               >  1:241692/1‑150 (MQ=255)
                                     gCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAg                                                                                                             >  2:1359461/1‑150 (MQ=255)
                                       gACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGAc                                                                                                           <  2:4794954/150‑1 (MQ=255)
                                        aCGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGAcc                                                                                                          <  2:474669/150‑1 (MQ=255)
                                          gcagcaGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCatctcga                                                                                                        <  2:2618389/150‑7 (MQ=255)
                                          gcagcaGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGc                                                                                                        >  2:5305828/1‑150 (MQ=255)
                                          gcagcaGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGc                                                                                                        >  2:2235241/1‑150 (MQ=255)
                                          gcagcaGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGc                                                                                                        >  2:1758032/1‑150 (MQ=255)
                                                ggtggtAAATCTGTGCCGGCGTCGTCGGTATGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGAc                                                                                                           >  2:3669008/1‑141 (MQ=255)
                                                ggtggtAAATCTGTGCCGGCATCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggc                                                                                                  <  2:3480881/150‑1 (MQ=255)
                                                  tggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGTGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcgg                                                                                                <  1:465018/150‑1 (MQ=255)
                                                  tggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcgg                                                                                                <  1:325479/150‑1 (MQ=255)
                                                  tggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcgg                                                                                                >  1:75125/1‑150 (MQ=255)
                                                   ggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcggc                                                                                               >  2:5163537/1‑150 (MQ=255)
                                                   ggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcggc                                                                                               <  1:2090555/150‑1 (MQ=255)
                                                   ggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcggc                                                                                               <  1:3060002/150‑1 (MQ=255)
                                                         tCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACt                                                                                         >  2:102926/1‑150 (MQ=255)
                                                           tgtgCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTtg                                                                                       >  1:3629441/1‑150 (MQ=255)
                                                             tgCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTtgtg                                                                                     >  1:5103233/1‑150 (MQ=255)
                                                                  gCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTcac                                                                                <  1:2854378/150‑1 (MQ=255)
                                                                       gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                           <  1:2969693/150‑1 (MQ=255)
                                                                       gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                           <  2:2180283/150‑1 (MQ=255)
                                                                       gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                           <  2:3884467/150‑1 (MQ=255)
                                                                        tcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATg                                                                                                                  >  1:4448240/1‑110 (MQ=255)
                                                                        tcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTgg                                                                          >  1:3002418/1‑150 (MQ=255)
                                                                         cgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGt                                                                         <  2:3737325/150‑1 (MQ=255)
                                                                         cgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGt                                                                         <  1:5193078/150‑1 (MQ=255)
                                                                          ggTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtg                                                                        >  2:1428109/1‑150 (MQ=255)
                                                                           gTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtga                                                                       <  1:586926/150‑1 (MQ=255)
                                                                               gCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGaa                                                                   <  2:4668851/150‑1 (MQ=255)
                                                                               gCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGaa                                                                   >  2:5112224/1‑150 (MQ=255)
                                                                                  gACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTg                                                                <  2:1774293/150‑1 (MQ=255)
                                                                                      tGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCg                                                            >  2:4169229/1‑150 (MQ=255)
                                                                                       ggATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGAc                                                                                                           <  2:3924034/102‑1 (MQ=255)
                                                                                        gATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAt                                                          >  2:1888992/1‑150 (MQ=255)
                                                                                          tgttgtGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATca                                                        <  2:2529786/150‑1 (MQ=255)
                                                                                          tgttgtGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATca                                                        >  1:1621973/1‑150 (MQ=255)
                                                                                              gtgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATcaccac                                                    >  1:3819988/1‑150 (MQ=255)
                                                                                               tgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACt                                                   >  1:3667384/1‑150 (MQ=255)
                                                                                               tgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACt                                                   >  1:3801463/1‑150 (MQ=255)
                                                                                                gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACtt                                                  <  1:1641719/150‑1 (MQ=255)
                                                                                                gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACtt                                                  <  1:1625625/150‑1 (MQ=255)
                                                                                                gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACtt                                                  <  2:4600210/150‑1 (MQ=255)
                                                                                                 cTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTg                                                 <  1:5282008/150‑1 (MQ=255)
                                                                                                 cTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTg                                                 >  2:1719661/1‑150 (MQ=255)
                                                                                                  tCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgc                                                <  2:1267132/150‑1 (MQ=255)
                                                                                                   cGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgcg                                               >  2:2347157/1‑150 (MQ=255)
                                                                                                        cacaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCg                                          >  1:4295423/1‑150 (MQ=255)
                                                                                                         acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCTTAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                         >  2:1183413/1‑150 (MQ=255)
                                                                                                         acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                         <  2:4674281/150‑1 (MQ=255)
                                                                                                         acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                         >  2:2291638/1‑150 (MQ=255)
                                                                                                             cTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGc                                     <  2:3305537/150‑1 (MQ=255)
                                                                                                                cAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGaa                                  <  2:3648755/150‑1 (MQ=255)
                                                                                                                   gTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTc                               <  1:696062/150‑1 (MQ=255)
                                                                                                                     aaCGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGc                             >  1:3539402/1‑150 (MQ=255)
                                                                                                                          tCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAAc                        >  2:1994080/1‑150 (MQ=255)
                                                                                                                            cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                      >  1:1552298/1‑150 (MQ=255)
                                                                                                                            cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                      <  1:4223139/150‑1 (MQ=255)
                                                                                                                            cAGACGCGCGGAGGACTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGACAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                      <  2:5398210/150‑1 (MQ=255)
                                                                                                                             cgACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCAGAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTgg                     <  1:3754091/149‑1 (MQ=255)
                                                                                                                              aggggccgGAGCACTCCTCCGCCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTggg                    >  2:4480612/7‑150 (MQ=255)
                                                                                                                                    cgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTccc              <  1:2699407/150‑1 (MQ=255)
                                                                                                                                    cgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTccc              <  1:2808173/150‑1 (MQ=255)
                                                                                                                                    cgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTccc              <  2:2079284/150‑1 (MQ=255)
                                                                                                                                     ggaggaGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCa             <  2:3539402/150‑1 (MQ=255)
                                                                                                                                         gaatGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGaga         <  2:4497488/147‑3 (MQ=255)
                                                                                                                                           gTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCAc       >  2:4985031/1‑150 (MQ=255)
                                                                                                                                              cTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGa           >  2:2914241/1‑143 (MQ=255)
                                                                                                                                              cTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGa           <  1:2914241/143‑1 (MQ=255)
                                                                                                                                               tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc   >  2:4010740/1‑150 (MQ=255)
                                                                                                                                               tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc   >  1:3880366/1‑150 (MQ=255)
                                                                                                                                               tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc   >  1:1861684/1‑150 (MQ=255)
                                                                                                                                                 cGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtgat                                                                      >  1:2192938/1‑81 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCG  >  NC_002947/4735253‑4735546

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: