Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,122,758 T→C 8.6% intergenic (+29/+118) PP_0983 → / ← PP_0984 membrane protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,122,7580TC 8.6% 197.1 / 5.6 69intergenic (+29/+118)PP_0983/PP_0984membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (3/3);  total (36/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 = 1.00e+00

GCCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAG  >  NC_002947/1122611‑1122905
                                                                                                                                                   |                                                                                                                                                    
gCCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTc                                                                                                                                                     <  5:114810/149‑1 (MQ=255)
 ccAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCa                                                                                                                                                    <  4:120487/149‑1 (MQ=255)
  cAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTACGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAg                                                                                                                                                   >  8:91754/1‑149 (MQ=255)
   aGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAgg                                                                                                                                                  <  7:36539/149‑1 (MQ=255)
     gACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGt                                                                                                                                                >  6:112553/1‑149 (MQ=255)
           cTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                          <  8:119864/148‑1 (MQ=255)
               ggCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtttct                                                                                                                                      <  8:290920/148‑3 (MQ=255)
                  cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCTGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCTTTCTtgtg                                                                                                                                  <  1:246754/149‑1 (MQ=255)
                         ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGtt                                                                                                                            >  2:113723/1‑148 (MQ=255)
                         ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGtt                                                                                                                            >  5:109395/1‑148 (MQ=255)
                         ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGtt                                                                                                                            >  4:134298/1‑148 (MQ=255)
                            ggTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATc                                                                                                                        >  8:124702/1‑149 (MQ=255)
                             gTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCt                                                                                                                       <  8:42755/149‑1 (MQ=255)
                             gTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCt                                                                                                                       <  3:134298/149‑1 (MQ=255)
                                  cGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGccc                                                                                                                  >  1:209050/1‑149 (MQ=255)
                                        cccGCCACTGCTGGCTGTGGTGATACCAGCGGGCATGTGTGCCTTGGCAGGTGTGTGGTTGGTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGTCCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCct                                                                                                             <  8:53680/148‑1 (MQ=255)
                                         ccGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCctct                                                                                                           <  4:162928/149‑1 (MQ=255)
                                         ccGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCctct                                                                                                           >  5:219191/1‑149 (MQ=255)
                                           gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                          >  4:262817/1‑116 (MQ=255)
                                           gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCt                                                                                                                                          <  3:262817/116‑1 (MQ=255)
                                           gCCACTGCTGGCGGTGGTGATACCAGCAGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTc                                                                                                         >  6:229363/1‑149 (MQ=255)
                                               ctgctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCgg                                                                                                     <  5:112553/149‑1 (MQ=255)
                                                        gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcg                                                                                            >  3:141321/1‑149 (MQ=255)
                                                             gATACCAGGGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtttcttgt                                                                                                                                   <  4:112874/105‑6 (MQ=38)
                                                             gATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTCGGCtttcttgt                                                                                                                                   >  3:112874/1‑100 (MQ=38)
                                                                aCCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGt                                                                                                                                   >  4:296610/1‑102 (MQ=255)
                                                                aCCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGt                                                                                                                                   <  3:296610/102‑1 (MQ=255)
                                                                   aGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAgg                                                                                  >  4:46621/1‑148 (MQ=255)
                                                                              gtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAgg                                                                       >  3:218576/1‑148 (MQ=255)
                                                                              gtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAgg                                                                       >  1:108962/1‑148 (MQ=255)
                                                                               tgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                         <  1:113723/145‑1 (MQ=255)
                                                                               tgCCTTGGCCGGTGTGTGGTTGTTACGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAg                                                                        <  7:91754/146‑1 (MQ=255)
                                                                                 ccTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATc                                                                   <  2:209050/149‑1 (MQ=255)
                                                                                   ttGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAg                                                                 >  3:232477/1‑149 (MQ=255)
                                                                                     ggCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAgg                                                                <  5:229363/148‑1 (MQ=255)
                                                                                      gCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                         <  8:20556/138‑1 (MQ=255)
                                                                                      gCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                         >  7:20556/1‑138 (MQ=255)
                                                                                      gCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGa                                                               <  7:124702/148‑1 (MQ=255)
                                                                                               tgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCg                                                      >  2:180795/1‑148 (MQ=255)
                                                                                                     ttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcgc                                                                                           <  1:138969/105‑1 (MQ=255)
                                                                                                     ttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGACCGGCCTCTTCGCGGGTAAGCCcgc                                                                                           >  2:138969/1‑105 (MQ=255)
                                                                                                        cgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACt                                             <  2:108962/148‑1 (MQ=255)
                                                                                                        cgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            >  5:198643/1‑149 (MQ=255)
                                                                                                           cGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCTTCAGGAAAGGGCCCGAAGGCAACTCCaa                                         <  1:180795/149‑1 (MQ=255)
                                                                                                             ggCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATc                                       <  4:218576/149‑1 (MQ=255)
                                                                                                                aGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCa                                                                  <  2:243271/119‑1 (MQ=255)
                                                                                                                aGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCa                                                                  >  1:243271/1‑119 (MQ=255)
                                                                                                                  ggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAgg                                                           <  6:205824/124‑1 (MQ=255)
                                                                                                                  ggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAgg                                                           >  5:205824/1‑124 (MQ=255)
                                                                                                                   gTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTgg                                 <  4:232477/149‑1 (MQ=255)
                                                                                                                      aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            <  3:291129/135‑1 (MQ=255)
                                                                                                                      aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            >  4:291129/1‑135 (MQ=255)
                                                                                                                        aGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAt                            >  1:63531/1‑149 (MQ=255)
                                                                                                                          ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTCAGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAtt                           >  1:202186/1‑149 (MQ=255)
                                                                                                                          ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAttt                          >  8:86178/1‑149 (MQ=255)
                                                                                                                          ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAtt                           >  7:204821/1‑148 (MQ=255)
                                                                                                                            tGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttct                        >  1:83825/1‑149 (MQ=255)
                                                                                                                              aGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc                      <  3:46621/149‑1 (MQ=255)
                                                                                                                                   cacgaaaaaGCCGACCCTGAGGTCGGCTTTCTTGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTATAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCtt                 >  6:160389/7‑149 (MQ=38)
                                                                                                                                   cacaAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGTAAAGGGCCCGAAGGCACCTCCAATCCTCTGGCTTATTTCTTCTGCtt                 >  2:226479/1‑149 (MQ=255)
                                                                                                                                      aaGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg              <  2:202186/149‑1 (MQ=255)
                                                                                                                                      aaGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg              >  3:157624/1‑149 (MQ=255)
                                                                                                                                        gAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCAGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCa            >  7:5886/1‑149 (MQ=255)
                                                                                                                                        aaaaGCCGACCCTGAGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCa            >  4:202621/2‑149 (MQ=255)
                                                                                                                                         aaaGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCAc           <  2:83825/149‑1 (MQ=255)
                                                                                                                                            gCCGGCCTCAGGGTCGGCTTTTTCGTGGGAGTTGATCTGGCCCGGCCTCTTCGCGCGTAACCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGCGCCCGAAGGCACCTCCAATCCTCTGGCGAATTTCTTATGCTTCGGGACCCg        >  7:120158/1‑149 (MQ=255)
                                                                                                                                                  cTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                       <  1:68360/128‑1 (MQ=255)
                                                                                                                                                  cTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                       >  2:68360/1‑128 (MQ=255)
                                                                                                                                                  cTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAg  <  8:204821/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTCGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAG  >  NC_002947/1122611‑1122905

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: