Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 930,220 C→A 5.2% K5135N (AAG→AAT PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947930,2200CA 5.2% 300.9 / 4.0 95K5135N (AAG→AATPP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (46/44);  new base A (3/2);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

TTGATGCCGTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACTGGCCGTTGTTCT  >  NC_002947/930071‑930368
                                                                                                                                                     |                                                                                                                                                    
ttGATGCCGTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGccc                                                                                                                                                      <  2:1459282/150‑1 (MQ=255)
  gATGCCGTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCtt                                                                                                                                                    >  1:2147742/1‑150 (MQ=255)
   aTGCCGTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTg                                                                                                                                                   >  2:1415372/1‑150 (MQ=255)
     gCCGTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCAACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGcc                                                                                                                                                 <  2:1113699/150‑1 (MQ=255)
      tgctgacggtggtgtccaccgatgatgccggcaacgcttcggtggctgtcccttaccagaccaacgacacctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                >  2:547568‑M1/128‑150 (MQ=255)
        gTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                <  2:1226932/148‑1 (MQ=255)
        gTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                >  1:1226932/1‑148 (MQ=255)
        gTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCATCGTCTCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGt                                                                                                                                              >  1:1620724/1‑150 (MQ=255)
         tgacggtggtgtccaccgatgatgccggcaacgcttcggtggctgtcccttaccagaccaacgacacctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGACCTTGCCAGTc                                                                                                                                             >  1:552450‑M1/125‑150 (MQ=255)
         tCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTc                                                                                                                                             <  1:1415372/150‑1 (MQ=255)
            gACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGc                                                                                                                                          <  2:1025379/150‑1 (MQ=255)
                   tCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGa                                                                                                                                   >  1:1285054/1‑150 (MQ=255)
                   tCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGa                                                                                                                                   >  2:2206858/1‑150 (MQ=255)
                     aCTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCCACGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 <  1:1234656/150‑1 (MQ=255)
                     aCTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 >  1:1239637/1‑150 (MQ=255)
                       ttCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  >  1:865756/1‑147 (MQ=255)
                       ttCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  <  2:865756/147‑1 (MQ=255)
                       ttCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 <  2:1386435/148‑1 (MQ=255)
                       ttCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 >  1:1386435/1‑148 (MQ=255)
                       ttCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 <  2:534279/148‑1 (MQ=255)
                       ttCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 >  1:534279/1‑148 (MQ=255)
                       accgatgatgccggcaacgcttcggtggctgtcccttaccagaccaacgacacctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                <  2:77285‑M1/23‑1 (MQ=255)
                       accgatgatgccggcaacgcttcggtggctgtcccttaccagaccaacgacacctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                >  1:77285‑M1/111‑133 (MQ=255)
                            gagTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGtgttgt                                                                                                                          <  2:1285054/150‑1 (MQ=255)
                               tGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGa                                                                                                                       >  1:189465/1‑150 (MQ=255)
                               tGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTCGCCAGCGCAATCCAACACGCCCGTACAGATCACGTCTCCGTCGAAGCGGATGGCTATCGGTGCATTGGCTTCGCCCTTGCCAGCCAGCACGTTGAAGGTGTTGTTGa                                                                                                                       >  2:1478068/1‑150 (MQ=255)
                                   gagcgagcCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 >  2:396830/1‑136 (MQ=255)
                                   gagcgagcCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 <  1:396830/136‑1 (MQ=255)
                                   gagcgagcCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGc                                                                                                                   >  2:2328278/1‑150 (MQ=255)
                                       gagcCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAgg                                                                                                               <  2:2376067/150‑1 (MQ=255)
                                         gcCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGtt                                                                                                             >  2:2088384/1‑150 (MQ=255)
                                         gcCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGtt                                                                                                             <  2:2147742/150‑1 (MQ=255)
                                           ccACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTgg                                                                                                           >  1:1900155/1‑150 (MQ=255)
                                              ccGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCa                                                                                                        <  1:2088384/150‑1 (MQ=255)
                                                      actaGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCg                                                                                                >  1:2402359/4‑150 (MQ=255)
                                                      aCCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                <  2:187153/102‑1 (MQ=255)
                                                      aCCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCa                                                                                                                                                >  1:187153/1‑102 (MQ=255)
                                                       ccAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCgg                                                                                               >  1:1434763/1‑150 (MQ=255)
                                                       ccAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCgg                                                                                               <  2:1049511/150‑1 (MQ=255)
                                                        cAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCggg                                                                                              >  2:2251196/1‑150 (MQ=255)
                                                          gtgtTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGtgttgt                                                                                                                          <  2:1511001/120‑1 (MQ=255)
                                                          gtgtTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGtgttgt                                                                                                                          >  1:1511001/1‑120 (MQ=255)
                                                               aTTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAg                                                                                       >  2:236165/1‑150 (MQ=255)
                                                                tttGGAAGTTGCCAGCGCCATCCACCACGCCGATGCCTGCCGCGACGCCGGCGAAGCTGCCAGTGACCGTCGCCTCGGCTACCTGATCGTCAGTCAGCACGATGAAGGTGTTGTTGATCGCCAGGTTGGAGACCAGGACGGGGTGCGAgg                                                                                      >  1:2350369/1‑150 (MQ=255)
                                                                 ttGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  >  1:1964095/1‑105 (MQ=255)
                                                                 ttGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  <  2:1964095/105‑1 (MQ=255)
                                                                 ttGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGccggtt                                                                                     <  2:2402359/150‑6 (MQ=255)
                                                                 tcggAAGTTGCCATCTTGATCTACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACTGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGgt                                                                                     <  2:1117961/148‑1 (MQ=255)
                                                                     acgacacctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTtgcta                                                                                                                                                <  2:1974401‑M1/23‑3 (MQ=255)
                                                                     acgacacctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTtgcta                                                                                                                                                >  1:1974401‑M1/65‑85 (MQ=255)
                                                                       gTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGtt                                                                               <  2:1076372/150‑1 (MQ=255)
                                                                            ctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTTGc                                                                                                                                                  >  2:568035‑M1/58‑78 (MQ=255)
                                                                            ctcgccacccgacctggtgaccaacctggcgatcaacaacaccttcaacgtgctgacTGGCATTGGCTTCGCCCTTGc                                                                                                                                                  <  1:568035‑M1/21‑1 (MQ=255)
                                                                            cAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  <  2:1548531/94‑1 (MQ=255)
                                                                            cAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  >  1:1548531/1‑94 (MQ=255)
                                                                             aGCGCCATCCACCACGCCGGTACCGATCACGTTTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTg                                                                         <  1:1217364/150‑1 (MQ=255)
                                                                             aGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTc                                                                                                 >  2:1978669/1‑126 (MQ=255)
                                                                             aGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTc                                                                                                 <  1:1978669/126‑1 (MQ=255)
                                                                                  cATCCACCCCGCCGGTACCGCTCACGCCTCCGTCGCCGCAGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAg                                                                    <  2:1330543/150‑1 (MQ=255)
                                                                                  cATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAg                                                                    <  1:2149630/150‑1 (MQ=255)
                                                                                   aTCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTAGGCCCGTGGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAgg                                                                   >  2:196142/1‑150 (MQ=255)
                                                                                     ccaccaCGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGa                                                                                        <  2:2067155/127‑1 (MQ=255)
                                                                                     ccaccaCGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGa                                                                                        >  1:2067155/1‑127 (MQ=255)
                                                                                       acctcGCCGGTACCGAGCACGTCTCCGTCGAAGCTGACGGTTACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACa                                                               <  2:181667/150‑1 (MQ=255)
                                                                                       accacGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACa                                                               >  1:1865092/1‑150 (MQ=255)
                                                                                       accacGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACa                                                               >  1:132999/1‑150 (MQ=255)
                                                                                          acGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACa                                                               >  1:1824999/1‑147 (MQ=255)
                                                                                          acGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACa                                                               <  2:1824999/147‑1 (MQ=255)
                                                                                          acGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGcc                                                            >  1:154481/1‑150 (MQ=255)
                                                                                            gCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCAc                                                          <  2:1865092/150‑1 (MQ=255)
                                                                                                gTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGaa                                                      <  2:1526517/150‑1 (MQ=255)
                                                                                                    cGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGt                                                  >  1:1754227/1‑150 (MQ=255)
                                                                                                       tCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGc                                               <  2:1239637/150‑1 (MQ=255)
                                                                                                            tctcCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 <  2:933452/63‑1 (MQ=255)
                                                                                                            tctcCGTCGAAGCTGACGGTGACAGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAg                                                                                                                                 >  1:933452/1‑63 (MQ=255)
                                                                                                             ctcCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  <  1:1092492/61‑1 (MQ=255)
                                                                                                             ctcCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGaa                                                                                                                                  >  2:1092492/1‑61 (MQ=255)
                                                                                                                cGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCATTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCgg                                                                                               >  1:1307136/1‑93 (MQ=255)
                                                                                                                cGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCATTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCgg                                                                                               <  2:1307136/93‑1 (MQ=255)
                                                                                                                 aacaccttcaacgtgctgacTGGCAAGGGCGAAGCCAATGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGa                                                                 <  1:1158428‑M1/102‑1 (MQ=255)
                                                                                                                 aacaccttcaacgtgctgacTGGCAAGGGCGAAGCCAATGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGa                                                                 >  2:1158428‑M1/21‑122 (MQ=255)
                                                                                                                     aaGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTg                                 <  2:154481/150‑1 (MQ=255)
                                                                                                                        cTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGCTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGac                              <  2:132999/150‑1 (MQ=255)
                                                                                                                          gACGGTGACCGTGGCATTGGCTTCGCCCTGGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGGCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGacac                            <  1:1478068/150‑1 (MQ=255)
                                                                                                                           aCGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGt                                                                                <  2:1772702/97‑1 (MQ=255)
                                                                                                                           aCGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGt                                                                                >  1:1772702/1‑97 (MQ=255)
                                                                                                                             ggTGACCTTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGAcaccac                         >  2:1658755/1‑150 (MQ=255)
                                                                                                                                   cGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAg                   >  2:2349084/1‑150 (MQ=255)
                                                                                                                                     TGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAg                                                                                       <  2:938912‑M1/80‑1 (MQ=255)
                                                                                                                                     TGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAg                                                                                       >  1:938912‑M1/1‑80 (MQ=255)
                                                                                                                                     TGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGAGCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCa                 >  2:1620724‑M1/1‑150 (MQ=255)
                                                                                                                                        CATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACt              <  1:1608397‑M1/150‑1 (MQ=255)
                                                                                                                                               ttCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACTGGCCgtt       <  2:454910/150‑1 (MQ=255)
                                                                                                                                               ttCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACTGGCCgtt       >  1:731677/1‑150 (MQ=255)
                                                                                                                                               ttCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACTGGCCgtt       >  1:2115340/1‑150 (MQ=255)
                                                                                                                                               ttCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGACACCGAAGCGTTGCCGGCATCATCGGGGGACACCATCGTCAGTATCTAGCCgtt       <  2:1900155/150‑1 (MQ=255)
                                                                                                                                                    ccTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACTGGCCGTTGTTCt  >  1:1808467/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGATGCCGTCGGACTGTGTCACTTCCAGAGTGGCGAGCGAGCCCACCGCTGGCACCAGTGTTATTTGGAAGTTGCCAGCGCCATCCACCACGCCGGTACCGATCACGTCTCCGTCGAAGCTGACGGTGACCGTGGCATTGGCTTCGCCCTTGCCAGTCAGCACGTTGAAGGTGTTGTTGATCGCCAGGTTGGTCACCAGGTCGGGTGGCGAGGTGTCGTTGGTCTGGTAAGGGACAGCCACCGAAGCGTTGCCGGCATCATCGGTGGACACCACCGTCAGCAACTGGCCGTTGTTCT  >  NC_002947/930071‑930368

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: