Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,806 C→T 5.6% V187V (GTG→GTA 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,8060CT 5.6% 223.4 / 2.3 71V187V (GTG→GTAyedIconserved 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 = 6.50e-01

TGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTCAAACT  >  NC_000913/2028658‑2028954
                                                                                                                                                    |                                                                                                                                                    
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)
                                                                                                                                                   ccACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTCAAACt  <  1:976798/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCCGACAATCGATAACGCTTTCATCAGCCAGGGCGCGATAATCAATAACCCTTTACCTAATGCCTGCATCAGCGCGCTGGATTTTTCCGCCAGCCAATACCCCAGGTCATCAATCTTAACGATAACCCCTACCAGACCGTAGACGCCCACGGTCACTACCAGCGCGATGCCTGAAAGCACCAGCACCTGATTAAGCAACGGCGCTTCGGCCACAATCCCCAGCGTGATGGCGACGATTTCCGCAGACAATATAAAATCGGTACGAATCGCCCCTTTTATTTTGTCCTTTTCAAACT  >  NC_000913/2028658‑2028954

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: