Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,174,596 A→G S190S (AGT→AGC)  ycfT ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,5960AG98.5% 195.2 / ‑4.3 68S190S (AGT→AGC) ycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (38/28);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

AGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAG  >  NC_000913/1174448‑1174737
                                                                                                                                                    |                                                                                                                                              
agaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTATTGCGTCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCt                                                                                                                                              <  1:385580/151‑1 (MQ=255)
   tAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGtt                                                                                                                                           >  2:249544/1‑151 (MQ=255)
   tAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGtt                                                                                                                                           >  2:542875/1‑151 (MQ=255)
    aGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGtt                                                                                                                                           <  2:693427/150‑1 (MQ=255)
       aaaCAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAGGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCAtt                                                                                                                                       >  2:251960/1‑151 (MQ=255)
          caacCAGGCACCGACCGCCAGAACGGAAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGACATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCAtt                                                                                                                                       >  2:363545/1‑148 (MQ=255)
             ccAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCccac                                                                                                                                >  1:118075/1‑151 (MQ=255)
                  cACCGCCCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCAGCTCTTTAACACAGGACATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACgg                                                                                                                           >  1:222136/1‑151 (MQ=255)
                     cgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGt                                                                                                                        >  1:412388/1‑151 (MQ=255)
                          gCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGaa                                                                                                                   <  2:573154/151‑1 (MQ=255)
                               aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGaaa                                                                                                              <  1:89672/151‑1 (MQ=255)
                                 cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAAtt                                                                                                            <  2:452654/151‑1 (MQ=255)
                                    tAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAAc                                                                                                         >  2:49092/1‑151 (MQ=255)
                                    tAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAAc                                                                                                         >  1:711944/1‑151 (MQ=255)
                                         aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCca                                                                                                    <  1:312657/151‑1 (MQ=255)
                                          aagaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCcac                                                                                                   >  2:208837/1‑151 (MQ=255)
                                           agaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCcac                                                                                                   >  1:783624/1‑150 (MQ=255)
                                             aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCacac                                                                                                 <  1:234964/150‑1 (MQ=255)
                                             aagCCATCAGCAGAGGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCAAACCCATGCGCCAGGGCTGTAATAAGGAAAATTGCGGACCCCGCTGTTCATC‑CCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACt                                                                                                <  1:630409/151‑1 (MQ=255)
                                                ccATCAGCAGATGGCGGCGCACCGGCCCCTTTTTAACACAGGTAATTATTGTTGCGCAAACCCATGCCCCAAGGCTGAAAAAAGGAAAATTGCGGATCCCGCTGTTCATT‑CCCCACCCCGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                               <  1:35794/149‑1 (MQ=255)
                                                  aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcag                                                                                             >  2:327678/1‑149 (MQ=255)
                                                   tCAGCAGAGGGCGGCGCAACGGCCCCTCTTTAACACGGTTCATTATTGTTGCGCAAACCCATGCCCCAAGGCTGTAATAAGGCAAATTGCGAATCACGCTGTTCATC‑CCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTcagca                                                                                           <  2:593055/150‑1 (MQ=255)
                                                   acagcagAGGGCGGCGAAACGGCCCCTCTTTAACACAGGTCATTTTTGTTGCGCCAAACCAGGCGCCAGGGCTGAAATAAGCCAAATGGCAAACCCCGCTGTTCTTC‑CCCCCCCACGGCGGGGGAAAGAAATTAACCGCCACACTcagcag                                                                                          <  1:723939/150‑1 (MQ=255)
                                                       cagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTac                                                                                       >  2:366165/1‑150 (MQ=255)
                                                        agagGGCGGCGCAACGGCCCCTTTTAAACACAGGTCATTTTTGTTGCCCCAACCCATGCGCAAAGGCTGAAATAAGGAAAATTGGGGATCACGCTTTTCATC‑CCCCCCCCCGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTaca                                                                                      <  2:447859/150‑1 (MQ=255)
                                                        agATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTaca                                                                                      >  1:831489/1‑150 (MQ=255)
                                                            ggcggcCCAACGGCCCCTTTTAAACACAGGTCTTTATTTTTGCGCAAACCCAGGCGCCAAGGCTGAAATAAGGCAAATTGCGGATCCCGCTGTTCATC‑CCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTacaaacaa                                                                                 <  2:803006/151‑1 (MQ=255)
                                                                 cgcAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                             >  2:416323/1‑150 (MQ=255)
                                                                          aCCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGTTGTAATAAGGCAAATTGCGGATCACGCTGTTCCTT‑CCCCACCACGGCGTGGGAACGAAAATAACCGCGACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                    >  1:71677/1‑150 (MQ=255)
                                                                             tctTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcg                                                                >  2:564859/1‑151 (MQ=255)
                                                                               tttAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                              >  2:769724/1‑151 (MQ=255)
                                                                               ttaaaaacaGTTCATTATTGTTGCGCCAAACCATGCGCCAGGGCTATAATAAGGAAAATGGCGGATCCCGCTGTTCATC‑CCCCACCCCGGCGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcgc                                                               <  1:117128/144‑1 (MQ=255)
                                                                                 tAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                              >  1:576332/1‑149 (MQ=255)
                                                                                 tAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                              >  1:290898/1‑149 (MQ=255)
                                                                                  aaCACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                              >  2:651819/1‑148 (MQ=255)
                                                                                  aaCACAGGTAATTATTGTTGCGCAAAACAATGCGCAAAGGTTGTAAAAAGGAAAATTGCGGATCACGTTTTTCATC‑CCCCACCCCGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                   <  1:293289/143‑1 (MQ=255)
                                                                                      caGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                       >  1:143232/1‑151 (MQ=255)
                                                                                        ggactttattGTTCCGCAAACCCATGCCCAAAGGTTGAAATAAGGAAATTTGCGTATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                     <  2:454354/146‑1 (MQ=255)
                                                                                        ggTCTTTATTGTTCCGCAAAACCCGCCGCCAGGGCTGTAATAGGGCAAGTTGCGTATCGCGCTG‑TCATTCCCCCACCCGGGCGTGGGAACGAAATTATCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                     <  1:100997/151‑1 (MQ=255)
                                                                                         gTCATTGTTGTCGGCCAAAACCATGCCCCAAGACAGTAAGAAGGCAAGTTGCGGCACCGGCTCTGCATC‑CCCCCCAACGGCGTGGAACAGAAAATAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                    <  1:53414/151‑1 (MQ=255)
                                                                                         gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGCCGGCTa                                                     >  1:605658/1‑150 (MQ=255)
                                                                                           cATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                   >  1:415852/1‑150 (MQ=255)
                                                                                            attattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                  >  2:441960/1‑150 (MQ=255)
                                                                                                ttgttgCGCCAAACCATGCGCCAAGGCTGAAAAAAGGAAAATTGCGGATCACGCTGTTCATC‑CCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                              <  1:49070/150‑1 (MQ=255)
                                                                                                ttgttgCCCAAACCCATGCGCCAAGGCTGTAATAAGGAAAATGGGGGATCACGCTGTTCATC‑CCCCCCCACGGCGTGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                              <  2:155940/150‑1 (MQ=255)
                                                                                                   ttgCGCAAAACCATGCGCAAGGGCTGTAATAAGGAAAATTGCGGATCACGCTGTTCATC‑CCCCCCCACGGGGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAc                                           <  1:162796/150‑1 (MQ=255)
                                                                                                     gcgcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                         >  1:551682/1‑150 (MQ=255)
                                                                                                      cgcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                       >  1:295303/1‑151 (MQ=255)
                                                                                                       gcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                      >  2:335900/1‑151 (MQ=255)
                                                                                                          aacccagcccccAGGGCTGTAAAAAGGCAATTTGCGAATCACGTTGTTCTTC‑CCCCCCCCCGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                    <  2:349857/140‑1 (MQ=255)
                                                                                                          aaaCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                   >  2:699742/1‑151 (MQ=255)
                                                                                                          aaaCCATGCGCCAAGGCTGTAAAAAGGAAAATGGCGGATCCCGTTGTTCATC‑CCCCCCCGCGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                    <  1:152104/150‑1 (MQ=255)
                                                                                                             ccATGCGCCAGGGCTGTAATAAGGAAAATTGCGGATCACGCTTTTCATC‑CCCCCCCACGGCGTGGGAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                <  1:233110/151‑1 (MQ=255)
                                                                                                             ccATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaa                                 >  2:805721/1‑150 (MQ=255)
                                                                                                              cATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                >  1:130402/1‑150 (MQ=255)
                                                                                                                   gcCAAGGCTGTAATAAGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACTGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                            >  1:491662/1‑149 (MQ=255)
                                                                                                                      aaGGCTGTAAAAAGGCAAATTGCGGATCCCGCTGTTCATC‑CCCCACCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                       <  2:248116/151‑1 (MQ=255)
                                                                                                                        ggCTTTAAAAAAGAAAATGGGGGACCACGTTGTTCATC‑CCCCCCCGCGGGGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                            <  1:4769/144‑1 (MQ=255)
                                                                                                                             taataaGGCAAATTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                 >  2:445628/1‑150 (MQ=255)
                                                                                                                                 aaGGAAAATGGGGGATCCCGCTGTTCATC‑CCCCCCCACGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg            <  1:551691/151‑1 (MQ=255)
                                                                                                                                  aGGCAAATTGCGGATCACACTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc              >  2:174600/1‑148 (MQ=255)
                                                                                                                                   ggaaaaTTGCGGACCACGCTTTTCATC‑CCCCACCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt          <  2:676093/148‑1 (MQ=255)
                                                                                                                                      aaaTTGCGGATCACGCTGATCATC‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCAAGACGGCTAAAAATTTTACATAACACGAAATAGACAATTAACGCATACAGATa         >  1:675899/1‑149 (MQ=255)
                                                                                                                                       aaTTGGGGATCACTCTTTTCATC‑CCCCCCCACGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca      <  1:335877/151‑1 (MQ=255)
                                                                                                                                       aaTTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc       >  2:177829/1‑150 (MQ=255)
                                                                                                                                        aTTGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac     >  1:850138/1‑151 (MQ=255)
                                                                                                                                          tGCGGATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcaca    >  2:362511/1‑150 (MQ=255)
                                                                                                                                              gATCACGCTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGCACAAAAAAGGCGAATAGTGGCAGCGCAAGACGGCTAAAAATTTTACATACCACGAAATATAAAATTAAAGCATACAGATACAACAAg  >  1:736849/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATT‑CCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAG  >  NC_000913/1174448‑1174737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑