Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,075,511 A→T 5.7% I181N (ATC→AAC)  PP_2685 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,075,5110AT 5.7% 220.2 / 2.3 70I181N (ATC→AAC) PP_2685hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base T (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTCAAACCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGTGTGGCCTGGAAGAAAT  >  NC_002947/3075370‑3075655
                                                                                                                                             |                                                                                                                                                
ttCAAACCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGAtggtg                                                                                                                                             >  8:202728/1‑147 (MQ=255)
 tCAACCCGTGGCTCCACTGGCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGa                                                                                                                                           <  8:245780/148‑1 (MQ=255)
   aaaCCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTc                                                                                                                                        >  8:200041/1‑149 (MQ=255)
    aaCCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGAtggtg                                                                                                                                             >  6:222779/1‑143 (MQ=255)
    aaCCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGAtggtg                                                                                                                                             <  5:222779/143‑1 (MQ=255)
    aaCCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTc                                                                                                                                        >  4:335652/1‑148 (MQ=255)
    aaCCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTc                                                                                                                                        >  2:79011/1‑148 (MQ=255)
      ccGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTcgaa                                                                                                                                     <  1:79011/149‑1 (MQ=255)
       cGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTcgaac                                                                                                                                    <  6:112773/149‑1 (MQ=255)
            cTCCACTGCCTGCGAATCTGTTGATAATACGCATCGTCGCTGGTAATCCGATTGCGAGACATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAATACCAGCGGAATGCCCACCGACAGCTTGCTGCGCATGGTGGAGTCGAACGaaa                                                                                                                                >  4:113441/1‑148 (MQ=255)
             tCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAAt                                                                                                                               >  1:346938/1‑148 (MQ=255)
                 ctgcctgcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGc                                                                                                                          <  8:195993/149‑1 (MQ=255)
                 ctgcctgcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGc                                                                                                                          <  2:54489/149‑1 (MQ=255)
                 ctgcctgcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGACCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGc                                                                                                                          <  3:113441/149‑1 (MQ=255)
                      tgcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCcgcag                                                                                                                     >  7:148627/1‑149 (MQ=255)
                       gcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTTGAGTCGAACGAAATCAGGCcgcagc                                                                                                                    >  2:274376/1‑149 (MQ=255)
                       gcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCcgcagc                                                                                                                    >  3:221866/1‑149 (MQ=255)
                       gcGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCcgcagc                                                                                                                    >  4:327950/1‑149 (MQ=255)
                                ttGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCtt                                                                                                           >  5:54805/1‑149 (MQ=255)
                                  gATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtc                                                                                                          <  2:346938/148‑1 (MQ=255)
                                      ataGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtca                                                                                                         >  6:10120/1‑145 (MQ=255)
                                      ataGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtca                                                                                                         <  5:10120/145‑1 (MQ=255)
                                           cATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTGCATCACgcggc                                                                                                >  6:266932/1‑149 (MQ=255)
                                            aTCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtcat                                                                                                        <  2:369393/140‑1 (MQ=255)
                                            aTCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtcat                                                                                                        >  1:369393/1‑140 (MQ=255)
                                            aTCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAgcggcg                                                                                               >  6:211542/1‑149 (MQ=255)
                                                tcgCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCg                                                                                           <  7:202728/149‑1 (MQ=255)
                                                tcgCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCg                                                                                           <  4:221866/149‑1 (MQ=255)
                                                  gCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCgg                                                                                          <  1:344497/148‑1 (MQ=255)
                                                      gTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCcgcagcaac                                                                                                                 >  4:366373/1‑118 (MQ=255)
                                                      gTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCcgcagcaac                                                                                                                 <  3:366373/121‑4 (MQ=255)
                                                      gTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATc                                                                                     >  7:134355/1‑149 (MQ=255)
                                                            cgcTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGtt                                                                               >  1:338215/1‑149 (MQ=255)
                                                                 gcgcGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGTTGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtcat                                                                                                        >  3:16492/1‑119 (MQ=255)
                                                                 gcgcGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGTTGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtcat                                                                                                        <  4:11474/119‑1 (MQ=255)
                                                                 gcgcGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGTTGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTtcat                                                                                                        <  4:16492/119‑1 (MQ=255)
                                                                 gcgcGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGTTGGTGGAGTCGAACGAAATCAGGCCGCAGAGCAGCGCTtcat                                                                                                        >  3:11474/1‑119 (MQ=255)
                                                                  cgcgTCATGCTTTCGAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGt                                                                          <  6:54805/148‑1 (MQ=255)
                                                                  cgcgTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGt                                                                          <  3:335652/148‑1 (MQ=255)
                                                                  cgcgTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTc                                                                         <  8:134355/149‑1 (MQ=255)
                                                                       cATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGACCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTa                                                                    >  8:321549/1‑149 (MQ=255)
                                                                        aTGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTAt                                                                   >  6:20100/1‑149 (MQ=255)
                                                                         tGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTAtt                                                                  >  3:232260/1‑149 (MQ=255)
                                                                          gCTTTCCAGGCTGTTATTGCGGTACACCAGCTAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCt                                                                                                            <  3:157085/106‑1 (MQ=255)
                                                                          gCTTTCCAGGCTGTTATTGCGGTACACCAGCTAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCt                                                                                                            >  4:157085/1‑106 (MQ=255)
                                                                               ccAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCt                                                             <  7:200041/148‑1 (MQ=255)
                                                                                             cGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGGGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCttt                                               <  3:326683/148‑1 (MQ=255)
                                                                                              ggTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCa                                             <  5:20100/149‑1 (MQ=255)
                                                                                                tACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGACCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCAcc                                           <  7:321549/149‑1 (MQ=255)
                                                                                                    ccAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAg                                                                                  <  6:173328/106‑1 (MQ=255)
                                                                                                    ccAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAg                                                                                  >  5:173328/1‑106 (MQ=255)
                                                                                                      aGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGcc                                                              >  3:63359/1‑124 (MQ=255)
                                                                                                      aGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGcc                                                              <  4:63359/124‑1 (MQ=255)
                                                                                                        cAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTg                                     >  3:369211/1‑147 (MQ=255)
                                                                                                          aGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGa                                 <  3:327950/149‑1 (MQ=255)
                                                                                                              cAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTACATCACGCGGCGTGCGGTAATCCAGGTTGCGGGCGAGTATTGGCCTGCCGTACTTGCTTGCACCCAGCTGCAGTAACg                             >  8:113319/1‑149 (MQ=255)
                                                                                                                gCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGTCGGGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACggg                           <  5:266932/149‑1 (MQ=255)
                                                                                                                gCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACggg                           <  8:148627/149‑1 (MQ=255)
                                                                                                                gCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACggg                           >  2:244819/1‑149 (MQ=255)
                                                                                                                 cggatctcccACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACggg                           <  1:258468/141‑1 (MQ=255)
                                                                                                                       gCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTg                                                   <  3:108538/118‑1 (MQ=255)
                                                                                                                       gCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGGTGCGGTCGAGTATTGGCCTGCCGTACTTg                                                   >  4:108538/1‑118 (MQ=255)
                                                                                                                         ccACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACggg                           >  2:258468/1‑140 (MQ=255)
                                                                                                                         ccACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCGTTCCCCCAGCTGCAGAAACGGGGTGTCCGgt                  >  6:308746/1‑149 (MQ=255)
                                                                                                                                 aGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGTGTGGCCTg          <  4:232260/149‑1 (MQ=255)
                                                                                                                                  ggTTGCTGCGGATGGT‑GAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGTGTGGCCTgg         >  6:64506/1‑148 (MQ=255)
                                                                                                                                    ttGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGTGTGGCCTgg         >  4:195910/1‑147 (MQ=255)
                                                                                                                                         gcGGATGGTGGAGTCGAACGAAATCATGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGTGTGGCCTGGAAGAAAt  >  1:112475/1‑149 (MQ=255)
                                                                                                                                          cGGATGGTGGAGTCGAACGAAATCGGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGtgtg               <  7:168828/135‑1 (MQ=255)
                                                                                                                                          cGGATGGTGGAGTCGAACGAAATCGGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGtgtg               >  8:168828/1‑135 (MQ=255)
                                                                                                                                             aTGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCAccc                                          <  4:283221/105‑1 (MQ=255)
                                                                                                                                             aTGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCAccc                                          >  3:283221/1‑105 (MQ=255)
                                                                                                                                             |                                                                                                                                                
TTCAAACCGTCGCTCCACTGCCTGCGAATCTGTTGATAATAGGCATCGTCGCTGGTAATGCGCTTGCGCGTCATGCTTTCCAGGCTGTTATTGCGGTACACCAGCAAGTCCAGCGGCATGCCCACCGACAGGTTGCTGCGGATGGTGGAGTCGAACGAAATCAGGCCGCAGCGCAGCGCTTCATCAAGCGGCGTGCGGTAATCCAGGTTGCGGTCGAGTATTGGCCTGCCGTACTTGCTTTCACCCAGCTGCAGAAACGGGGTGTCCGGTGTGGCCTGGAAGAAAT  >  NC_002947/3075370‑3075655

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