Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,766,621 T→C 6.2% E1417E (GAA→GAG)  PP_5060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,766,6210TC 6.2% 187.2 / 2.1 64E1417E (GAA→GAG) PP_5060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AGAAAATTCAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGCCAGCATGGCCTCGGCCAGCTCCCGGTT  >  NC_002947/5766478‑5766764
                                                                                                                                               |                                                                                                                                               
agaAAATTCAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGCTCACg                                                                                                                                            >  1:141334/1‑149 (MQ=255)
agaAAATTCAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGCTCACg                                                                                                                                            >  1:136498/1‑149 (MQ=255)
 gaAAATTCAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACg                                                                                                                                            >  6:35422/1‑148 (MQ=255)
     aTTCAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACt                                                                                                                                       >  5:102966/1‑149 (MQ=255)
        cAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCg                                                                                                                                     <  6:328269/148‑1 (MQ=255)
         aGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGa                                                                                                                                   >  7:235856/1‑149 (MQ=255)
            ccTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGTCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCTTGCTCTGCGTTGGCATCGCTTAGCTGTTCACGCAACTCGGaca                                                                                                                                 <  5:46712/148‑1 (MQ=255)
             ctctGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACt                                                                                                                               >  2:286906/1‑149 (MQ=255)
              tctGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACtt                                                                                                                              >  2:173510/1‑149 (MQ=255)
               ctGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTc                                                                                                                             >  6:185598/1‑149 (MQ=255)
                  gCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGtt                                                                                                                          <  5:78878/149‑1 (MQ=255)
                  gCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGtt                                                                                                                          <  2:49824/149‑1 (MQ=255)
                  gCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGt                                                                                                                           <  1:211848/148‑1 (MQ=255)
                  gCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGt                                                                                                                           <  8:75849/148‑1 (MQ=255)
                  gCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGCTCACGCAACTCGGACACTTCGtt                                                                                                                          <  2:136498/149‑1 (MQ=255)
                  gCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGCTCACGCAACTCGGACACTTCGtt                                                                                                                          <  2:141334/149‑1 (MQ=255)
                     ccTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGAt                                                                                                                       >  8:248391/1‑149 (MQ=255)
                           tAGCGCGTCTTGCTGTCAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGcc                                                                                                                 <  5:109716/149‑1 (MQ=255)
                           tAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGcc                                                                                                                 >  7:232742/1‑149 (MQ=255)
                           tAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGcc                                                                                                                 <  5:166062/149‑1 (MQ=255)
                              cgcgTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCg                                                                                                                <  7:8825/147‑1 (MQ=255)
                              cgcgTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCg                                                                                                                >  8:8825/1‑147 (MQ=255)
                               gcgTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGcc                                                                                                                 <  4:274684/145‑1 (MQ=255)
                               gcgTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGcc                                                                                                                 >  3:274684/1‑145 (MQ=255)
                                cgTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCa                                                                                                            <  1:173510/149‑1 (MQ=255)
                                cgTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCa                                                                                                            <  5:185598/149‑1 (MQ=255)
                                   cTTGCTGTAAAGCGCTACCAAACACTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCACCAGCGCCTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGGAACTCGGACACTTCGTTGATCTGCCCGTCCAGc                                                                                                          >  3:232584/1‑148 (MQ=255)
                                    ttGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCAt                                                                                                        >  3:240060/1‑149 (MQ=255)
                                    ttGCAGTAAAGCGCTACCAAACAGTCGCGGCGGGCCAGACGCAACTGCATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTACGTTGGCATCGCTGAGCTGTTAACGCAACTCGGACACTTCGTTTATATGTCAGTCCAGCAt                                                                                                        >  6:182158/1‑149 (MQ=255)
                                      gCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCAGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCg                                                                                                      >  7:281936/1‑149 (MQ=255)
                                          tAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACAt                                                                                                  >  7:147942/1‑149 (MQ=255)
                                          tAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCAGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACAt                                                                                                  >  3:158728/1‑149 (MQ=255)
                                             aGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAg                                                                                               >  5:253406/1‑149 (MQ=255)
                                                    ccAAACAGTCGCGGCGGGCCTGACGCAACTGAATCTGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCata                                                                                        >  8:92576/1‑149 (MQ=255)
                                                         cAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTc                                                                                   <  6:102966/149‑1 (MQ=255)
                                                         cAGTCGCGGCGGGCCTGACGAAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCTTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCTTTGTTCTGGACGTCCAGCATCGACATCAGGCCCATATCGTc                                                                                   <  5:182158/149‑1 (MQ=255)
                                                          aGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCa                                                                                  <  1:171907/149‑1 (MQ=255)
                                                                  cggGGCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGGCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATGTGGGCGTCCAGCATCGAAATCAGGCCCATATCGTCAGCATCCtt                                                                          <  4:158728/149‑1 (MQ=255)
                                                                  cggGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCtt                                                                          <  2:185265/149‑1 (MQ=255)
                                                                   gggCCTGACGCAACTGAATCTGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTg                                                                         <  7:92576/149‑1 (MQ=255)
                                                                    ggCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATcc                                                                            <  5:70214/145‑1 (MQ=255)
                                                                    ggCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATcc                                                                            >  6:70214/1‑145 (MQ=255)
                                                                       cTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCa                                                                                                <  1:150465/122‑1 (MQ=255)
                                                                       cTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCa                                                                                                >  2:150465/1‑122 (MQ=255)
                                                                          aCGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAAc                                                                    <  5:35422/147‑1 (MQ=255)
                                                                            gCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCCTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAgg                                                                >  4:214283/1‑149 (MQ=255)
                                                                              aaCTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGc                                                              >  6:193280/1‑149 (MQ=255)
                                                                                 tGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGctct                                                           <  8:235856/149‑1 (MQ=255)
                                                                                   aaTCAGCGCATCGTCCATGCCCGCCAGCAGGGCGTTCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTgg                                                         <  4:232584/149‑1 (MQ=255)
                                                                                        gcgcATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGa                                                    <  3:95777/149‑1 (MQ=255)
                                                                                        gcgcATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGa                                                    <  8:281936/149‑1 (MQ=255)
                                                                                           cATCGTCAATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGc                                                  >  1:222904/1‑148 (MQ=255)
                                                                                                tCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGa                                             >  1:171073/1‑148 (MQ=255)
                                                                                                 ccATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACa                                           >  2:290896/1‑149 (MQ=255)
                                                                                                   aTGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCCTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACAtt                                         <  3:214283/149‑1 (MQ=255)
                                                                                                    tGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACAttt                                        <  7:248391/149‑1 (MQ=255)
                                                                                                       ccGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCg                                     >  1:172759/1‑149 (MQ=255)
                                                                                                         gCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGtt                                   <  4:240060/149‑1 (MQ=255)
                                                                                                               agcGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGc                             <  1:286906/149‑1 (MQ=255)
                                                                                                               agcGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGc                             <  2:171073/149‑1 (MQ=255)
                                                                                                                             gTTGGCATCGCTGAGCTGTTCACTCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACGAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGCCAGCATGGCCTCgg               <  6:253406/149‑1 (MQ=255)
                                                                                                                                   aTCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGCCAGCATGGCCTCGGCCAGCt         <  2:279647/149‑1 (MQ=255)
                                                                                                                                          aGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGTAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGCCAGCATGGCCTCGGCCAGCTCCCGgtt  >  3:125328/1‑149 (MQ=255)
                                                                                                                                          aGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGCCAGCATGGCCTCGGCCAGCTCCCGgtt  >  3:223890/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
AGAAAATTCAGCCCTCTGGCACCTGGATAGCGCGTCTTGCTGTAAAGCGCTACCAAACAGTCGCGGCGGGCCTGACGCAACTGAATCAGCGCATCGTCCATGCCCGCCAGCAGCGCGTGCTCTGCGTTGGCATCGCTGAGCTGTTCACGCAACTCGGACACTTCGTTGATCTGGCCGTCCAGCATCGACATCAGGCCCATATCGTCAGCATCCTTGGCAACCAGGGCTCTGGCTTGAGCGATGACATTTTCGTTCTGTGCCAGCATGGCCTCGGCCAGCTCCCGGTT  >  NC_002947/5766478‑5766764

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