Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,490,819 C→A 8.8% D349Y (GAT→TAT) ‡ PP_3982 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,490,8190CA 8.8% 203.4 / 8.3 68D349Y (GAT→TAT) ‡PP_3982hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base A (4/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

TCCGCTCCAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTACCTCTTGGAGC  >  NC_002947/4490672‑4490967
                                                                                                                                                   |                                                                                                                                                    
tCCGCTCCAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGAGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTCGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCg                                                                                                                                                     >  2:248870/1‑149 (MQ=255)
   gCTCCAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCtg                                                                                                                                                  >  3:176800/1‑149 (MQ=255)
   gCTCCAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCtg                                                                                                                                                  <  6:357124/149‑1 (MQ=255)
       cAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCtgtg                                                                                                                                                <  5:207874/147‑1 (MQ=255)
       cAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGtt                                                                                                                                              >  5:125032/1‑149 (MQ=255)
        aaGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTg                                                                                                                                             <  8:186044/149‑1 (MQ=255)
          gCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACACAGCGATGTTGCt                                                                                                                                           <  2:302125/149‑1 (MQ=255)
                   aGCCTTTTGGTAAGGGTCGACATGTACGTGATGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTTGTAAGCAACATCGCTGTGTTGCTGTGCAAg                                                                                                                                    <  7:232069/147‑1 (MQ=255)
                   aGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAgg                                                                                                                                   <  4:353546/148‑1 (MQ=255)
                      cTTTTGGTAAGGGTCGACATGTACGTTAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGtt                                                                                                                               >  1:5301/1‑149 (MQ=255)
                                   tCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCtt                                                                                                                   >  5:85559/1‑148 (MQ=255)
                                      aCATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAAc                                                                                                               <  6:125032/149‑1 (MQ=255)
                                      aCATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAAc                                                                                                               <  1:248870/149‑1 (MQ=255)
                                                      ccccTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCa                                                                                               >  6:63407/1‑149 (MQ=255)
                                                       cccTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGc                                                                                                   <  2:276985/144‑1 (MQ=255)
                                                       cccTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCAt                                                                                              >  5:28536/1‑149 (MQ=255)
                                                       cccTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCAt                                                                                              <  4:322361/149‑1 (MQ=255)
                                                       cccTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTAAACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGc                                                                                                   >  1:276985/1‑144 (MQ=255)
                                                        ccTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCAt                                                                                              >  5:290324/1‑148 (MQ=255)
                                                                         gCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTTCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAggcg                                                                             <  5:296387/148‑1 (MQ=255)
                                                                         gCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAggcgg                                                                            >  3:358164/1‑149 (MQ=255)
                                                                         gCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAggcgg                                                                            <  7:103592/149‑1 (MQ=255)
                                                                          cACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAggcggc                                                                           >  8:23866/1‑149 (MQ=255)
                                                                             gCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCa                                                                         >  1:100734/1‑148 (MQ=255)
                                                                               tGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCaaaa                                                                      <  2:5301/149‑1 (MQ=255)
                                                                                       gATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGc                                                              >  8:293297/1‑149 (MQ=255)
                                                                                        aTCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCg                                                             <  4:176800/149‑1 (MQ=255)
                                                                                           gggTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGc                                                          <  2:100734/149‑1 (MQ=255)
                                                                                             gTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCt                                                         >  8:296440/1‑148 (MQ=255)
                                                                                              tGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCaaaa                                                                      >  6:370208/1‑134 (MQ=255)
                                                                                              tGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCaaaa                                                                      <  5:370208/134‑1 (MQ=255)
                                                                                              tGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTaa                                                       >  5:18304/1‑149 (MQ=255)
                                                                                                aaaTGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATg                                                     <  6:18304/149‑1 (MQ=255)
                                                                                                aaaTGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATg                                                     >  5:134095/1‑149 (MQ=255)
                                                                                                   tGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCg                                                  <  1:340417/149‑1 (MQ=255)
                                                                                                     gCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCaa                                                                                                                                     <  3:122313/64‑1 (MQ=255)
                                                                                                     gCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCaa                                                                                                                                     >  4:122313/1‑64 (MQ=255)
                                                                                                       cTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTa                                                >  3:352024/1‑147 (MQ=255)
                                                                                                          tACACGCTGAACAGCATTGAACAGGGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCa                                           >  8:14083/1‑149 (MQ=255)
                                                                                                              cGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATggg                                       <  6:85559/149‑1 (MQ=255)
                                                                                                              cGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATggg                                       <  6:290324/149‑1 (MQ=255)
                                                                                                               gCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGc                                      >  5:203705/1‑149 (MQ=255)
                                                                                                                   aaCAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAg                                  >  8:202966/1‑149 (MQ=255)
                                                                                                                     cAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCAGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCTGGCTAATGcgg                                                  >  8:334474/1‑129 (MQ=255)
                                                                                                                        cATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCaaaa                                                                      >  2:254717/1‑108 (MQ=255)
                                                                                                                        cATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGa                             >  8:287795/1‑149 (MQ=255)
                                                                                                                        cATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTACAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCaaaa                                                                      <  1:254717/108‑1 (MQ=255)
                                                                                                                         aTTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGAt                            <  4:180088/149‑1 (MQ=255)
                                                                                                                          ttGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGAtt                           >  8:243152/1‑149 (MQ=255)
                                                                                                                          ttGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGAtt                           <  7:296440/149‑1 (MQ=255)
                                                                                                                             aaCAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAg                        >  1:342512/1‑149 (MQ=255)
                                                                                                                                aaaGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAgcg                      >  5:9567/1‑148 (MQ=255)
                                                                                                                                  aGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCaa                   <  4:161444/149‑1 (MQ=255)
                                                                                                                                  aGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACAGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCGGTAACGCATGGGCGGAGGCTGAGTGAgcgc                     >  8:322142/1‑147 (MQ=255)
                                                                                                                                    cttgcacaGCAACACAGCGATGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGc                 >  3:334469/8‑149 (MQ=255)
                                                                                                                                         tAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTAc            <  6:28536/149‑1 (MQ=255)
                                                                                                                                         tAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTAc            <  4:352024/149‑1 (MQ=255)
                                                                                                                                         acaGCAACACAGCGATGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTAc            >  8:218683/3‑149 (MQ=255)
                                                                                                                                         acaGCAACACAGCGATGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTAc            >  1:45041/3‑149 (MQ=255)
                                                                                                                                          caGCAACACAGCGATGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAgc                       <  8:239004/136‑1 (MQ=255)
                                                                                                                                          caGCAACACAGCGATGTTGCTGTGCAAGGTGTTGCGACAGAGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAgc                       >  7:239004/2‑137 (MQ=255)
                                                                                                                                             cAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGc                                                                                                   <  8:163849/58‑1 (MQ=255)
                                                                                                                                             cAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGc                                                                                                   >  7:163849/1‑58 (MQ=255)
                                                                                                                                              aaCATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAg                                  <  3:194234/122‑1 (MQ=255)
                                                                                                                                              aaCATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAg                                  >  4:194234/1‑122 (MQ=255)
                                                                                                                                               aCATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTACCTCTTg      <  6:134095/149‑1 (MQ=255)
                                                                                                                                               aCATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTACCTCTTg      >  6:145388/1‑149 (MQ=255)
                                                                                                                                                 aTCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTACCTCTTGGa    <  7:14083/149‑1 (MQ=255)
                                                                                                                                                   cGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTACCTCTTGGAGc  >  2:253725/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCGCTCCAAGCCTTGCTCAGCCTTTTGGTAAGGGTCGACATGTACGTGAGGAACCCCTTGCAGTGGCGTGTCGCACGCTGGGCAGCGATCGGGTGAAATGGCCTGTACACGCTGAACAGCATTGAACAAAGCTTGGTAAGCAACATCGCTGTGTTGCTGTGCAAGGTGTTGCGACAGCGCTTCAACCGCAGCGTGGGCCGCATCAGCTGCCTGATAGGCGGCCAAAAGTTTTTGCGGGCTAATGCCGTACAGCATGGGCGGAGGCTGATTGAGCGCAAGCTCTACCTCTTGGAGC  >  NC_002947/4490672‑4490967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑