Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,963,048 T→C 5.2% Q668R (CAG→CGG)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,963,0480TC 5.2% 239.3 / 2.4 77Q668R (CAG→CGG) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base C (2/2);  total (36/41)
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

GCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTGGCTCCAG  >  NC_000913/1962902‑1963193
                                                                                                                                                  |                                                                                                                                                 
gcctgccACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTg                                                                                                                                                  <  5:146384/148‑1 (MQ=255)
 cctgccACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGc                                                                                                                                                 <  5:355423/148‑1 (MQ=255)
 cctgccACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCg                                                                                                                                                >  3:287681/1‑149 (MQ=255)
     ccACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGc                                                                                                                                              >  3:62520/1‑147 (MQ=255)
       acaTCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTAGACAGGACCACTAACTGCGGCAAGc                                                                                                                                          >  8:168469/1‑149 (MQ=255)
         aTCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTg                                                                                                                                        >  3:7948/1‑149 (MQ=255)
         aTCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTg                                                                                                                                        >  7:254250/1‑149 (MQ=255)
          tcccctccccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTgc                                                                                                                                       >  5:335417/1‑149 (MQ=255)
               tccccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTgcggcg                                                                                                                                   >  4:25747/1‑148 (MQ=255)
               tccccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTgcggcg                                                                                                                                   >  3:16417/1‑148 (MQ=255)
               tccccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTgcggcg                                                                                                                                   >  4:41550/1‑148 (MQ=255)
                ccccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTgcggcg                                                                                                                                   >  5:148610/1‑147 (MQ=255)
                 cccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGgcgc                                                                                                                                  <  8:342314/147‑1 (MQ=255)
                  ccGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAg                                                                                                                                >  1:306752/1‑148 (MQ=255)
                    ggcggcATGCACCACCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTAATTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGaa                                                                                                                             <  2:111658/149‑1 (MQ=255)
                      cggcTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGc                                                                                                                           <  2:309688/149‑1 (MQ=255)
                        gcTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGGTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCga                                                                                                                         <  6:235843/149‑1 (MQ=255)
                          ttGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCgaga                                                                                                                       >  1:218439/1‑149 (MQ=255)
                           tGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCgaga                                                                                                                       >  1:119685/1‑148 (MQ=255)
                                cagcagCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACaataa                                                                                                                 >  2:319940/1‑149 (MQ=255)
                                  gcagcGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATg                                                                                                               >  5:94419/1‑149 (MQ=255)
                                  gcagcGGAAACAATAATATTGCTAATAAGGTGCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATg                                                                                                               <  2:259967/149‑1 (MQ=255)
                                       ggAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACCGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGc                                                                                                          >  8:353299/1‑149 (MQ=255)
                                       ggAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACCGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGc                                                                                                          >  8:358327/1‑149 (MQ=255)
                                                 atatTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGt                                                                                                >  8:41832/1‑149 (MQ=255)
                                                    ttGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCa                                                                                             <  4:7948/149‑1 (MQ=255)
                                                          ataaGGTTCTCATTCTTTGCCGCCAATTGTCGCCGTCATGCGGATATGGCGGTTATCAGACAGGTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACcaa                                                                                         <  6:22916/147‑1 (MQ=255)
                                                          ataaGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGCCCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACcaac                                                                                        >  6:31032/1‑148 (MQ=255)
                                                          ataaGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACcaac                                                                                        <  2:119685/148‑1 (MQ=255)
                                                             aGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATa                                                                                    <  8:254250/149‑1 (MQ=255)
                                                             aGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATa                                                                                    <  2:218439/149‑1 (MQ=255)
                                                               gTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACt                                                                                  <  3:75400/149‑1 (MQ=255)
                                                                ttCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTg                                                                                 <  6:148610/149‑1 (MQ=255)
                                                                    cATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTggcgg                                                                             >  8:307322/1‑149 (MQ=255)
                                                                     attattTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTggcggc                                                                            <  6:94419/149‑1 (MQ=255)
                                                                     attattTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACCGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTggc                                                                               <  7:358327/146‑1 (MQ=255)
                                                                     attattTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACCGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTggc                                                                               <  7:353299/146‑1 (MQ=255)
                                                                         tttGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCa                                                                         <  7:378590/148‑1 (MQ=255)
                                                                           tGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCAcc                                                                       <  6:335417/148‑1 (MQ=255)
                                                                              cgccAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAg                                                                    >  1:251860/1‑148 (MQ=255)
                                                                                ccAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCAt                                                                 <  2:169106/149‑1 (MQ=255)
                                                                                ccAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCAt                                                                 <  3:157615/149‑1 (MQ=255)
                                                                                ccAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCAt                                                                 >  3:325071/1‑149 (MQ=255)
                                                                                  aaTTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCAtct                                                               >  2:148016/1‑149 (MQ=255)
                                                                                       tCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGGGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGa                                                          <  1:319940/149‑1 (MQ=255)
                                                                                         gCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACg                                                        <  8:109869/149‑1 (MQ=255)
                                                                                          ccGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGAGAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACgg                                                       <  4:286167/149‑1 (MQ=255)
                                                                                          ccGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACgg                                                       <  2:306752/149‑1 (MQ=255)
                                                                                            gTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGa                                                     <  1:148016/149‑1 (MQ=255)
                                                                                             tCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGa                                                     >  6:347882/1‑148 (MQ=255)
                                                                                               aTGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAg                                                  >  8:319670/1‑149 (MQ=255)
                                                                                                 gCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAgc                                                 >  4:66145/1‑148 (MQ=255)
                                                                                                  cGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAgcgc                                               <  6:123961/149‑1 (MQ=255)
                                                                                                         gTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCt                                          >  6:322952/1‑147 (MQ=255)
                                                                                                                 tCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTgcgc                                <  2:203191/149‑1 (MQ=255)
                                                                                                                     aCAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGa                                        >  8:191625/1‑137 (MQ=255)
                                                                                                                     aCAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGa                                        <  7:191625/137‑1 (MQ=255)
                                                                                                                        gTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACaataat                                                                                                                >  2:239901/1‑62 (MQ=255)
                                                                                                                        gTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACAAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTaaa                         <  7:243494/149‑1 (MQ=255)
                                                                                                                         ttCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAAc                        <  3:25747/149‑1 (MQ=255)
                                                                                                                           ccAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGa                      >  7:11040/1‑149 (MQ=255)
                                                                                                                            ccgATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCTCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTTAAGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGa                      <  2:301848/146‑1 (MQ=255)
                                                                                                                            cAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGAt                     <  4:325071/149‑1 (MQ=255)
                                                                                                                            cAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGGGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGa                      <  4:92644/148‑1 (MQ=255)
                                                                                                                              gATTCGACAGGCCCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCGTCCTGAGTTTGCGCCAGTAAACGATcc                   <  5:84234/149‑1 (MQ=255)
                                                                                                                              gATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATcc                   >  3:240962/1‑149 (MQ=255)
                                                                                                                                   gACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCa               <  4:240962/148‑1 (MQ=255)
                                                                                                                                   gACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCAACTCCTGACGGGATACCGCTTCCTGAGATTGCGCCAGTAAACGATCCGCCAg              >  4:75578/1‑149 (MQ=255)
                                                                                                                                     cAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGcc            <  5:322952/149‑1 (MQ=255)
                                                                                                                                      aGGACCACTAACTGCGGAAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACAAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCAGCCAGccc           >  1:7118/1‑149 (MQ=255)
                                                                                                                                        gACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTg         >  2:37408/1‑149 (MQ=255)
                                                                                                                                        gACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTg         >  4:195978/1‑149 (MQ=255)
                                                                                                                                             cTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGc                 >  7:213037/1‑136 (MQ=255)
                                                                                                                                             cTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGc                 <  8:213037/136‑1 (MQ=255)
                                                                                                                                              tAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTGGCTCCa   <  4:62520/149‑1 (MQ=255)
                                                                                                                                                aCTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTGGCTCCAg  <  5:347882/148‑1 (MQ=255)
                                                                                                                                                aCTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTGGCTCCAg  <  6:157740/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCCAGCATCTCCTGACGGGATAGCGCTTCCTGAGTTTGCGCCAGTAAACGATCCGCCAGCCCTGGCTCCAG  >  NC_000913/1962902‑1963193

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: