Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,438,548 T→C V200A (GTG→GCG)  flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,5480TC100.0% 284.5 / NA 80V200A (GTG→GCG) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/40);  total (40/40)

AGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAA  >  NC_000913/2438413‑2438667
                                                                                                                                       |                                                                                                                       
agcagcGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGaa                                                                                                                      >  1:152879/1‑139 (MQ=255)
       cgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               <  2:373438/139‑1 (MQ=255)
              accgaccgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTaa                                                                                                        >  1:102160/1‑139 (MQ=255)
                      tCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCg                                                                                                >  2:1363762/1‑139 (MQ=255)
                                ccGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTa                                                                                      >  2:484948/1‑139 (MQ=255)
                                 cggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTaa                                                                                     >  2:1322813/1‑139 (MQ=255)
                                 cggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTaa                                                                                     >  1:1658144/1‑139 (MQ=255)
                                  ggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                <  2:229514/111‑1 (MQ=255)
                                  ggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                >  1:229514/1‑111 (MQ=255)
                                       aGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               <  2:715753/107‑1 (MQ=255)
                                       aGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               >  1:715753/1‑107 (MQ=255)
                                       aGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCt                                                                               >  2:514997/1‑139 (MQ=255)
                                             aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                         <  1:430836/139‑1 (MQ=255)
                                             aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                         <  2:148649/139‑1 (MQ=255)
                                             aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                         <  2:1658144/139‑1 (MQ=255)
                                             aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                         <  1:1322813/139‑1 (MQ=255)
                                              tGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                               <  1:152698/116‑1 (MQ=255)
                                              tGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                               >  2:152698/1‑116 (MQ=255)
                                              tGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAcc                                                                        >  1:1202374/1‑139 (MQ=255)
                                                 gcTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTgg                                                                     >  2:10641/1‑139 (MQ=255)
                                                   tAAAACAGCTGGTGACCAAACTTTCGGCGCCAACCGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGCCGAAAGAAGGTGAGTTAATTCCAGCGAAACTGTTTAACCATCAGGTGACCTGGCt                                                                   >  2:310026/1‑139 (MQ=255)
                                                          gCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                            >  2:879174/1‑139 (MQ=255)
                                                           cTGGTGACCAAACTTGCTGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               >  2:966861/1‑87 (MQ=255)
                                                           cTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               <  1:966861/87‑1 (MQ=255)
                                                           cTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcgc                                                           >  2:1032096/1‑139 (MQ=255)
                                                             ggTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                             >  1:862148/1‑103 (MQ=255)
                                                             ggTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                             <  2:862148/103‑1 (MQ=255)
                                                             ggTGACCAAACATGCGGCGGCAACGGGGAAACCCAGCAAACAGATCTGGCAACCGATGCTGGAACTCTCCGGGGCGGAAGATGGCGAGTTAATTCCAGCGAACCTGTTTATCCATCTGGTGACCTGGCTCCAGCCGCGt                                                         >  2:113120/1‑139 (MQ=255)
                                                              gTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTc                                                        <  1:39213/139‑1 (MQ=255)
                                                                       cTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTaa                                               <  1:514997/139‑1 (MQ=255)
                                                                         tGCGGCGGCAGCGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTGTCCGGGGCGAAAGATGGCGAGTTAATGCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGc                                             <  1:113120/139‑1 (MQ=255)
                                                                           cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaa                                                                                              >  1:318534/1‑88 (MQ=255)
                                                                           cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaa                                                                                              <  2:318534/88‑1 (MQ=255)
                                                                                gcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                >  2:1453256/1‑65 (MQ=255)
                                                                                gcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                <  1:1453256/65‑1 (MQ=255)
                                                                                       gggAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCTGACGCTAAGCCAGCAAAATACGcc                               <  1:1032096/139‑1 (MQ=255)
                                                                                       gggAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGcc                               <  2:152879/139‑1 (MQ=255)
                                                                                            cccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGc                          <  1:1363762/139‑1 (MQ=255)
                                                                                             ccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                         <  1:431624/139‑1 (MQ=255)
                                                                                              cAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAg                                          <  2:584605/121‑1 (MQ=255)
                                                                                              cAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAg                                          >  1:584605/1‑121 (MQ=255)
                                                                                               aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTgg                       >  1:88088/1‑139 (MQ=255)
                                                                                               aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTgg                       <  2:1581793/139‑1 (MQ=255)
                                                                                               aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTgg                       >  1:1581793/1‑139 (MQ=255)
                                                                                               aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTgg                       >  2:797950/1‑139 (MQ=255)
                                                                                                gCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                               <  2:755056/66‑1 (MQ=255)
                                                                                                gCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                               >  1:755056/1‑66 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGa                                                                                                             >  1:526570/1‑51 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGa                                                                                                             <  2:526570/51‑1 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                         >  1:669265/1‑87 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                         <  2:669265/87‑1 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAg                                          <  2:1338079/118‑1 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAg                                          >  1:1338079/1‑118 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                     <  1:879174/139‑1 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                     <  2:88088/139‑1 (MQ=255)
                                                                                                 cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCNACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                     >  1:1390825/1‑139 (MQ=255)
                                                                                                  aaaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                    <  1:484948/139‑1 (MQ=255)
                                                                                                   aaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGt                                                                            >  2:693527/1‑82 (MQ=255)
                                                                                                   aaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGt                                                                            <  1:693527/82‑1 (MQ=255)
                                                                                                      aGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               <  2:986363/44‑1 (MQ=255)
                                                                                                      aGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                               >  1:986363/1‑44 (MQ=255)
                                                                                                       gATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                            >  2:463770/1‑94 (MQ=255)
                                                                                                       gATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                            <  1:463770/94‑1 (MQ=255)
                                                                                                       gATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                      <  2:886207/132‑1 (MQ=255)
                                                                                                       gATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                      >  1:886207/1‑132 (MQ=255)
                                                                                                           tGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                             <  1:1086471/57‑1 (MQ=255)
                                                                                                           tGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                             >  2:1086471/1‑57 (MQ=255)
                                                                                                            ggCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGAt          <  2:102160/139‑1 (MQ=255)
                                                                                                             gCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGcc                               <  2:1132136/117‑1 (MQ=255)
                                                                                                             gCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGcc                               >  1:1132136/1‑117 (MQ=255)
                                                                                                                  cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAcc                                                                                   <  1:164596/60‑1 (MQ=255)
                                                                                                                  cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAcc                                                                                   >  2:164596/1‑60 (MQ=255)
                                                                                                                  cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGcc                               >  1:385277/1‑112 (MQ=255)
                                                                                                                  cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGcc                               <  2:385277/112‑1 (MQ=255)
                                                                                                                    aTGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTaa  >  1:265716/1‑139 (MQ=255)
                                                                                                                       cTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                             >  2:1076399/1‑45 (MQ=255)
                                                                                                                       cTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                             <  1:1076399/45‑1 (MQ=255)
                                                                                                                          gAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                         <  2:565629/110‑1 (MQ=255)
                                                                                                                          gAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                         >  1:565629/1‑110 (MQ=255)
                                                                                                                                  cGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                     >  1:48866/1‑106 (MQ=255)
                                                                                                                                  cGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGaa                     <  2:48866/106‑1 (MQ=255)
                                                                                                                                       |                                                                                                                       
AGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAA  >  NC_000913/2438413‑2438667

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: