Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,796,855 C→G 7.9% intergenic (+39/‑93) PP_t43 → / → PP_2449 tRNA‑Val/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,796,8550CG 7.9% 203.4 / 7.4 63intergenic (+39/‑93)PP_t43/PP_2449tRNA‑Val/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/22);  new base G (2/3);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCTGTGAGAACG  >  NC_002947/2796707‑2797002
                                                                                                                                                    |                                                                                                                                                   
cGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGc                                                                                                                                                     <  5:364412/149‑1 (MQ=255)
  cATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAt                                                                                                                                                   <  3:313345/149‑1 (MQ=255)
    tCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAttt                                                                                                                                                 <  4:125645/149‑1 (MQ=255)
      gggCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTgg                                                                                                                                               <  4:226387/149‑2 (MQ=17)
           aCTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGccc                                                                                                                                           >  4:354778/1‑148 (MQ=255)
              aaTATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAag                                                                                                                                       >  3:126319/1‑149 (MQ=255)
                         cacaCGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc                                                                                                                            >  8:273528/1‑149 (MQ=255)
                         cacaCGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc                                                                                                                            >  7:35253/1‑149 (MQ=255)
                           cacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGt                                                                                                                          <  5:144166/149‑1 (MQ=255)
                           cacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGt                                                                                                                          <  8:35253/149‑1 (MQ=255)
                            acGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCATACCAAATCGAAAAAGCGCGATCAGAAATGATCGCCCTTATTGCATTTACGCCCGAACCAGCGTGTAACGta                                                                                                                         >  4:95426/1‑148 (MQ=18)
                             cgtcAGGCACATAGCTCAGTTGGTTAGATCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGttt                                                                                                                        <  8:90914/146‑1 (MQ=255)
                                   gCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTa                                                                                                                              <  5:332814/137‑1 (MQ=255)
                                   gCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTa                                                                                                                              >  6:332814/1‑137 (MQ=255)
                                       aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAt                                                                                                                                                   <  6:9115/112‑1 (MQ=255)
                                       aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAt                                                                                                                                                   >  5:9115/1‑112 (MQ=255)
                                              aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg                                                                                                                                                <  7:66082/108‑1 (MQ=14)
                                              aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg                                                                                                                                                <  7:71196/108‑1 (MQ=14)
                                              aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg                                                                                                                                                >  8:66082/1‑108 (MQ=14)
                                              aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg                                                                                                                                                >  8:71196/1‑108 (MQ=14)
                                                   gTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATAcccc                                                                                                  >  3:255305/1‑149 (MQ=255)
                                                       gagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGc                                                                                               <  3:95426/148‑1 (MQ=255)
                                                            ccaccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                                         <  7:77126/149‑1 (MQ=255)
                                                               ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGtt                                                                                                                         >  4:253248/1‑114 (MQ=255)
                                                               ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa                                                                                      >  5:131179/1‑149 (MQ=255)
                                                                 ttGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATc                                                                                    >  1:275250/1‑149 (MQ=255)
                                                                   gACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATcc                                                                                   <  4:126319/148‑1 (MQ=255)
                                                                        ggtggGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCt                                                                             >  2:156751/1‑149 (MQ=255)
                                                                         gtggGGGTCGTTGGTTCGAGTCCAATTGTGCATACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATAccc                                                                                                   >  7:201293/1‑126 (MQ=255)
                                                                               gTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGCAACGTTTGGTCGCGGTCCg                                                                                                            >  1:196694/1‑111 (MQ=255)
                                                                               gTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGCAACGTTTGGTCGCGGTCCg                                                                                                            <  2:196694/111‑1 (MQ=255)
                                                                                  gTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTgg                                                                   <  4:255305/149‑1 (MQ=255)
                                                                                   ttGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGa                                                                  <  3:354778/149‑1 (MQ=255)
                                                                                               cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTc                                                                           >  1:193155/1‑128 (MQ=255)
                                                                                               cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTc                                                                           <  2:193155/128‑1 (MQ=255)
                                                                                               cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTg                                                      >  7:22174/1‑149 (MQ=255)
                                                                                                      gCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAAc                                               >  5:122552/1‑149 (MQ=255)
                                                                                                         tACCAAATCGAAAAAGGGCGATCAGAAATGCTCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATACGCACCAAAGAAACAGGAGGAAAACGCGTGCGCAAAcaa                                             >  7:232701/1‑146 (MQ=255)
                                                                                                           ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCt                                          >  2:166985/1‑149 (MQ=255)
                                                                                                            cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtg                                         >  8:262263/1‑149 (MQ=255)
                                                                                                            cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtg                                         >  3:312315/1‑149 (MQ=255)
                                                                                                              aaTCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtgtg                                       <  7:273528/149‑1 (MQ=255)
                                                                                                                   aaaaaGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt                                  >  8:226564/1‑149 (MQ=255)
                                                                                                                   aaaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt                                  >  4:194556/1‑149 (MQ=255)
                                                                                                                    aaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATg                                 >  4:91086/1‑149 (MQ=255)
                                                                                                                     aaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg                                >  6:14162/1‑149 (MQ=255)
                                                                                                                      aaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg                                >  7:113307/1‑148 (MQ=255)
                                                                                                                      aaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg                                >  6:333147/1‑148 (MQ=255)
                                                                                                                        gggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAg                             >  3:64973/1‑149 (MQ=255)
                                                                                                                          gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg                                                                    <  6:181046/108‑1 (MQ=255)
                                                                                                                          gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg                                                                    >  5:181046/1‑108 (MQ=255)
                                                                                                                              tcatttctGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCg                       <  7:226564/142‑1 (MQ=255)
                                                                                                                               cAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGt                                                               >  4:20682/1‑108 (MQ=255)
                                                                                                                                 gAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGccc                    >  8:151169/1‑149 (MQ=255)
                                                                                                                                    aTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTgg                 <  3:91086/149‑1 (MQ=255)
                                                                                                                                          gCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCt           >  8:186428/1‑149 (MQ=255)
                                                                                                                                           cccTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTacac                                                           >  5:355892/1‑100 (MQ=255)
                                                                                                                                           cccTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTacac                                                           <  6:355892/100‑1 (MQ=255)
                                                                                                                                            ccTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCtgt         <  5:14162/149‑1 (MQ=255)
                                                                                                                                             cTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCGGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCtgtg        <  3:194556/149‑1 (MQ=255)
                                                                                                                                                tttGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTacac                                                           <  5:16506/95‑1 (MQ=255)
                                                                                                                                                tttGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTacac                                                           >  6:16506/1‑95 (MQ=255)
                                                                                                                                                tttGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCTGTgag      >  3:252959/1‑148 (MQ=255)
                                                                                                                                                tttGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCTGTgag      >  8:302273/1‑148 (MQ=255)
                                                                                                                                                    cATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCTGTGAGAAc   >  1:203558/1‑147 (MQ=255)
                                                                                                                                                    cATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCTGTGAGAACg  <  7:302273/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCTGGCGGGCTGTGAGAACG  >  NC_002947/2796707‑2797002

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: