Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 850,994 C→G 5.6% intergenic (‑574/‑231) pcs ← / → hemA phosphatidylcholine synthase/glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947850,9940CG 5.6% 287.1 / 7.0 90intergenic (‑574/‑231)pcs/hemAphosphatidylcholine synthase/glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base G (3/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

TCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAT  >  NC_002947/850850‑851142
                                                                                                                                                |                                                                                                                                                    
tctGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGc                                                                                                                                                 >  1:438779/1‑150 (MQ=255)
     cTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                            <  2:2344738/150‑1 (MQ=255)
              ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGc                                                                                                                                                 <  2:573294/136‑1 (MQ=255)
              ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGc                                                                                                                                                 >  1:573294/1‑136 (MQ=255)
              ttACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTg                                                                                                                                   >  2:1362322/1‑150 (MQ=255)
                  tGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCa                                                                                                                               >  2:2243612/1‑150 (MQ=255)
                  tGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCa                                                                                                                               >  2:111816/1‑150 (MQ=255)
                   gTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCaa                                                                                                                              >  1:714866/1‑150 (MQ=255)
                    tCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAg                                                                                                                             >  1:1185538/1‑150 (MQ=255)
                      gCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                           >  2:1383271/1‑150 (MQ=255)
                          cTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGc                                                                                                                       <  2:2028144/150‑1 (MQ=255)
                              gcgcgcGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc                                                                                                                   <  1:2207396/150‑1 (MQ=255)
                                 tgatcttGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGc                                                                                                                >  2:2390423/6‑150 (MQ=255)
                                    ggTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                          <  2:2186632/121‑2 (MQ=255)
                                    ggTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                                          >  1:2186632/1‑120 (MQ=255)
                                     gTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctcc                                                                                                            >  1:2087738/1‑150 (MQ=255)
                                      ttGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCCTCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGATCGGcctcct                                                                                                           >  1:975297/1‑150 (MQ=255)
                                          tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATGGATGCCCGt                                                                                                                                             >  1:1303049/1‑112 (MQ=255)
                                          tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATGGATGCCCGt                                                                                                                                             <  2:1303049/112‑1 (MQ=255)
                                          tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                       >  2:1336904/1‑150 (MQ=255)
                                          tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                       >  2:124098/1‑150 (MQ=255)
                                              gcgcTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTAACGGTCATCCCTGCCCGTCACATCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                   <  1:1261041/150‑1 (MQ=255)
                                                 cTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCGCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                <  1:1968809/150‑1 (MQ=255)
                                                  tGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaag                                                                                               >  1:202121/1‑150 (MQ=255)
                                                            gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctc                                                                                                             <  2:2309523/126‑1 (MQ=255)
                                                            gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctc                                                                                                             >  1:2309523/1‑126 (MQ=255)
                                                            gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAagatcaagatca                                                                                     <  2:1051381/150‑1 (MQ=255)
                                                            gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAagatcaagatca                                                                                     <  1:2390423/150‑1 (MQ=255)
                                                            gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAagatcaagatca                                                                                     <  2:2087738/150‑1 (MQ=255)
                                                              cTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAgg                                                                                   <  2:438779/150‑1 (MQ=255)
                                                               ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                           >  1:243057/1‑109 (MQ=255)
                                                               ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc                                                                                                                           <  2:243057/109‑1 (MQ=255)
                                                               ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGa                                                                                  <  1:380870/150‑1 (MQ=255)
                                                                tCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGAt                                                                                 <  2:714866/150‑1 (MQ=255)
                                                                         cGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa                                                                        >  1:1671851/1‑150 (MQ=255)
                                                                          gAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAg                                                                       <  1:124098/150‑1 (MQ=255)
                                                                          gAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGCACGGTTTGCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAg                                                                       >  1:760179/1‑150 (MQ=21)
                                                                           aGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGa                                                                      <  2:567067/150‑1 (MQ=255)
                                                                           aGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGa                                                                      <  1:2243612/150‑1 (MQ=255)
                                                                            gCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGAt                                                                     <  1:1336904/150‑1 (MQ=255)
                                                                                aGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCg                                                                                                                  >  1:127222/1‑101 (MQ=255)
                                                                                aGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCg                                                                                                                  <  2:127222/101‑1 (MQ=255)
                                                                                aGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaaa                                                                 >  1:1862737/1‑150 (MQ=255)
                                                                                  ggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAg                                                               >  2:43683/1‑150 (MQ=255)
                                                                                    cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                             >  1:236563/1‑150 (MQ=255)
                                                                                       gCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGc                                                          <  2:236563/150‑1 (MQ=255)
                                                                                        cTCCCGGGAGAACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCa                                                         >  1:2133605/1‑150 (MQ=255)
                                                                                         tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATTAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                        >  1:508949/1‑150 (MQ=255)
                                                                                           ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGGCCTGc                                                                                                                       <  2:2394188/85‑1 (MQ=255)
                                                                                           ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGc                                                                                                                       >  1:2394188/1‑85 (MQ=255)
                                                                                           ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                      >  2:52954/1‑150 (MQ=255)
                                                                                           ccGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                      >  2:2155918/1‑150 (MQ=255)
                                                                                           ccGGCAGCACCGCTCTGCAGGCCCAATCCTCGGGCACAGCTCTTGACGGGCATCCATGCCCATCACCTCCCCGCGCAAGGCCTGCGTTCGACCTCCTGAAGTGGCGAAGATCAAGATGAGGATCAAAAGCCAGATCAAAAGCCGGCAAGa                                                      <  2:1570163/150‑1 (MQ=21)
                                                                                              ggaggaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCt                                                   <  2:1862737/150‑1 (MQ=255)
                                                                                                     cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATNCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatca                                                                                           >  1:810264/1‑103 (MQ=255)
                                                                                                     cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatca                                                                                           <  2:810264/103‑1 (MQ=255)
                                                                                                      gcgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                           <  2:1185538/150‑1 (MQ=255)
                                                                                                       cgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACAcc                                          <  1:52954/150‑1 (MQ=255)
                                                                                                        gcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCg                                         <  2:202121/150‑1 (MQ=255)
                                                                                                          tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa                                              <  2:290934/143‑1 (MQ=255)
                                                                                                          tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa                                              >  1:290934/1‑143 (MQ=255)
                                                                                                          tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCt                                       >  1:1338843/1‑150 (MQ=255)
                                                                                                                ccccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                      >  2:546858/1‑97 (MQ=255)
                                                                                                                ccccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                                                                                      <  1:546858/97‑1 (MQ=255)
                                                                                                                ccccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCa                                 <  1:2155918/150‑1 (MQ=255)
                                                                                                                  ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTa                                      >  1:409599/1‑143 (MQ=255)
                                                                                                                  ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTa                                      <  2:409599/143‑1 (MQ=255)
                                                                                                                   cATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatc                                                                                            <  2:2335528/88‑1 (MQ=255)
                                                                                                                   cATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatc                                                                                            >  1:2335528/1‑88 (MQ=255)
                                                                                                                     tCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaaga                                                                                        >  1:1558186/1‑90 (MQ=255)
                                                                                                                     tCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaaga                                                                                        <  2:1558186/90‑1 (MQ=255)
                                                                                                                       cTCGGGCACAGCGCTTGACGGGCATCCATGCGCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGa                                                                      >  2:1486681/1‑106 (MQ=255)
                                                                                                                       cTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGa                                                                      <  1:1486681/106‑1 (MQ=255)
                                                                                                                       cTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACg                          >  1:1953176/1‑150 (MQ=255)
                                                                                                                        tCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCgg                                                           >  2:1530780/1‑116 (MQ=255)
                                                                                                                        tCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGAGCAAAAGCCAGATCAAAAGCCgg                                                           <  1:1530780/116‑1 (MQ=255)
                                                                                                                             cacaGCGCTTGACGGGCATCCATGCCCCTCACCTCCCTGCTCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg                    <  1:1362322/150‑1 (MQ=255)
                                                                                                                              acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                   >  2:1158973/1‑70 (MQ=255)
                                                                                                                              acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                   <  1:1158973/70‑1 (MQ=255)
                                                                                                                              acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGc                   >  1:1668795/1‑150 (MQ=255)
                                                                                                                                aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                 <  1:1383271/150‑1 (MQ=255)
                                                                                                                                   gcTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAg              >  2:1489268/1‑150 (MQ=255)
                                                                                                                                      tGACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCg           >  2:651354/1‑150 (MQ=255)
                                                                                                                                          gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                           <  2:2396335/114‑1 (MQ=255)
                                                                                                                                          gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                           >  1:2396335/1‑114 (MQ=255)
                                                                                                                                          gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc       >  2:678511/1‑150 (MQ=255)
                                                                                                                                            gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGACCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcgc     >  1:1588822/1‑150 (MQ=255)
                                                                                                                                               tCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc       <  2:400109/145‑1 (MQ=255)
                                                                                                                                               tCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc       >  1:400109/1‑145 (MQ=255)
                                                                                                                                               tCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAt  <  2:1671851/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TCTGGCTTTTGCTCTTACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAT  >  NC_002947/850850‑851142

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