Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,799 A→G 5.5% Y190H (TAC→CAC)  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,7990AG 5.5% 227.4 / 2.6 73Y190H (TAC→CAC) yedIconserved inner membrane protein YedI
Reads supporting (aligned to +/- strand):  ref base A (40/29);  new base G (2/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTT  >  NC_000913/2028651‑2028946
                                                                                                                                                    |                                                                                                                                                   
gCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAg                                                                                                                                                    <  1:1220397/150‑1 (MQ=255)
 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                                                                                        >  2:1175172/1‑150 (MQ=255)
                                                            ttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgc                                                                                        <  1:818290/150‑1 (MQ=255)
                                                            ttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgc                                                                                        <  1:823299/150‑1 (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:636012/1‑150 (MQ=255)
                                                                                                                                                  gTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCtt  >  2:1317303/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCAGCGTGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTT  >  NC_000913/2028651‑2028946

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: