Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,759,133 A→T 6.7% Q99L (CAA→CTA)  PP_1565 → HK97 family portal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,759,1330AT 6.7% 181.9 / 2.6 60Q99L (CAA→CTA) PP_1565HK97 family portal protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base T (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCTCGACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAACATGTCGCCGCCGCGGCGCCT  >  NC_002947/1758986‑1759274
                                                                                                                                                   |                                                                                                                                             
tctcGACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGTCGAGCCGGTGGTGGCCAGTTATCACCCGGTCAACTTCGTCATAAGCCAGCAGCCGAAAGCCGATCa                                                                                                                                               >  5:343242/1‑148 (MQ=255)
 ctcGACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCaaa                                                                                                                                             >  3:35870/1‑149 (MQ=255)
  tcGACGGTGTGGGCTTGCGTGCGCTTGAAAGCAGAGACCATAGCAACGATTCCGATAGCAAAGTATGAAGAAAAGAACGGCGCGCCGGTGGTGGCCAGTTCTAACCCGGTAAACGTCGTAATAAGCCAGGAGCCGAACGGCGATCAAAc                                                                                                                                            >  7:272941/1‑149 (MQ=255)
    gACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCaaa                                                                                                                                             <  5:361989/146‑1 (MQ=255)
    gACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAAccc                                                                                                                                          <  4:35870/149‑1 (MQ=255)
    gACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGACCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAAccc                                                                                                                                          >  6:272129/1‑149 (MQ=255)
       ggcgggggCTTGCGGGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGAGCGCTCTGTATGAGGACAAGTACGGCGCGCCTGGGGTGGCCAGTTCTCACGCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCgg                                                                                                                                       <  4:419769/144‑1 (MQ=255)
                ttGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCTAACCCCGGTGGAGTTCt                                                                                                                              >  5:381424/1‑149 (MQ=255)
                ttGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCAACAGCCAGCAGCCGAACGCCGATCTAACCCCGGTGGAGTTCt                                                                                                                              >  5:374527/1‑149 (MQ=255)
                  gcgtgcgCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTg                                                                                                                             <  7:378996/148‑1 (MQ=255)
                   cgtgcgCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTgg                                                                                                                            <  2:15382/148‑1 (MQ=255)
                       cgcTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGaa                                                                                                                       <  7:29895/149‑1 (MQ=255)
                         cTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACg                                                                                                                     >  4:133900/1‑149 (MQ=255)
                            gATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTAGGAGGCCAAGAGCGGCGCGCCCGTGGTGGCCAGTTCTCCCCCGGTCAACTGCGTCATCCGCCAGCAGCCGAGGGCGGGTCAAACCCCGGTGGAGTTCTGGGAGTACGTCa                                                                                                                  >  6:363197/1‑149 (MQ=255)
                              tCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCaaa                                                                                                                                             >  3:248352/1‑120 (MQ=255)
                              tCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCaaa                                                                                                                                             <  4:248352/120‑1 (MQ=255)
                              tCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATg                                                                                                                <  6:324879/149‑1 (MQ=255)
                                 ccGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGc                                                                                                              >  4:54942/1‑148 (MQ=255)
                                   gagaCCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAg                                                                                                           >  7:259470/1‑149 (MQ=255)
                                        cATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCctgc                                                                                                      <  8:238226/149‑1 (MQ=255)
                                             cAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTc                                                                                                 <  7:168879/149‑1 (MQ=255)
                                             cAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTc                                                                                                 >  5:252969/1‑149 (MQ=255)
                                                 gCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAgg                                                                                             <  3:54942/149‑1 (MQ=255)
                                                   tGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAggg                                                                                            >  8:259909/1‑148 (MQ=255)
                                                      cGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAAc                                                                                        <  3:133900/149‑1 (MQ=255)
                                                            ctctGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAg                                                                                              <  8:219513/137‑1 (MQ=255)
                                                            ctctGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAg                                                                                              >  7:219513/1‑137 (MQ=255)
                                                                gTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGc                                                                               >  7:11359/1‑148 (MQ=255)
                                                                  aTGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAg                                                                            >  4:237778/1‑149 (MQ=255)
                                                                  aTGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAg                                                                            >  3:361725/1‑149 (MQ=255)
                                                                   tGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGc                                                                           <  8:259470/149‑1 (MQ=255)
                                                                   tGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGc                                                                           >  2:17738/1‑149 (MQ=255)
                                                                   tGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGc                                                                           >  6:307374/1‑149 (MQ=255)
                                                                        aCAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCAc                                                                      <  3:366950/149‑1 (MQ=255)
                                                                         cAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCAcc                                                                     <  6:252969/149‑1 (MQ=255)
                                                                            gAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGc                                                                                      <  4:412307/129‑1 (MQ=255)
                                                                            gAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGc                                                                                      >  3:412307/1‑129 (MQ=255)
                                                                             aaCGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAgg                                                                                             >  2:167499/1‑121 (MQ=255)
                                                                             aaCGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAgg                                                                                             <  1:167499/121‑1 (MQ=255)
                                                                             aaCGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCagag                                                                 <  8:164714/149‑1 (MQ=255)
                                                                               cGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCg                                                               <  5:267415/149‑1 (MQ=255)
                                                                                              ggCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCa                                                >  8:142530/1‑149 (MQ=255)
                                                                                               gCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAg                                               >  7:58353/1‑149 (MQ=255)
                                                                                               gCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCACCAGCCGAAAGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCa                                                <  8:272941/148‑1 (MQ=255)
                                                                                                 cAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGcc                                             >  2:199932/1‑149 (MQ=255)
                                                                                                  aGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCt                                            <  5:114118/149‑1 (MQ=255)
                                                                                                  aGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCt                                            <  3:214174/149‑1 (MQ=255)
                                                                                                   gTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTg                                           >  2:68425/1‑149 (MQ=255)
                                                                                                      ctcACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCTAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGa                                         <  6:381424/148‑1 (MQ=255)
                                                                                                      ctcACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCTAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGa                                         <  6:374527/148‑1 (MQ=255)
                                                                                                         aCCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGaa                                        <  5:272129/146‑1 (MQ=255)
                                                                                                         aCCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGaa                                        <  1:199932/146‑1 (MQ=255)
                                                                                                          cccGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGc                                              <  6:82445/139‑1 (MQ=255)
                                                                                                          cccGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGc                                              >  5:82445/1‑139 (MQ=255)
                                                                                                          cccGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTc                                    <  5:360228/149‑1 (MQ=255)
                                                                                                           ccGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCt                                   <  8:11359/149‑1 (MQ=255)
                                                                                                                  aCTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCg                            <  4:361725/149‑1 (MQ=255)
                                                                                                                   cTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGc                           <  1:17738/149‑1 (MQ=255)
                                                                                                                       gTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGaa                       <  7:142530/149‑1 (MQ=255)
                                                                                                                        tcatcaGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAAc                      <  3:237778/149‑1 (MQ=255)
                                                                                                                          atcaGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAACAt                    >  5:364660/1‑149 (MQ=255)
                                                                                                                          atcaGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAACAt                    >  3:282252/1‑149 (MQ=255)
                                                                                                                                      ccGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAACATGTCGCCGCcgcg        >  7:5126/1‑149 (MQ=255)
                                                                                                                                            gCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAACATGTCGCCGCCGCGGCGCCt  >  6:404661/1‑149 (MQ=255)
                                                                                                                                            gCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCATGCTGCTCCAGGGGAACGCCTTCTGAGAGCCTAACCAGAGCGGCCGGACGCAGACCAGACTGGATATTCGAAGGAGGGCGAACATGTCGCCGCAgcgtcggag  >  4:319155/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TCTCGACGGTGTGGGCTTGCGTGCGCTTGATCGCCGAGACCATCGCAACGCTGCCGATCGCTCTGTATGAGGACAAGAACGGCGCGCCGGTGGTGGCCAGTTCTCACCCGGTCAACTTCGTCATCAGCCAGCAGCCGAACGCCGATCAAACCCCGGTGGAGTTCTGGGAGAACGTCATGGCTAGCCTGCTGCTCCAGGGGAACGCCTTCTGCGAGCCTCACCAGAGCGGCCGGACGCTGACCAGCCTGGAATTTCTACTGCCGCAGAACATGTCGCCGCCGCGGCGCCT  >  NC_002947/1758986‑1759274

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