Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 554,939 T→A 7.0% R51R (CGT→CGA rplP → 50S ribosomal protein L16

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947554,9390TA 7.0% 197.7 / 5.0 71R51R (CGT→CGArplP50S ribosomal protein L16
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base A (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TGCAACCAAAGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTCAGC  >  NC_002947/554791‑555081
                                                                                                                                                     |                                                                                                                                               
tGCAACCAAAGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCAcgt                                                                                                                                                 >  6:149863/1‑149 (MQ=255)
      cAAAGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGTAGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCGGAAAGATGTTGATCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTgc                                                                                                                                            >  1:83721/1‑149 (MQ=255)
      cAAAGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTgcg                                                                                                                                           >  5:161824/1‑149 (MQ=255)
         aGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGACGTgcggac                                                                                                                                        <  6:50943/149‑4 (MQ=255)
         aGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGACGTgcggac                                                                                                                                        <  8:100126/149‑4 (MQ=255)
             tACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGAccc                                                                                                                                    >  5:43028/1‑149 (MQ=255)
               cAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGt                                                                                                                                  <  3:240537/149‑1 (MQ=255)
                aaaaTTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTc                                                                                                                                 >  2:228016/1‑149 (MQ=255)
                  aaTTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGATTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCACATCGAGTCCGAACGTCGTGCGCTGACCCGTccc                                                                                                                               >  2:255101/1‑147 (MQ=255)
                       cGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGtt                                                                                                                            >  8:89347/1‑147 (MQ=255)
                         cAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGc                                                                                                                        <  1:228016/149‑1 (MQ=255)
                          aaGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTgc                                                                                                                                            >  1:192008/1‑128 (MQ=255)
                               gATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGACGTGCGCTGACCCGTCACGTTAAGCGTGGCg                                                                                                                 >  5:234744/1‑149 (MQ=255)
                                   aCTGGCCACAACCGTTGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTaa                                                                                                             <  3:171315/149‑1 (MQ=255)
                                   aCTGGCCACAACCGTGGTCTGGCACTGCGCGTTAGCAAAGTCTGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTaa                                                                                                             <  1:255101/149‑1 (MQ=255)
                                   aCTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGt                                                                                                               >  2:233979/1‑147 (MQ=255)
                                   aCTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTaa                                                                                                             <  1:233979/149‑1 (MQ=255)
                                   aCTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTaa                                                                                                             <  6:43028/149‑1 (MQ=255)
                                     tGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGGCACGGTAAGCGTGGCGGTAAg                                                                                                            >  3:167845/1‑148 (MQ=255)
                                        ccACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCt                                                                                                        <  6:125962/149‑1 (MQ=255)
                                              ccGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATcc                                                                                                  >  3:296515/1‑149 (MQ=255)
                                               cGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGATGCTCGCGCTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCg                                                                                                 >  5:218537/1‑149 (MQ=255)
                                                gTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATccat                                                                                                >  3:137932/1‑147 (MQ=255)
                                                   gTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCgtgt                                                                                              <  5:130127/148‑1 (MQ=255)
                                                    tCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCgt                                                                                                <  3:297602/145‑1 (MQ=255)
                                                    tCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCgt                                                                                                >  4:297602/1‑145 (MQ=255)
                                                        gCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTccc                                                                                        >  8:43089/1‑149 (MQ=255)
                                                        gCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTccc                                                                                        <  6:132955/149‑1 (MQ=255)
                                                                 cgGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGc                                                                                <  1:308906/149‑1 (MQ=255)
                                                                         aaGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTg                                                                                                       >  5:87647/1‑118 (MQ=255)
                                                                         aaGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTg                                                                                                       <  6:87647/118‑1 (MQ=255)
                                                                         aaGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGAtctc                                                                         >  3:260489/1‑148 (MQ=255)
                                                                         aaGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTcc                                                                        <  3:37714/149‑1 (MQ=255)
                                                                         aaGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTcc                                                                        <  5:149863/149‑1 (MQ=255)
                                                                              aGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCaaaa                                                                    <  8:58448/148‑1 (MQ=255)
                                                                              aGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCaaaa                                                                    >  4:264879/1‑148 (MQ=255)
                                                                                cTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGACGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAAc                                                                 <  6:234744/149‑1 (MQ=255)
                                                                                 ttCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAAcc                                                                >  7:161565/1‑149 (MQ=255)
                                                                                     gCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCtt                                                            >  6:110104/1‑149 (MQ=255)
                                                                                       gAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTGCGCTCGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGa                                                          <  2:305039/149‑1 (MQ=255)
                                                                                        aaTTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCATATCGAGTCCGCACGGCGTGCGCTGACCCGTCACGTTAAGCGTGGCCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGa                                                          <  1:238757/149‑1 (MQ=255)
                                                                                            cGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGtt                                                      <  7:89347/148‑1 (MQ=255)
                                                                                             gaccTGAAACCAGTTGCTCGAGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCCTGCGCTGACCCGTCACGTTAACCGTGGCGGAAAGATCTGGATCCGTGTGTTCCCGGACAAGCAGAACTCCAAAAAACCTCTTGAGGTTCg                                                    >  6:99569/3‑149 (MQ=255)
                                                                                             gCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCaaaa                                                                    >  4:187720/1‑133 (MQ=255)
                                                                                             gCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCaaaa                                                                    <  3:187720/133‑1 (MQ=255)
                                                                                               ccTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTa                                                  >  2:58977/1‑149 (MQ=255)
                                                                                                cTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTTAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTGCGTAt                                                 >  8:18105/1‑149 (MQ=255)
                                                                                                cTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTAt                                                 >  1:166460/1‑149 (MQ=255)
                                                                                                   aaaGCTGTTGCTCGCGGTCGCCTCACTGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGGTAAGGGTGGCGGTAAGATATGGATCCGTGTGGTCCCGGACAAGCGGATGTCAAAAACGCGTCTTGAGGTTCGTCTggg                                              >  8:222688/1‑149 (MQ=255)
                                                                                                   aaaGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATgg                                               <  1:11544/148‑1 (MQ=255)
                                                                                                   aaaGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATgg                                               <  5:110104/148‑1 (MQ=255)
                                                                                                   aaaGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATgg                                               <  7:43089/148‑1 (MQ=255)
                                                                                                   aaaGCTGTTGCTCGCGGTCGCCGCACCGCTCGCCAGATCGAGTCCGCACGTCGGGCGCTGACCCGTCACGTTAAGCGGGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATgg                                               <  6:218537/148‑1 (MQ=255)
                                                                                                        tGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAg                                         >  4:276277/1‑149 (MQ=255)
                                                                                                          ttGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGt                                       <  3:264879/149‑1 (MQ=255)
                                                                                                            gCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGt                                       >  5:181521/1‑147 (MQ=255)
                                                                                                             cTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAg                                    <  6:161824/149‑1 (MQ=255)
                                                                                                                gcgGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAggg                                  <  3:276277/148‑1 (MQ=255)
                                                                                                                 cgGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCt                                >  8:235227/1‑149 (MQ=255)
                                                                                                                  ggTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTc                               >  3:134212/1‑149 (MQ=255)
                                                                                                                      gCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTg                           >  1:176924/1‑149 (MQ=255)
                                                                                                                       ccTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTgg                          <  1:58977/149‑1 (MQ=255)
                                                                                                                             cGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACt                    <  2:83721/149‑1 (MQ=255)
                                                                                                                                     aGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGcc            >  4:216377/1‑149 (MQ=255)
                                                                                                                                     aGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGcc            >  1:57504/1‑149 (MQ=255)
                                                                                                                                          gAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGAt       >  7:13653/1‑149 (MQ=255)
                                                                                                                                              gcGCACGACGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTCAg   >  6:147264/2‑149 (MQ=255)
                                                                                                                                               cGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTCAGc  >  1:46495/1‑149 (MQ=255)
                                                                                                                                               cGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTGGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTc     >  1:273937/1‑146 (MQ=255)
                                                                                                                                               cGCACGTCGTGCGCTGACCCGTCACGTTAACCGTGGCGGTAATATCTGGTTCCGTGTTTTTCCGGAAAATCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTc     <  2:273937/146‑1 (MQ=255)
                                                                                                                                                gCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTCAGc  >  6:111560/1‑148 (MQ=255)
                                                                                                                                                 cACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTCAGc  >  6:130028/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TGCAACCAAAGCGTACAAAATTCCGCAAGCAGATGACTGGCCACAACCGTGGTCTGGCACTGCGCGGTAGCAAAGTCAGCTTCGGCGAATTCGCCCTGAAAGCTGTTGCTCGCGGTCGCCTCACCGCTCGCCAGATCGAGTCCGCACGTCGTGCGCTGACCCGTCACGTTAAGCGTGGCGGTAAGATCTGGATCCGTGTGTTCCCGGACAAGCCGATCTCCAAAAAACCTCTTGAGGTTCGTATGGGTAAAGGTAAGGGCTCCGTGGAATACTGGGTTGCCCAGATTCAGC  >  NC_002947/554791‑555081

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: