Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 211,900 T→A 7.3% V5802V (GTT→GTA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947211,9000TA 7.3% 274.5 / 12.4 96V5802V (GTT→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base A (4/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

CGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAG  >  NC_002947/211753‑212049
                                                                                                                                                    |                                                                                                                                                      
cGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTAc                                                                                                                                                      <  2:416830/150‑1 (MQ=35)
    ggTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACcgac                                                                                                                                                  >  2:551663/1‑150 (MQ=35)
        cTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGc                                                                                                                                              <  2:1911513/150‑1 (MQ=35)
         tGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGGCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGTTCAGCACCACGGTTACCGACGAGcc                                                                                                                                             >  1:2367261/1‑150 (MQ=255)
              cTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTa                                                                                                                                        <  1:2380626/150‑1 (MQ=35)
                  cACTGCCCCTGATAACGTCCACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCGACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGcc                                                                                                                                    <  1:256159/150‑1 (MQ=255)
                  cACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGcc                                                                                                                                    <  2:1381894/150‑1 (MQ=35)
                   aCTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTATCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCg                                                                                                                                   >  1:2391084/1‑150 (MQ=25)
                   aCTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCg                                                                                                                                   >  1:683293/1‑150 (MQ=35)
                    cTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGTCTCCCAGTTACCGACGAGCCAGGTACGCCgg                                                                                                                                  <  1:62843/150‑1 (MQ=11)
                     tGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCggg                                                                                                                                 <  1:825961/150‑1 (MQ=35)
                       ccccGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCa                                                                                                                               >  1:1387480/1‑150 (MQ=35)
                             aTAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGTACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAg                                                                                                                         <  1:297913/150‑1 (MQ=21)
                             aTAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAg                                                                                                                         <  2:719694/150‑1 (MQ=35)
                              tAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAgg                                                                                                                        >  1:299130/1‑150 (MQ=35)
                                  gTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCAGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGAc                                                                                                                    <  2:2038140/150‑1 (MQ=21)
                                   tCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTATCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGAcc                                                                                                                   <  2:2391084/150‑1 (MQ=25)
                                          cGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGCCCCCGGTCAGCACGTCGGTTACAGACGGTCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCaa                                                                                                            <  1:1069078/150‑1 (MQ=255)
                                           ggCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGCCCCCGGTCAGCACGTCGGTTACCGACGAGCAAGGTACGCCGGGCAACGAATGCTATCTGTTCAAg                                                                                                           <  2:1956064/150‑1 (MQ=11)
                                                 aaCGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAg                                                                                                                                           >  1:442252/1‑112 (MQ=37)
                                                 aaCGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAg                                                                                                                                           <  2:442252/112‑1 (MQ=38)
                                                        cTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTAACcgacgtgc                                                                                                                                              >  2:1319379/1‑99 (MQ=255)
                                                        cTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTAACCGACGTGc                                                                                                                                              <  1:1319379/102‑1 (MQ=255)
                                                         tGTCTTGATGTCGATCGCAACCGTTGGCGGTGCGTATGTCAGCAAGTTCGAGTAACTGCCACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTAGTCAAGGTCACCATCACGGc                                                                                            <  2:187543/150‑1 (MQ=11)
                                                           tCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCACGGACGGGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGACGACCTGGTCAAGGTCACCGTCACGGccc                                                                                          >  1:1354897/1‑149 (MQ=11)
                                                           tCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCg                                                                                          >  2:421012/1‑150 (MQ=39)
                                                               gATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCa                                                                                      <  1:2079081/150‑1 (MQ=39)
                                                               gATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCa                                                                                      <  2:981260/150‑1 (MQ=38)
                                                                     gATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTc                                                                                >  2:1846446/1‑150 (MQ=40)
                                                                     gATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTc                                                                                >  2:1318478/1‑150 (MQ=255)
                                                                           tACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGaa                                                                                                                          >  1:595726/1‑103 (MQ=38)
                                                                           tACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGaa                                                                                                                          <  2:595726/103‑1 (MQ=38)
                                                                           tACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTa                                                                            <  1:96932/148‑1 (MQ=39)
                                                                           tACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTa                                                                            >  2:96932/1‑148 (MQ=40)
                                                                           tACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGc                                                                          >  1:1710143/1‑150 (MQ=255)
                                                                           tACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGc                                                                          >  2:1072491/1‑150 (MQ=255)
                                                                            aCCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCATGTGGCATGGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGcc                                                                         <  2:2372256/150‑1 (MQ=14)
                                                                            aCCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGc                                                                                                                       <  2:1279021/105‑1 (MQ=38)
                                                                            aCCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGc                                                                                                                       >  1:1279021/1‑105 (MQ=35)
                                                                            aCCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGcc                                                                         <  1:1801785/150‑1 (MQ=40)
                                                                              cGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCGCGGCCGACCAGACTTCGGTAGCCga                                                                       >  2:1382926/1‑150 (MQ=255)
                                                                              cGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCga                                                                       <  2:299130/150‑1 (MQ=255)
                                                                               gTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCga                                                                       >  1:2188538/1‑149 (MQ=40)
                                                                               gTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCga                                                                       <  2:2188538/149‑1 (MQ=40)
                                                                                             tGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGCGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGAc                                                        <  2:1652824/150‑1 (MQ=25)
                                                                                               tCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCt                                                      >  1:439868/1‑150 (MQ=35)
                                                                                                 agCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTc                                                    <  1:698601/149‑1 (MQ=25)
                                                                                                   cAAGTTCAAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCAc                                                  >  1:1316318/1‑150 (MQ=25)
                                                                                                      gTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGt                                               >  2:546887/1‑150 (MQ=35)
                                                                                                      gTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGt                                               >  2:1903235/1‑150 (MQ=35)
                                                                                                       ttCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGtt                                              >  1:1756168/1‑150 (MQ=35)
                                                                                                       ttCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGtt                                              >  1:1820425/1‑150 (MQ=35)
                                                                                                          gAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           >  1:316060/1‑150 (MQ=35)
                                                                                                             cAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGAACTGGTCAAGGTCACCATCACGGGAGACCAGGCGGCGTTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTc                                        <  1:421012/150‑1 (MQ=11)
                                                                                                                ttGACAGCGGAAAAGACCTCGGTCAGCACGTCGGTGACCGACGAGCAAGGTACGCCGGGCGATGAGGGCGACCTAGTGAAGGTCACCATCACGGCCGACCAGTCTTCGGTAGCCGAGAACGGCAAACCGACCTTCACTGTTCGCGTCAAc                                     <  2:2367261/149‑1 (MQ=11)
                                                                                                                 tGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTc                                                                >  2:131811/1‑122 (MQ=35)
                                                                                                                 tGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTc                                                                <  1:131811/122‑1 (MQ=37)
                                                                                                                  gACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCa                                   <  1:546887/150‑1 (MQ=35)
                                                                                                                     acTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGc                                                                                            <  1:337000/90‑1 (MQ=33)
                                                                                                                     acTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGc                                                                                            >  2:337000/1‑90 (MQ=33)
                                                                                                                     acTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGcc                                >  2:1114300/1‑150 (MQ=35)
                                                                                                                     acTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGcc                                >  1:1571129/1‑150 (MQ=35)
                                                                                                                      cTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCg                               >  2:1645720/1‑150 (MQ=35)
                                                                                                                       tGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCa                                                                                      <  2:1034222/94‑1 (MQ=34)
                                                                                                                       tGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCa                                                                                      >  1:1034222/1‑94 (MQ=34)
                                                                                                                       tGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           <  2:237328/137‑1 (MQ=37)
                                                                                                                       tGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           >  1:237328/1‑137 (MQ=37)
                                                                                                                       tGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           <  2:236836/137‑1 (MQ=37)
                                                                                                                       tGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           >  1:236836/1‑137 (MQ=37)
                                                                                                                               aCCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAc                      >  2:95806/1‑150 (MQ=34)
                                                                                                                               aCCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAc                      >  2:2269535/1‑150 (MQ=35)
                                                                                                                               aCCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAc                      >  1:94903/1‑150 (MQ=34)
                                                                                                                                 cccGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCGGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCCACCAACCGCTGGCTCACGa                    >  2:379202/1‑150 (MQ=255)
                                                                                                                                 cccGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCga                                                                       >  1:560433/1‑99 (MQ=39)
                                                                                                                                 cccGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCga                                                                       <  2:560433/99‑1 (MQ=39)
                                                                                                                                  ccGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTgg                           >  1:1567549/1‑142 (MQ=35)
                                                                                                                                  ccGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTgg                           <  2:1567549/142‑1 (MQ=35)
                                                                                                                                  ccGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGAc                   >  2:1264427/1‑150 (MQ=35)
                                                                                                                                    ggTCAGCAAGTCCGTTCCCGACGCTCTAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCAACATCACTGCCGACCAGACTTCGGTAGCCGACAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCt                 <  2:559941/150‑1 (MQ=11)
                                                                                                                                     gTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTg                <  1:2269535/150‑1 (MQ=35)
                                                                                                                                     gTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTg                <  2:1387480/150‑1 (MQ=35)
                                                                                                                                       cAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGt              <  2:1571129/150‑1 (MQ=35)
                                                                                                                                         cgCACGTCGGTTACCGACGAGCCAGGTACGCCGAGCAATGTAGGCTACCGGGTCAAGGTCACCATCACCGGCCGGCACACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTc             <  1:539862/149‑1 (MQ=1)
                                                                                                                                            tcGTCGGTAACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           <  1:1065917/116‑1 (MQ=255)
                                                                                                                                            tcGTCGGTAACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCAc                                           >  2:1065917/2‑117 (MQ=255)
                                                                                                                                            tcGTCGGTAACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTc             >  2:527317/2‑147 (MQ=25)
                                                                                                                                            tcGTCGGTAACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTc             <  1:527317/146‑1 (MQ=25)
                                                                                                                                            tcGTCGGTAACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTg          >  2:1418300/2‑150 (MQ=25)
                                                                                                                                             cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa         <  2:1820425/150‑1 (MQ=35)
                                                                                                                                             cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa         <  1:1114300/150‑1 (MQ=35)
                                                                                                                                              gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGag        >  2:1824307/1‑149 (MQ=25)
                                                                                                                                              gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGATAGCCGAGAAc                                                                   <  1:1439510/91‑1 (MQ=255)
                                                                                                                                              gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGATAGCCGAGAAc                                                                   >  2:1439510/1‑91 (MQ=255)
                                                                                                                                                cGGTTACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAccc      <  1:1072491/150‑1 (MQ=25)
                                                                                                                                                  gTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGa                                                                                    <  2:1066093/70‑1 (MQ=37)
                                                                                                                                                  gTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGa                                                                                    >  1:1066093/1‑70 (MQ=37)
                                                                                                                                                   ttACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa   <  1:95806/150‑1 (MQ=35)
                                                                                                                                                   ttACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa   <  1:1903235/150‑1 (MQ=35)
                                                                                                                                                    tACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGagccagcg  <  1:1824307/150‑8 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAG  >  NC_002947/211753‑212049

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: