Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,801 A→G 5.7% V189A (GTC→GCC)  yedI ← conserved inner membrane protein YedI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,8010AG 5.7% 201.6 / 2.6 71V189A (GTC→GCC) yedIconserved inner membrane protein YedI
Reads supporting (aligned to +/- strand):  ref base A (39/27);  new base G (2/2);  total (42/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 = 7.12e-01

CCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTC  >  NC_000913/2028652‑2028949
                                                                                                                                                     |                                                                                                                                                    
ccAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGa                                                                                                                                                      <  2:97474/150‑1 (MQ=255)
      tGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGa                                                                                                                                                      >  1:706775/1‑144 (MQ=255)
      tGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGa                                                                                                                                                      <  2:706775/144‑1 (MQ=255)
      tGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCa                                                                                                                                                <  1:543318/150‑1 (MQ=255)
        ccGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGCCGCCCACg                                                                                                                                              >  1:1121938/1‑150 (MQ=255)
         cGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACgg                                                                                                                                             <  2:938976/150‑1 (MQ=255)
              aTCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATGAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACt                                                                                                                                        <  1:564009/150‑1 (MQ=255)
               tCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTGGGCGTCTACGGTctggt                                                                                                                                       <  1:17063/150‑5 (MQ=255)
               tCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTGGGCGTCTACGGTctggt                                                                                                                                       <  1:1048050/150‑5 (MQ=255)
                 gATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCACCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTAcc                                                                                                                                     >  1:1227882/1‑150 (MQ=255)
                     aCGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAgcg                                                                                                                                 >  2:1278087/1‑150 (MQ=255)
                      cGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAgcgc                                                                                                                                >  1:46082/1‑150 (MQ=255)
                        cTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGa                                                                                                                              >  1:188382/1‑150 (MQ=255)
                        cTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGa                                                                                                                              >  1:74802/1‑150 (MQ=255)
                                aGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGaa                                                                                                                      <  2:1065438/150‑1 (MQ=255)
                                aGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGaa                                                                                                                      >  1:835826/1‑150 (MQ=255)
                                 gCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGaaa                                                                                                                     >  2:976798/1‑150 (MQ=255)
                                   cAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAagc                                                                                                                   <  1:1354724/150‑1 (MQ=255)
                                      ggCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAagcacc                                                                                                                <  1:292779/150‑1 (MQ=255)
                                          cgATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAagcaccagca                                                                                                            >  1:126623/1‑150 (MQ=255)
                                          cgATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAagcaccagca                                                                                                            >  1:586453/1‑150 (MQ=255)
                                             tAATCATTAACGCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTAAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAgcacctgcacct                                                                                                         <  2:996983/150‑1 (MQ=255)
                                               aTCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGa                                                                                                       >  1:1067512/1‑150 (MQ=255)
                                                          tttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACgg                                                                                            >  1:589003/1‑150 (MQ=255)
                                                           ttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgc                                                                                           <  1:818290/150‑1 (MQ=255)
                                                           ttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgc                                                                                           <  1:823299/150‑1 (MQ=255)
                                                           ttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgc                                                                                           >  2:1175172/1‑150 (MQ=255)
                                                             aCCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgcgc                                                                                         >  2:1171207/1‑150 (MQ=255)
                                                                 aaTGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCg                                                                                     <  1:837908/150‑1 (MQ=255)
                                                                    gCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCCCGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGcc                                                                                  <  2:572001/150‑1 (MQ=255)
                                                                     ccTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCCCGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGccc                                                                                 >  2:1157368/1‑149 (MQ=255)
                                                                      cTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCcac                                                                                <  2:707611/150‑1 (MQ=255)
                                                                       tGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCcaca                                                                               >  1:1116656/1‑150 (MQ=255)
                                                                            cAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATccc                                                                          <  2:1011445/150‑1 (MQ=255)
                                                                                 cgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATACCCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                     >  2:1195853/1‑150 (MQ=255)
                                                                                 cgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                     <  2:188382/150‑1 (MQ=255)
                                                                                 cgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                     >  2:15944/1‑150 (MQ=255)
                                                                                 cgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                     <  1:1175172/150‑1 (MQ=255)
                                                                                   cTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTg                                                                   >  2:371162/1‑150 (MQ=255)
                                                                                        tttttCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGc                                                              >  1:495597/1‑150 (MQ=255)
                                                                                         ttttCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGcg                                                             >  2:1055545/1‑150 (MQ=255)
                                                                                           ttCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGcgac                                                           <  2:710921/150‑1 (MQ=255)
                                                                                                    ccaATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGc                                                  <  2:74802/150‑1 (MQ=255)
                                                                                                     caATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCa                                                 >  1:815515/1‑150 (MQ=255)
                                                                                                          ccccAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGc                                                  >  1:1431310/1‑144 (MQ=255)
                                                                                                          ccccAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGc                                                  <  2:1431310/144‑1 (MQ=255)
                                                                                                          ccccAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACaa                                            >  2:866295/1‑150 (MQ=255)
                                                                                                              aGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAAtata                                        >  2:46079/1‑150 (MQ=255)
                                                                                                                gTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATaaa                                      >  1:625413/1‑150 (MQ=255)
                                                                                                                  catcaATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAAt                                    >  2:164494/1‑150 (MQ=255)
                                                                                                                   atcaATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATc                                   >  1:796222/1‑150 (MQ=255)
                                                                                                                        tCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTAc                              <  2:126623/150‑1 (MQ=255)
                                                                                                                             aCGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATc                         >  2:882882/1‑150 (MQ=255)
                                                                                                                             aCGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATc                         >  1:640594/1‑150 (MQ=255)
                                                                                                                              cGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCg                        >  1:1260260/1‑150 (MQ=255)
                                                                                                                                aTAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGcc                      <  2:1227882/150‑1 (MQ=255)
                                                                                                                                  aaCCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGcccc                    <  1:882882/150‑1 (MQ=255)
                                                                                                                                   aCCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAAtata                                        <  1:680489/129‑1 (MQ=255)
                                                                                                                                   aCCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAAtata                                        >  2:680489/1‑129 (MQ=255)
                                                                                                                                    ccccTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCGCCtt                  >  2:1330565/1‑150 (MQ=255)
                                                                                                                                     cccTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCttt                 >  2:896566/1‑150 (MQ=255)
                                                                                                                                      tggtagtgACCGTGGGCGTCTACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCCGACAATATAAAATCGGTACGAATCGCCCCtttt                >  2:1345110/8‑150 (MQ=255)
                                                                                                                                      ccTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCtttt                >  2:1415132/1‑150 (MQ=255)
                                                                                                                                         tagtgACCGTGGGCGTCTACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTAtt             >  1:249512/5‑150 (MQ=255)
                                                                                                                                          ccAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTAttt            >  2:474384/1‑150 (MQ=255)
                                                                                                                                          ccAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCATGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTAttt            <  2:848666/150‑1 (MQ=255)
                                                                                                                                           cAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTAtttt           <  1:1195853/150‑1 (MQ=255)
                                                                                                                                             gACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGt         <  1:371162/150‑1 (MQ=255)
                                                                                                                                               ccGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTcc       <  2:135571/150‑1 (MQ=255)
                                                                                                                                                cGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCt      >  2:53083/1‑150 (MQ=255)
                                                                                                                                                 gTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCtt     >  2:1317303/1‑150 (MQ=255)
                                                                                                                                                 gTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCtt     >  2:636012/1‑150 (MQ=255)
                                                                                                                                                    gACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTc  <  2:589003/150‑1 (MQ=255)
                                                                                                                                                    gACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTc  >  1:750292/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTC  >  NC_000913/2028652‑2028949

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: