Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,622,205 A→T 6.9% K80I (AAA→ATA)  PP_5515 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,622,2050AT 6.9% 172.2 / 2.8 58K80I (AAA→ATA) PP_5515hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base T (2/2);  total (32/26)
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.92e-01

CCTGACCTGGCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCG  >  NC_002947/2622059‑2622352
                                                                                                                                                  |                                                                                                                                                   
cctgacctgGCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCaaa                                                                                                                                                    <  1:210204/148‑1 (MQ=255)
      ctgGCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGTCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGTtttt                                                                                                                                             >  2:306955/1‑149 (MQ=255)
      ctgGCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGAtttt                                                                                                                                             <  1:25824/149‑1 (MQ=255)
         gCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTttgtt                                                                                                                                          <  5:164713/149‑1 (MQ=255)
         gCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTttgtt                                                                                                                                          <  2:34711/149‑1 (MQ=255)
         gCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTttg                                                                                                                                            <  1:271201/147‑1 (MQ=255)
           acaTGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTttgttg                                                                                                                                         >  5:150086/1‑148 (MQ=255)
             aTGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCa                                                                                                                                       >  2:64012/1‑148 (MQ=255)
                     ccTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGTCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGTTTTTGTTGCAAAAAttt                                                                                                                                <  1:306955/147‑1 (MQ=255)
                     ccTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTga                                                                                                                              <  7:46920/149‑1 (MQ=255)
                     ccTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTg                                                                                                                               <  4:227512/148‑1 (MQ=255)
                      cTCACGGCTTACGTGCGGTACGGCTCCAATGAATCCACCTCGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGAttt                                                                                                                                              <  2:297174/132‑1 (MQ=255)
                      cTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAAttt                                                                                                                                >  5:373073/1‑146 (MQ=255)
                      cTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATACCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGAttt                                                                                                                                              >  1:297174/1‑132 (MQ=255)
                       tCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCAATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTga                                                                                                                              >  1:72369/1‑147 (MQ=255)
                         aCGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTTGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCCAATGATTTTGTTGCAAAAATTTgaga                                                                                                                            <  1:10746/147‑1 (MQ=255)
                         aCGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTgaga                                                                                                                            <  7:97284/147‑1 (MQ=255)
                           ggCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGccc                                                                                                                        <  6:373073/149‑1 (MQ=255)
                                  gTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCg                                                                                                                  >  5:265954/1‑148 (MQ=255)
                                      ggTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTga                                                                                                                              <  5:2371/132‑1 (MQ=255)
                                      ggTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAAtt                                                                                                             <  8:30872/149‑1 (MQ=255)
                                      ggTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTCGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTga                                                                                                                              >  6:2371/1‑132 (MQ=255)
                                       gTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAAtt                                                                                                             >  7:66299/1‑148 (MQ=255)
                                        tACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGa                                                                                                           <  1:64012/149‑1 (MQ=255)
                                                 aaTGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCcgc                                                                                                   >  4:323880/1‑148 (MQ=255)
                                                        ccaccTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTAccc                                                                                           >  8:145659/1‑149 (MQ=255)
                                                          accTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGTAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGc                                                                                                      <  8:344722/136‑1 (MQ=255)
                                                          accTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGTAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGc                                                                                                      >  7:344722/1‑136 (MQ=255)
                                                            cTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCgg                                                                                       <  8:66299/149‑1 (MQ=255)
                                                               ggCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGc                                                                                    <  6:265954/149‑1 (MQ=255)
                                                                  cATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTc                                                                                 >  5:55132/1‑149 (MQ=255)
                                                                         aGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGag                                                                          >  3:364773/1‑149 (MQ=255)
                                                                            ctctTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCGGCGACTACCCCCGGTGCCTCCAGCGagcag                                                                       >  4:222617/1‑149 (MQ=255)
                                                                                   tGGGTTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGCGCAGCCaaact                                                                >  4:9595/1‑147 (MQ=255)
                                                                                      gcgggccgggcGGCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCg                                                             <  6:55132/146‑1 (MQ=255)
                                                                                       ctggcctggcGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGc                                                            >  6:295930/1‑149 (MQ=255)
                                                                                             tggcGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCATATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCctct                                                       >  6:72614/1‑148 (MQ=255)
                                                                                             tggcGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCATATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCctct                                                       >  6:67244/1‑148 (MQ=255)
                                                                                                     cagTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCGACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCt                                              >  8:263158/3‑149 (MQ=255)
                                                                                                      tGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTc                                             >  5:156118/1‑149 (MQ=255)
                                                                                                        tAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCg                                           >  1:268202/1‑149 (MQ=255)
                                                                                                               gtgGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGc                                    <  3:222617/149‑1 (MQ=255)
                                                                                                               gtgGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGc                                    >  3:249960/1‑149 (MQ=255)
                                                                                                                tgGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCATATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGcc                                   <  5:72614/149‑1 (MQ=255)
                                                                                                                tgGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCATATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGcc                                   <  5:67244/149‑1 (MQ=255)
                                                                                                                tgGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGcc                                   >  6:16302/1‑149 (MQ=255)
                                                                                                                tgGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGcc                                   <  6:150086/149‑1 (MQ=255)
                                                                                                                  gTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGccc                                  >  6:270097/1‑148 (MQ=255)
                                                                                                                      gTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAg                              <  2:72369/148‑1 (MQ=255)
                                                                                                                       tCCTGGGTCTCAAGGGCAAATGACCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTg                            >  3:353839/1‑149 (MQ=255)
                                                                                                                        ccTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTgg                           >  7:1988/1‑149 (MQ=255)
                                                                                                                                 cAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGc                  >  3:3224/1‑149 (MQ=255)
                                                                                                                                  aaGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCa                 >  3:127353/1‑149 (MQ=255)
                                                                                                                                    gggCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGTCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGc               <  7:314338/149‑1 (MQ=255)
                                                                                                                                     ggCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGACCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGcc              >  2:333474/1‑149 (MQ=255)
                                                                                                                                         aaTGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAAt          <  6:156118/149‑1 (MQ=255)
                                                                                                                                         aaTGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAAt          <  3:323880/149‑1 (MQ=255)
                                                                                                                                               ccAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTcc    >  7:70113/1‑149 (MQ=255)
                                                                                                                                                cAAATGATTTTGTTGCAAAAATTTGAGAGCCGACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAACTCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAACAGGCCCAGCGGCAGCCCAAAGTGTccc   >  2:129212/1‑149 (MQ=255)
                                                                                                                                                 aaaTGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATCTGTCCCg  >  3:135787/1‑149 (MQ=255)
                                                                                                                                                  aaTGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCg  >  7:6677/1‑148 (MQ=255)
                                                                                                                                                  aaTGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCg  >  7:12062/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCTGACCTGGCACATGGACTCCCTCACGGCTTACGTGCGGTACGGCTCCAATGAAACCACCTGGGCCATAGCGAGCCTCTTAGTGGGCTGGCCTGGCGTCCCTGTAGTGCTGTGGTTCGTCCTGGGTCTCAAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCG  >  NC_002947/2622059‑2622352

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