Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,104,698 C→T 6.6% E62E (GAG→GAA)  PP_1879 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,104,6980CT 6.6% 174.7 / 2.0 60E62E (GAG→GAA) PP_1879hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

TAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCT  >  NC_002947/2104553‑2104846
                                                                                                                                                 |                                                                                                                                                       
tAAGCTTGTCGACCGGTACCAATATGGATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCAGGGTGCGAATTTGCCAACGGGTGGAACTc                                                                                                                                                       >  2:166591/1‑148 (MQ=255)
tAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTc                                                                                                                                                       >  7:208409/1‑148 (MQ=255)
  aGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAc                                                                                                                                                    <  7:96802/149‑1 (MQ=255)
   gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAcc                                                                                                                                                   >  4:339539/1‑149 (MQ=255)
   gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAcc                                                                                                                                                   >  2:120595/1‑149 (MQ=255)
       gTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGc                                                                                                                                                <  5:145116/148‑1 (MQ=255)
        tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATACCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAg                                                                                                                                              >  7:73945/1‑149 (MQ=255)
        tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGAAAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgtgg                                                                                                                                              <  6:90637/149‑4 (MQ=255)
           aCCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAgggg                                                                                                                                           >  7:25024/1‑149 (MQ=255)
            ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                          <  4:109653/149‑1 (MQ=255)
               gTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTACTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGaa                                                                                                                                       <  7:131974/149‑1 (MQ=255)
                   cAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAgg                                                                                                                                   <  3:339539/149‑1 (MQ=255)
                               ctctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑Ct                                                                                                                     <  7:223814/149‑1 (MQ=255)
                               ctctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑Ct                                                                                                                     >  3:255792/1‑149 (MQ=255)
                                tctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGc                                                                                                                               >  7:204254/1‑140 (MQ=255)
                                tctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGc                                                                                                                               <  8:204254/140‑1 (MQ=255)
                                  ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                          >  5:21477/1‑127 (MQ=255)
                                  ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                          <  6:21477/127‑1 (MQ=255)
                                      ggTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTACGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCgt                                                                                                             >  8:1449/1‑149 (MQ=255)
                                           gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTgg                                                                                                        >  6:321415/1‑149 (MQ=255)
                                           gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATACCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTg                                                                                                         <  8:73945/148‑1 (MQ=255)
                                            aTAGGCTTGCTGAATGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTggg                                                                                                       <  8:208409/149‑1 (MQ=255)
                                               ggCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTgggggt                                                                                                    <  4:146487/149‑4 (MQ=255)
                                                gCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTc                                                                                                   >  4:115500/1‑149 (MQ=255)
                                                 cTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTACTCATCTAATTACGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGCGGCAAGAGGTTGCCCGGAG‑‑AT‑CGCGCGAGTTGGGTTTcc                                                                                                  >  4:74731/1‑149 (MQ=255)
                                                     gtGAAGGTCTTCTTGTCGACCAGTAAGTCATTCGAGTCTAGGTATTGGTATTGCTCATCTATTTCCGGGTGCGAATTAGCACACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCGCGGAA‑‑CT‑CGAGCGTGTTGGGTTTcccccc                                                                                              <  5:90538/148‑1 (MQ=255)
                                                     cTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTcccccc                                                                                              <  3:115500/149‑1 (MQ=255)
                                                       gAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCtt                                                                                            >  8:190250/1‑149 (MQ=255)
                                                        aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTggg                                                                                                       <  5:65746/137‑1 (MQ=255)
                                                        aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTggg                                                                                                       >  6:65746/1‑137 (MQ=255)
                                                               tcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                                    >  1:415/1‑149 (MQ=255)
                                                                cttGTCGACCAGTACGTGATTCCATTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGTTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                                    <  1:166591/148‑1 (MQ=255)
                                                                cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATc                                                                                     <  5:321415/147‑1 (MQ=255)
                                                                 ttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTc                                                                                  >  4:365179/1‑149 (MQ=255)
                                                                    tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTccac                                                                               >  1:300907/1‑149 (MQ=255)
                                                                           gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCtt                                                                        >  2:192927/1‑149 (MQ=255)
                                                                               gTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGCCCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCCCGCTTGGTa                                                                    >  4:345920/1‑149 (MQ=255)
                                                                                  aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACgg                                                                 <  7:30343/149‑1 (MQ=255)
                                                                                      cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgg                                                             <  2:300907/149‑1 (MQ=255)
                                                                                         tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                          <  7:190250/149‑1 (MQ=255)
                                                                                         tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGGTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                          <  4:255792/149‑1 (MQ=255)
                                                                                                         agtCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGCTACGGTTGGGGGGACAAGAGAGGCCGAg                                          <  1:210481/147‑1 (MQ=255)
                                                                                                           tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGTAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACAGTGTTGGGTTTCCCCCCTTCAACATCGTCGACAACGCTTGGTACGGTTGGGGGGACAAGAGAGCCCGAGTc                                        >  4:9535/1‑149 (MQ=255)
                                                                                                             aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCa                                      >  7:234909/1‑149 (MQ=255)
                                                                                                              tCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAg                                     >  2:203771/1‑149 (MQ=255)
                                                                                                                   ttCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCATGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                                >  7:182691/1‑149 (MQ=255)
                                                                                                                    tCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                                <  2:415/148‑1 (MQ=255)
                                                                                                                      cGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                              <  1:192927/148‑1 (MQ=255)
                                                                                                                         gTGCGAATTTGCCCACGGGTGGAACTCCACCCGCACGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                          <  8:182691/149‑1 (MQ=255)
                                                                                                                               aTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATATGGTGCCCGGAATCCTAAGACCTCATACAAATACCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAg                       <  1:120595/149‑1 (MQ=255)
                                                                                                                                  cccctgcGGGTGGAGTTCCACCCGTAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                                >  2:208803/7‑134 (MQ=255)
                                                                                                                                  cccctgcGGGTGGAGTTCCACCCGTAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                                <  1:208803/128‑1 (MQ=255)
                                                                                                                                     ctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt              >  7:299969/4‑149 (MQ=255)
                                                                                                                                       aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTa            >  7:90083/1‑149 (MQ=255)
                                                                                                                                         gggTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                              <  7:115762/129‑1 (MQ=255)
                                                                                                                                         gggTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                              >  8:115762/1‑129 (MQ=255)
                                                                                                                                            tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCg         >  6:296039/1‑147 (MQ=255)
                                                                                                                                              gAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgc     >  1:252649/1‑149 (MQ=255)
                                                                                                                                              gAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgc     <  2:252649/149‑1 (MQ=255)
                                                                                                                                              gAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgc     >  3:358369/1‑149 (MQ=255)
                                                                                                                                                 cTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCt  >  3:90998/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
TAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCGCT  >  NC_002947/2104553‑2104846

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: ‑