Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,804 C→T 5.6% G188D (GGC→GAC)  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,8040CT 5.6% 227.7 / 2.1 71G188D (GGC→GAC) yedIconserved inner membrane protein YedI
Reads supporting (aligned to +/- strand):  ref base C (40/27);  new base T (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 = 9.85e-01

TGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTC  >  NC_000913/2028658‑2028949
                                                                                                                                                  |                                                                                                                                                 
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:1048050/150‑5 (MQ=255)
         tCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTGGGCGTCTACGGTctggt                                                                                                                                       <  1:17063/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:823299/150‑1 (MQ=255)
                                                     ttACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGgc                                                                                           <  1:818290/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                                                                     <  1:1175172/150‑1 (MQ=255)
                                                                           cgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                     <  2:188382/150‑1 (MQ=255)
                                                                           cgcTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCg                                                                     >  2:15944/1‑150 (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                         >  1:640594/1‑150 (MQ=255)
                                                                                                                       aCGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATc                         >  2:882882/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)
                                                                                                                                              gACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTc  <  2:589003/150‑1 (MQ=255)
                                                                                                                                              gACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTc  >  1:750292/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTC  >  NC_000913/2028658‑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: