Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,306 T→C 19.1% N604N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3060TC19.1% 138.2 / 24.8 68N604N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base C (4/9);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

GCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGC  >  NC_002947/196161‑196455
                                                                                                                                                 |                                                                                                                                                     
gCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTg                                                                                                                                                   >  2:700256/1‑150 (MQ=255)
       tttCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGACGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGAc                                                                                                                                            >  1:364164/1‑150 (MQ=255)
                aCTGACCCCGAACCCGACGCCGCCTCAGACCACGATCAACGACTCGGGCCACACCACTACCGTGACCCTGACGGCGAGCCCGTCCGTCACCGAACGCCGCGTGATCACTTACACCGTGCCCCGGAGCAATCCTGCCCAGACGGCGCTGAc                                                                                                                                   >  2:192842/1‑150 (MQ=11)
                            cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTc                                                                                                                       <  1:446011/150‑1 (MQ=255)
                            cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTc                                                                                                                       >  1:72538/1‑150 (MQ=255)
                              cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCa                                                                                                                     >  1:669560/1‑150 (MQ=255)
                                 cGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACACTGACGGCGAGCCCGTCGGTCACCGGAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                                                  <  1:317516/150‑1 (MQ=255)
                                   ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGc                                                                                                                <  2:72538/150‑1 (MQ=255)
                                    cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc                                                                                                                            <  2:65078/137‑1 (MQ=255)
                                    cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc                                                                                                                            >  1:65078/1‑137 (MQ=255)
                                             ccACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                                                                                                                                        >  1:28799/1‑116 (MQ=38)
                                             ccACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                                                                                                                                        <  2:28799/116‑1 (MQ=38)
                                                 gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGAc                                                                                                                                            <  1:839897/108‑1 (MQ=38)
                                                 gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGAc                                                                                                                                            >  2:839897/1‑108 (MQ=38)
                                                 gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTTACCCTGTCCAACGGCCAAACCATTACCGt                                                                                                  >  2:375070/1‑150 (MQ=255)
                                                 gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGt                                                                                                  >  1:186317/1‑150 (MQ=40)
                                                 gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCCGTCCAACGGCCAAACCATTACCGt                                                                                                  >  1:618058/1‑150 (MQ=255)
                                                  aTCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACCAAGTGACCCTGAGCAATCCTGCCCAGATCTCGTTGACATTGTCCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                 <  1:700256/150‑1 (MQ=14)
                                                     aaCGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGaa                                                                                              >  1:846868/1‑150 (MQ=255)
                                                     aaCGACTCGGTCGACACCACTACCACGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGaa                                                                                              >  2:622375/1‑150 (MQ=255)
                                                      aCGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAg                                                                                             >  1:542177/1‑150 (MQ=255)
                                                       cGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCAGTCACCGAAGGCGGCGTGATCACTTACACCGAGCACATGACGATTGATCGACTTACGCCGAATACAGTGACCCTGTCCAACGGCCAAACCATTACGGTTGAAGc                                                                                            >  1:101321/1‑150 (MQ=14)
                                                           tCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGc                                                                                        >  2:289476/1‑150 (MQ=40)
                                                           tCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGc                                                                                        >  1:678089/1‑150 (MQ=255)
                                                              gTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaaa                                                                                     >  1:259953/1‑150 (MQ=40)
                                                              gTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaaa                                                                                     >  1:92113/1‑150 (MQ=255)
                                                               tCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaaaa                                                                                    >  2:436447/1‑150 (MQ=255)
                                                                  acacCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAAccc                                                                                 >  2:5246/1‑150 (MQ=40)
                                                                   caccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCa                                                                                >  1:291885/1‑150 (MQ=40)
                                                                     ccaccaCCGCGTCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGTTTACTTACACCGTGACCCTGAGCAAGCCTTCCACGCCGCCGGTGACCGTGACCCTGTCCAACGCTCTACCGATTACCGTTGAAGCCGGCAAAACCCAgg                                                                              <  2:658319/150‑1 (MQ=11)
                                                                     ccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGACCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAgg                                                                              <  1:60662/150‑1 (MQ=255)
                                                                      cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCaaa                                                                                                            <  2:842953/119‑1 (MQ=38)
                                                                      cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCaaa                                                                                                            >  1:842953/1‑119 (MQ=37)
                                                                      cacTACCACGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAggg                                                                             <  1:622375/150‑1 (MQ=255)
                                                                       acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                                               >  1:128423/1‑147 (MQ=17)
                                                                       acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                                               <  2:128423/147‑1 (MQ=14)
                                                                       acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGc                                                                            <  1:436447/150‑1 (MQ=38)
                                                                        cTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGc                                                                                                                                                  <  2:191630/79‑1 (MQ=39)
                                                                        cTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGc                                                                                                                                                  >  1:191630/1‑79 (MQ=39)
                                                                             gcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTc                                                                      <  1:289476/150‑1 (MQ=35)
                                                                             gcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTc                                                                      <  1:818013/150‑1 (MQ=255)
                                                                                   cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTc                                                                <  2:678089/150‑1 (MQ=35)
                                                                                     gACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCa                                                              <  2:31849/150‑1 (MQ=35)
                                                                                             gCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAgg                                                                              >  1:429769/1‑126 (MQ=255)
                                                                                             gCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAgg                                                                              <  2:429769/126‑1 (MQ=255)
                                                                                               ccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCg                                                    <  2:714006/150‑1 (MQ=11)
                                                                                                 gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                 <  2:549761/103‑1 (MQ=32)
                                                                                                 gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                 >  1:549761/1‑103 (MQ=32)
                                                                                                     tcgaCCGAAGGCAGCGTGATCACTTACCCCGTGCCCCTGAGCAATCCTTCCCAGTCGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                              <  2:201544/147‑1 (MQ=11)
                                                                                                       cACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGACGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGt                                            <  1:450170/150‑1 (MQ=11)
                                                                                                        aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAgg                                                                              <  2:416634/115‑1 (MQ=255)
                                                                                                        aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAgg                                                                              >  1:416634/1‑115 (MQ=11)
                                                                                                          cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTa                                         <  2:752050/150‑1 (MQ=11)
                                                                                                                 ggcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac                                  <  1:761637/150‑1 (MQ=35)
                                                                                                                  gcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                 <  2:542177/150‑1 (MQ=35)
                                                                                                                   cgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACgg                                <  2:92113/150‑1 (MQ=35)
                                                                                                                    gTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                   <  2:45526/66‑1 (MQ=37)
                                                                                                                    gTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                   >  1:45526/1‑66 (MQ=37)
                                                                                                                     tGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                              >  2:639539/1‑150 (MQ=35)
                                                                                                                      gATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                             <  2:291885/150‑1 (MQ=35)
                                                                                                                        tCACTTACCCCGCGACCCTGACCAACCCTGCCCAGACGCCGGTCACCGCGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGAGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                           <  1:581168/150‑1 (MQ=11)
                                                                                                                            ttACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                       >  2:234910/1‑150 (MQ=11)
                                                                                                                                        ccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTac                                        >  1:759246/1‑121 (MQ=32)
                                                                                                                                        ccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTac                                        <  2:759246/121‑1 (MQ=33)
                                                                                                                                        ccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                 <  2:593618/128‑1 (MQ=11)
                                                                                                                                        ccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                 >  1:593618/1‑128 (MQ=11)
                                                                                                                                         cTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATc          <  1:305960/150‑1 (MQ=38)
                                                                                                                                           gAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAAGACAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCga        >  1:104084/1‑150 (MQ=255)
                                                                                                                                                 tCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGc  >  1:755029/1‑150 (MQ=38)
                                                                                                                                                 |                                                                                                                                                     
GCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGC  >  NC_002947/196161‑196455

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: