Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,226,509 C→T *254* (TAG→TAA yohF ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,226,5090CT100.0% 275.3 / NA 84*254* (TAG→TAAyohFputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/44);  total (40/44)

CGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAA  >  NC_000913/2226375‑2226641
                                                                                                                                      |                                                                                                                                    
cGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                  <  2:132231/139‑1 (MQ=255)
  tAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                  <  2:851994/137‑1 (MQ=255)
  tAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                  >  1:851994/1‑137 (MQ=255)
         aTTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                  <  1:814164/130‑1 (MQ=255)
         aTTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                  >  2:814164/1‑130 (MQ=255)
           tGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAAc                                                                                                                       >  1:1214722/1‑139 (MQ=255)
                 ggATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                 <  1:121782/139‑1 (MQ=255)
                       ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                           >  1:851806/1‑139 (MQ=255)
                       ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                           >  2:412869/1‑139 (MQ=255)
                                gccgcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                                                  <  1:718273/139‑1 (MQ=255)
                                gccgcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                                                  <  1:1581674/139‑1 (MQ=255)
                                 ccgcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgc                                                                                                 <  2:851806/139‑1 (MQ=255)
                                    cATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGTCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccgc                                                                                              <  2:1088348/139‑1 (MQ=255)
                                          aCATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        <  2:800084/139‑1 (MQ=255)
                                          aCATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        <  2:487991/139‑1 (MQ=255)
                                            aTCCGGATTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACt                                                                                      <  2:1623212/139‑1 (MQ=255)
                                            aTCCGACGTACCATGCCTGATACGACGCTTGCGCGTCTTTTCCTGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACt                                                                                      <  2:1657775/139‑1 (MQ=255)
                                                         tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                               >  2:249118/1‑117 (MQ=255)
                                                         tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                               <  1:249118/117‑1 (MQ=255)
                                                         tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                                         >  1:776542/1‑139 (MQ=255)
                                                           ccTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       <  1:640559/139‑1 (MQ=255)
                                                           ccTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       <  1:185262/139‑1 (MQ=255)
                                                            cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCgg                                                                      >  1:479622/1‑139 (MQ=255)
                                                             tGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGTCTGCCCggt                                                                     <  1:1110932/139‑1 (MQ=255)
                                                              gATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtg                                                                    >  2:39903/1‑139 (MQ=255)
                                                               aTACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtgg                                                                   >  1:1332796/1‑139 (MQ=255)
                                                                tACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtggt                                                                  <  1:412869/139‑1 (MQ=255)
                                                                    acgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTaa                                                              >  2:1335819/1‑139 (MQ=255)
                                                                     cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGtt                                                                                                                           >  1:935046/1‑77 (MQ=255)
                                                                     cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGtt                                                                                                                           <  2:935046/77‑1 (MQ=255)
                                                                     cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCTACTATCACCGACTGCCCGGTGGTGTAAt                                                             >  2:1259411/1‑139 (MQ=255)
                                                                     cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAAt                                                             >  2:1463453/1‑139 (MQ=255)
                                                                        ttGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCgg                                                                      >  2:1639157/1‑127 (MQ=255)
                                                                        ttGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCgg                                                                      <  1:1639157/127‑1 (MQ=255)
                                                                          gtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCAt                                                                                           >  1:1077996/1‑104 (MQ=255)
                                                                          gtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCAt                                                                                           <  2:1077996/104‑1 (MQ=255)
                                                                          gtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgcg                                                        <  1:716630/139‑1 (MQ=255)
                                                                                          cGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGcc                                        >  1:1045932/1‑139 (MQ=255)
                                                                                          cGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCACCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgc                                                         >  2:830493/1‑122 (MQ=255)
                                                                                          cGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCACCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgc                                                         <  1:830493/122‑1 (MQ=255)
                                                                                           gACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       <  2:601218/139‑1 (MQ=255)
                                                                                           gACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       <  2:1045932/139‑1 (MQ=255)
                                                                                           gACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAGCcac                                      >  1:1204822/1‑139 (MQ=255)
                                                                                            aCCAAAATATTATCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                      >  2:450692/1‑139 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       <  1:1656930/106‑1 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       >  2:1656930/1‑106 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       >  2:1656929/1‑106 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       <  1:1656929/106‑1 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                      >  1:651341/1‑139 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       <  1:1327469/138‑1 (MQ=255)
                                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       >  2:1327469/1‑138 (MQ=255)
                                                                                                     ttACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAgg                             <  2:1272159/139‑1 (MQ=255)
                                                                                                       aCCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGt                                                                                                                            <  1:155694/42‑1 (MQ=255)
                                                                                                       aCCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGt                                                                                                                            >  2:155694/1‑42 (MQ=255)
                                                                                                        cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                              >  2:33406/1‑119 (MQ=255)
                                                                                                        cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                              <  1:33406/119‑1 (MQ=255)
                                                                                                        cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       <  2:1181858/126‑1 (MQ=255)
                                                                                                        cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       >  1:1181858/1‑126 (MQ=255)
                                                                                                           aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                           >  2:1577863/1‑55 (MQ=255)
                                                                                                           aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                           <  1:1577863/55‑1 (MQ=255)
                                                                                                           aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa                       <  1:1463453/139‑1 (MQ=255)
                                                                                                           aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa                       <  2:1214722/139‑1 (MQ=255)
                                                                                                             tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                          <  1:278294/38‑1 (MQ=255)
                                                                                                             tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                          <  2:592242/38‑1 (MQ=255)
                                                                                                             tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                          >  2:278294/1‑38 (MQ=255)
                                                                                                             tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                          >  1:592242/1‑38 (MQ=255)
                                                                                                                  tatCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        <  2:798791/67‑1 (MQ=255)
                                                                                                                  tatCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        >  1:798791/1‑67 (MQ=255)
                                                                                                                  tatCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCAt                >  1:1059635/1‑139 (MQ=255)
                                                                                                                   atCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                 <  2:1357622/41‑1 (MQ=255)
                                                                                                                   atCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                 >  1:1357622/1‑41 (MQ=255)
                                                                                                                    tCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGc              >  1:1032624/1‑139 (MQ=255)
                                                                                                                       aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGtt           >  2:791640/1‑139 (MQ=255)
                                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATc                                                                                   >  2:227175/1‑65 (MQ=255)
                                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATc                                                                                   <  1:227175/65‑1 (MQ=255)
                                                                                                                           aCCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTgcgc       <  1:1335819/139‑1 (MQ=255)
                                                                                                                              tCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       >  2:348034/1‑72 (MQ=255)
                                                                                                                              tCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                       <  1:348034/72‑1 (MQ=255)
                                                                                                                               cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGgtgt                                                                >  2:176904/1‑78 (MQ=255)
                                                                                                                               cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGgtgt                                                                <  1:176904/78‑1 (MQ=255)
                                                                                                                               cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCaaa   >  1:484825/1‑139 (MQ=255)
                                                                                                                                aCTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCaaaa  <  1:450692/139‑1 (MQ=255)
                                                                                                                                  tGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                               <  1:1302742/44‑1 (MQ=255)
                                                                                                                                  tGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                               >  2:1302742/1‑44 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAA  >  NC_000913/2226375‑2226641

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