Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 335,166 T→A 19.5% intergenic (‑144/‑114) yahC ← / → yahD putative inner membrane protein/ankyrin repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913335,1660TA19.5% 372.0 / 89.8 174intergenic (‑144/‑114)yahC/yahDputative inner membrane protein/ankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base T (83/57);  new base A (15/19);  total (98/76)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

Only 100 of 187 total aligned reads displayed.

GTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAA  >  NC_000913/335067‑335265
                                                                                                   |                                                                                                   
  gCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAatat                                                                                                     >  2:3978526/1‑98 (MQ=255)
   cAGACTAGAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtattat                                                                                                  <  2:3655534/100‑1 (MQ=255)
    aGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATattatt                                                                                                 <  1:3052244/100‑1 (MQ=255)
    aGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATattatt                                                                                                 >  2:39345/1‑100 (MQ=255)
     gACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTc                                                                                                <  2:835246/100‑1 (MQ=255)
     gACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTc                                                                                                >  2:1467521/1‑100 (MQ=255)
      aCTATAGCGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTca                                                                                               <  2:2221342/100‑1 (MQ=255)
      aCTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTca                                                                                               <  1:2891744/100‑1 (MQ=255)
      aCTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGACAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGATGAAtatt                                                                                                    >  2:3366011/1‑95 (MQ=255)
       cTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTTCCGCACACTGCGAAAAGATTGTTGAATATTATTca                                                                                               <  2:905295/99‑1 (MQ=255)
        tataGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaac                                                                                             <  1:2109227/100‑1 (MQ=255)
        tataGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaa                                                                                              <  1:472191/99‑1 (MQ=255)
        tataGAGAAATAAAACCTGATTGAGTGAAAAGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaac                                                                                             <  2:370861/100‑1 (MQ=255)
         ataGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaaca                                                                                            >  1:3822935/1‑100 (MQ=255)
         ataGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaaca                                                                                            >  1:906661/1‑100 (MQ=255)
          taGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaacaa                                                                                           >  2:1414214/1‑100 (MQ=255)
          taGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaacaa                                                                                           >  2:3076450/1‑100 (MQ=255)
           agagAAATAAAACCTGATTGAGTGAAAAGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaa                                                                                           <  2:2308764/99‑1 (MQ=255)
            gagaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatc                                                                                         <  1:3017117/100‑3 (MQ=255)
             agaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaaa                                                                                        <  2:1217000/100‑1 (MQ=255)
             agaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaaa                                                                                        >  2:535647/1‑100 (MQ=255)
             agaAATAAAACCTAATTAAGTGAACGGGTATCGGGCCAAAGAAACAAATATTGCACTACCGTACACTGCGAAAAGATTGTTGAATATTATaaaaaaaaaa                                                                                        <  2:2537390/100‑1 (MQ=255)
              gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAg                                                                                       >  2:2776305/1‑100 (MQ=255)
              gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAg                                                                                       >  2:1774164/1‑100 (MQ=255)
              gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatctt                                                                                       <  2:3555319/100‑5 (MQ=255)
               aaaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttt                                                                                      <  2:3858662/100‑6 (MQ=255)
               aaaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttt                                                                                      <  1:2088421/100‑6 (MQ=255)
                aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatctttt                                                                                     <  1:1596241/100‑7 (MQ=255)
                aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatctttt                                                                                     <  1:3208202/100‑7 (MQ=255)
                 aTAAAACCTGCTTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCg                                                                                    >  2:4124279/1‑100 (MQ=255)
                 aTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCg                                                                                    <  1:2569550/100‑1 (MQ=255)
                 aTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCg                                                                                    >  1:1574538/1‑100 (MQ=255)
                 aTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttttc                                                                                    <  1:1616219/100‑8 (MQ=39)
                 aTAAAAACTGATTGAGTGAAAGGGGATCGGGTCAAATAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatcttttc                                                                                    <  2:3579076/100‑8 (MQ=38)
                  tAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCga                                                                                   <  2:2807339/100‑1 (MQ=255)
                   aaaaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCgag                                                                                  <  2:1295578/100‑1 (MQ=255)
                     aaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGAt                                                                                >  1:1453388/1‑100 (MQ=255)
                      aCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATa                                                                               >  1:1119367/1‑100 (MQ=255)
                       ccTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATAc                                                                              >  2:1753995/1‑100 (MQ=255)
                        cTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACGGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACg                                                                             >  2:2691327/1‑100 (MQ=255)
                         tGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACgg                                                                            >  2:3808803/1‑100 (MQ=255)
                           aTTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCg                                                                          <  2:2425602/100‑1 (MQ=255)
                           aTTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCg                                                                          <  1:779439/100‑1 (MQ=255)
                           aTTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCg                                                                          <  1:1632501/100‑1 (MQ=255)
                              gAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAt                                                                       >  2:1094870/1‑100 (MQ=255)
                               aGTGAAAGGGTATCGGGACAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAt                                                                       >  2:1141318/1‑99 (MQ=255)
                                gTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTa                                                                     >  2:868948/1‑100 (MQ=255)
                                 tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaa                                                                    >  1:1769432/1‑100 (MQ=255)
                                  gAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaaa                                                                   >  1:1690597/1‑100 (MQ=255)
                                       ggTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTaa                                                              >  2:3341149/1‑100 (MQ=255)
                                        gTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAg                                                             <  1:379711/100‑1 (MQ=255)
                                            cGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCggg                                                         <  1:3314837/100‑1 (MQ=255)
                                              ggTCAAAGAAACAAATATTGCACTAGCGCACGCTGCGAAAAGATTGTTAAATATTATTCAAAAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGt                                                        <  2:443627/99‑1 (MQ=255)
                                                tCAAAGAAACAAATATTGCAATACCGCACACTGCGAAAAGATTGTTGAATAAAATTAAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCt                                                     <  2:3165107/100‑1 (MQ=255)
                                                 cAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtt                                                    >  1:3114593/1‑100 (MQ=255)
                                                  aaaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCttt                                                   >  1:3904340/1‑100 (MQ=255)
                                                    aGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTa                                                 >  1:604949/1‑100 (MQ=255)
                                                      aaaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGGATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATc                                               >  2:2260663/1‑100 (MQ=255)
                                                       aaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCa                                              >  1:2869817/1‑100 (MQ=255)
                                                          aaaTATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAt                                           <  1:1067913/100‑1 (MQ=255)
                                                             tatTGCTCTACCGCACACTGCGAAAAGAGTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCTGGTTCTTTTATCAGAt                                           <  2:2853892/97‑1 (MQ=255)
                                                             tatTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTAt                                        >  1:3823243/1‑100 (MQ=255)
                                                              atTACACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTAt                                        <  2:2030941/99‑1 (MQ=255)
                                                                 gCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGtt                                    >  1:3910531/1‑100 (MQ=255)
                                                                   aCTACCGCACACTGCGAAAAGATTGTTGAATATTATTTAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  <  2:855995/100‑1 (MQ=255)
                                                                   aCTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  >  1:2344704/1‑100 (MQ=255)
                                                                      aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTAt                               >  1:730048/1‑100 (MQ=255)
                                                                      aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTAtt                              >  2:1441936/1‑100 (MQ=255)
                                                                      aCCGCACACTGCGAAAAGATTGTTGAATATTATCCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTAt                               >  1:3122059/1‑100 (MQ=255)
                                                                         gCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAATCGGGTTCTTTTATCAGATTATCGTTTCTAtt                              >  2:1417019/1‑98 (MQ=255)
                                                                          cacaCTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATtctc                           >  1:1723323/1‑100 (MQ=255)
                                                                             acTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAAt                        <  2:1176006/100‑1 (MQ=255)
                                                                               tGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTc                      <  2:1616160/100‑1 (MQ=255)
                                                                                gCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCa                     >  1:1190034/1‑100 (MQ=255)
                                                                                  gAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATAAGGCGAATTAAATTTAAGAGGATTCTTTTATCAGATTATCGTTGCTATTATCAATTCATc                   >  2:1994399/1‑100 (MQ=255)
                                                                                   gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  >  1:1965570/7‑100 (MQ=255)
                                                                                   gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  >  1:950006/7‑100 (MQ=255)
                                                                                   gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  >  1:3631776/7‑100 (MQ=255)
                                                                                   aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  <  2:567292/100‑1 (MQ=255)
                                                                                    cctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGa                 >  1:1343303/6‑100 (MQ=255)
                                                                                    aaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGTCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGa                 <  1:167483/100‑1 (MQ=255)
                                                                                     ctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaa                >  1:1226827/5‑100 (MQ=255)
                                                                                       gATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaaaa              >  2:3055383/1‑100 (MQ=255)
                                                                                       gATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaaaa              <  1:3603528/100‑1 (MQ=255)
                                                                                         ttgttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACa            <  2:2111828/100‑1 (MQ=255)
                                                                                         ttgttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACa            <  1:2477250/100‑1 (MQ=255)
                                                                                         ttgttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAACGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAACCa            >  2:114578/1‑100 (MQ=255)
                                                                                          tgttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAg           <  2:3137226/100‑1 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  2:3674515/1‑100 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  1:3744057/1‑100 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  1:150409/1‑100 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          <  1:4129264/100‑1 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  2:495997/1‑100 (MQ=255)
                                                                                            ttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGcc         <  1:2595185/100‑1 (MQ=255)
                                                                                             tgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCca        >  1:2203756/1‑100 (MQ=255)
                                                                                              gAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac       >  1:1924520/1‑100 (MQ=255)
                                                                                              gAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac       <  1:2243078/100‑1 (MQ=255)
                                                                                              gAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac       <  1:4005692/100‑1 (MQ=255)
                                                                                                atatTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaa     >  1:2841546/1‑100 (MQ=255)
                                                                                                  attattCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaaa   >  2:2386093/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAA  >  NC_000913/335067‑335265

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