Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,451,123 C→A V22V (GTC→GTA ytfR → putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,451,1230CA98.9% 346.0 / ‑3.8 88V22V (GTC→GTAytfRputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (41/46);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAG  >  NC_000913/4450975‑4451270
                                                                                                                                                    |                                                                                                                                                   
ggCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTaaa                                                                                                                                                   <  1:136206/151‑1 (MQ=255)
   gATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGc                                                                                                                                                <  2:265011/151‑1 (MQ=255)
   gATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGc                                                                                                                                                <  2:684985/151‑1 (MQ=255)
     tGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGt                                                                                                                                              >  2:1370000/1‑151 (MQ=255)
         gCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGa                                                                                                                                          <  2:53640/151‑1 (MQ=255)
         gCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGa                                                                                                                                          <  1:381199/151‑1 (MQ=255)
           aCTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGAc                                                                                                                                         >  1:1202899/1‑150 (MQ=255)
            cTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACa                                                                                                                                        >  1:406359/1‑150 (MQ=255)
              gAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACg                                                                                                                                     >  1:479203/1‑151 (MQ=255)
              gAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACg                                                                                                                                     <  2:254926/151‑1 (MQ=255)
              gAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACg                                                                                                                                     <  1:1311449/151‑1 (MQ=255)
               aaGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGt                                                                                                                                    >  2:293320/1‑151 (MQ=255)
               aaGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGt                                                                                                                                    <  2:1139948/151‑1 (MQ=255)
                aGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGtt                                                                                                                                   <  1:719016/151‑1 (MQ=255)
                  ccTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTg                                                                                                                                  >  2:842914/1‑150 (MQ=255)
                      aTGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACGCCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTc                                                                                                                             <  1:1083004/151‑1 (MQ=255)
                        gCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCa                                                                                                                            <  1:105693/150‑1 (MQ=255)
                         cGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGc                                                                                                                          <  2:463982/151‑1 (MQ=255)
                              cTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTgcg                                                                                                                     >  2:41870/1‑151 (MQ=255)
                                 aCCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCg                                                                                                                  <  1:1268281/151‑1 (MQ=255)
                                    tcgTCTTATCATTCCTACTGGGAGCTCGCCTTACACCGGGGGAACCCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTCAGTAAATTTTTCTTCGGCGTAAAAGCGTTAGACATCGTTGATTTCAGCCTGCTCCGt                                                                                                                 >  1:626853/2‑149 (MQ=255)
                                     cgTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGc                                                                                                              >  1:241306/1‑151 (MQ=255)
                                       tCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCg                                                                                                             >  2:287669/1‑150 (MQ=255)
                                        cTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGa                                                                                                            >  2:1271303/1‑150 (MQ=255)
                                          tATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGaaa                                                                                                          <  1:439604/150‑1 (MQ=255)
                                          tATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGaa                                                                                                           <  1:633261/149‑1 (MQ=255)
                                                 ccTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAGGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGGATAAGTAAATTTTTACCCGGCGTAACAGCGTTAGACAACGTTGATTTCAGCCTGCCCCGTGGCGAAATCACGgc                                                                                                   >  1:1320068/1‑150 (MQ=255)
                                                    aCTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCgctg                                                                                               >  1:984602/1‑151 (MQ=255)
                                                      tGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCgctgct                                                                                             >  1:624668/1‑151 (MQ=255)
                                                       gggAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCgctgct                                                                                             <  2:830992/150‑1 (MQ=255)
                                                       gggAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTc                                                                                            >  2:125674/1‑151 (MQ=255)
                                                         gAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCgg                                                                                          <  2:406472/151‑1 (MQ=255)
                                                          aGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGt                                                                                         <  2:826179/151‑1 (MQ=255)
                                                           gCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTg                                                                                        <  1:667007/151‑1 (MQ=255)
                                                              cGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGaa                                                                                      >  1:380328/1‑150 (MQ=255)
                                                                 tttACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACg                                                                                  >  2:254024/1‑151 (MQ=255)
                                                                  ttACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAAc                                                                                   <  2:1373406/149‑1 (MQ=255)
                                                                  ttACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAAc                                                                                   <  1:100959/149‑1 (MQ=255)
                                                                   tACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACggg                                                                                >  2:999552/1‑151 (MQ=255)
                                                                   tACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACgg                                                                                 <  1:1322085/150‑1 (MQ=255)
                                                                    acacCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACgggg                                                                               >  2:935045/1‑151 (MQ=255)
                                                                    acacCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACgggg                                                                               <  2:487384/151‑1 (MQ=255)
                                                                     cacCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGc                                                                              <  2:37521/151‑1 (MQ=255)
                                                                          ggggAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGaa                                                                         <  1:45778/151‑1 (MQ=255)
                                                                              aaaCCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATc                                                                     >  1:920027/1‑151 (MQ=255)
                                                                              aaaCCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATc                                                                     >  2:686817/1‑151 (MQ=255)
                                                                               aaCCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATc                                                                     >  2:1025289/1‑150 (MQ=255)
                                                                                  cATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAAc                                                                  >  1:1218139/1‑150 (MQ=255)
                                                                                  cATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACg                                                                 <  1:1332188/151‑1 (MQ=255)
                                                                                    tgtcgaccgacCACCACCAGGAGATCATCCGCACCGAAGGATTAAGTAAATTTTTCCCCAGCGTAAACGCGTAAGACAACGTTGATTTCAACCTGCGCCGTGGCGAAATCCTGGCGCTGCTCGGTGAAAACGGGGCGGGaaaagtaccgct                                                               >  1:1382409/4‑143 (MQ=255)
                                                                                          ccgaccAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCaa                                                         <  1:1116436/151‑1 (MQ=255)
                                                                                           cgaccAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCaaa                                                        <  2:58804/151‑1 (MQ=255)
                                                                                           cgaccAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCaa                                                         >  1:1160185/1‑150 (MQ=255)
                                                                                              ccAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCa                                                     <  2:140110/151‑1 (MQ=255)
                                                                                               cAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCAt                                                    >  2:346522/1‑151 (MQ=255)
                                                                                                aaCACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCAtt                                                   <  1:1013390/151‑1 (MQ=255)
                                                                                                   acCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCTTAAAAGTGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGTCGCGAAAATCATCGCTAATCAAAGCATTaa                                                 <  2:865777/150‑1 (MQ=255)
                                                                                                   acCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAAc                                                <  2:1079069/151‑1 (MQ=255)
                                                                                                    ccAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACt                                               <  2:1132454/151‑1 (MQ=255)
                                                                                                     cAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTg                                              <  1:1120433/151‑1 (MQ=255)
                                                                                                      aGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCAATAAATgg                                             >  1:1338641/1‑151 (MQ=255)
                                                                                                       ggAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTgg                                             >  2:240279/1‑150 (MQ=255)
                                                                                                       ggAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgt                                            <  2:1228973/151‑1 (MQ=255)
                                                                                                       ggAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgt                                            <  2:1076603/151‑1 (MQ=255)
                                                                                                          gaTCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGgtgt                                          >  2:660649/1‑150 (MQ=255)
                                                                                                                ccGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCa                                    <  2:1364539/150‑1 (MQ=255)
                                                                                                                   cACCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGc                                 <  1:1391209/150‑1 (MQ=255)
                                                                                                                    aCCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCg                               <  1:607603/151‑1 (MQ=255)
                                                                                                                    aCCGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCg                               >  2:932352/1‑151 (MQ=255)
                                                                                                                     ccGAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCg                               >  1:869438/1‑150 (MQ=255)
                                                                                                                     ccGAAGGATTAAGCAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATAAAAGCATTAACTGGGGTATAACACGCCg                               >  1:901953/1‑150 (MQ=255)
                                                                                                                       gAAGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATc                            >  2:666898/1‑151 (MQ=255)
                                                                                                                         aGGATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcgc                          >  2:54831/1‑151 (MQ=255)
                                                                                                                          ggATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcgcg                         <  2:1253776/151‑1 (MQ=255)
                                                                                                                          ggATTAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATcgcg                         >  2:1034619/1‑151 (MQ=255)
                                                                                                                             ttAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGGGGATACCACGCCGATCGCGGc                       >  2:26743/1‑150 (MQ=255)
                                                                                                                              tAAGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCAc                     >  1:1318981/1‑151 (MQ=255)
                                                                                                                               aaGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCAc                     <  2:1352344/150‑1 (MQ=255)
                                                                                                                                aGTAAATTTTTCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCa                   <  1:1369115/151‑1 (MQ=255)
                                                                                                                                          tCCCCGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATctggctgg          <  2:508716/150‑1 (MQ=255)
                                                                                                                                           ccccGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGa         >  2:739717/1‑150 (MQ=255)
                                                                                                                                             ccGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAgg      >  2:897926/1‑151 (MQ=255)
                                                                                                                                              cGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGc     >  2:631718/1‑151 (MQ=255)
                                                                                                                                              cGGCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGc     >  2:1291532/1‑151 (MQ=255)
                                                                                                                                               ggCGTACAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCCCCATCTGGCTGGAAGGcc    >  1:202532/1‑151 (MQ=255)
                                                                                                                                               ggCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGcc    <  2:1011507/151‑1 (MQ=255)
                                                                                                                                                gCGTAAAAGCGTTAGACAGCGTTGATTTCAGTCTGCGCCGTGGCGAAATCATGGCGATGCTCGGTGAAAACGGGGCGGGAAAAGCACCGGTAATCGAAGCATTAACTGGTGTATACCACGCCGATCGTGGCATCATCTGGCTGGAAGGCCa   <  1:1164839/151‑1 (MQ=255)
                                                                                                                                                gCGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCa   <  1:1266767/151‑1 (MQ=255)
                                                                                                                                                 cGTAAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAg  <  1:757448/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAG  >  NC_000913/4450975‑4451270

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