Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 335,166 T→A 15.1% 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,1660TA15.1% 547.2 / 73.5 211intergenic (‑144/‑114)yahC/yahDputative inner membrane protein/ankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base T (105/74);  new base A (12/20);  total (117/94)
Fisher's exact test for biased strand distribution p-value = 3.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 228 total aligned reads displayed.

GTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAA  >  NC_000913/335067‑335265
                                                                                                   |                                                                                                   
  gCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatta                                                                                                   <  1:1540570/100‑1 (MQ=255)
  gCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatta                                                                                                   >  2:1193533/1‑100 (MQ=255)
  gCAGACTATAGAGAAATAAAACCTGAGTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAataata                                                                                                   <  2:83851/100‑1 (MQ=255)
       cTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaa                                                                                              >  2:3888227/1‑100 (MQ=255)
        tataAAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaac                                                                                             <  1:1644931/100‑1 (MQ=255)
         ataGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTcaa                                                                                              >  2:2532717/1‑98 (MQ=255)
          taGAGAAATAAAACCTGATTTAGTGAAAGCGTATCGAGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGAGTGTTGAATATTATTcaacaa                                                                                           <  2:3679453/100‑1 (MQ=255)
          taGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAGAAGATTGTTGAATAATATTcaacaa                                                                                           <  2:1269772/100‑1 (MQ=255)
           agagAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaa                                                                                          >  2:953284/1‑100 (MQ=255)
           agagAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaat                                                                                          <  2:2390164/100‑2 (MQ=255)
           agagAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaat                                                                                          <  1:2451115/100‑2 (MQ=255)
           agagAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaat                                                                                          <  1:4057725/100‑2 (MQ=255)
           agagAAATAAAACCTGAGTGAGTGAAAGGGTATCGTGTCAAAGAAACAAATATTGCTCTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaat                                                                                          <  2:2070195/100‑2 (MQ=255)
            gagaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaa                                                                                         >  2:4174805/1‑100 (MQ=255)
              gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAg                                                                                       >  1:2791812/1‑100 (MQ=255)
              gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatctt                                                                                       <  2:1854118/100‑5 (MQ=255)
                aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGc                                                                                     >  2:1993907/1‑100 (MQ=255)
                aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGc                                                                                     >  2:1540880/1‑100 (MQ=255)
                aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGc                                                                                     >  1:3523605/1‑100 (MQ=255)
                aaTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTcaacaatctttt                                                                                     <  2:608508/100‑7 (MQ=255)
                 aTAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCg                                                                                    >  2:1062778/1‑100 (MQ=255)
                  tAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCga                                                                                   >  2:3771100/1‑100 (MQ=255)
                   aaaaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCgag                                                                                  <  2:1633618/100‑1 (MQ=255)
                     aaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGAt                                                                                >  1:2663515/1‑100 (MQ=255)
                      aCCTGAGTGTGTGTAAGGGTCTCGCGTGTTAGAAACAAATATTGGACAACGGCACCGGGCGTAAAGATTGTTAAAAATTATTcaac                                                                                             >  2:671321/1‑86 (MQ=255)
                       ccTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATAc                                                                              >  1:622789/1‑100 (MQ=255)
                       ccTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATAc                                                                              >  2:112217/1‑100 (MQ=255)
                          gATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGc                                                                           <  1:142133/100‑1 (MQ=255)
                          gATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGc                                                                           <  2:2239648/100‑1 (MQ=255)
                           aTTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCg                                                                          >  1:1994752/1‑100 (MQ=255)
                            ttGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGa                                                                         >  1:3858090/1‑100 (MQ=255)
                             tGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGaa                                                                        >  2:3020360/1‑100 (MQ=255)
                               aGTGAAAGGGTATCGGGTCAAAGAAAGAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAtt                                                                      >  1:2534248/1‑100 (MQ=255)
                               aGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAAtt                                                                      >  1:3937528/1‑100 (MQ=255)
                                 tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaa                                                                    >  2:874248/1‑100 (MQ=255)
                                  gAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaaa                                                                   <  1:3451406/100‑1 (MQ=255)
                                   aaaGGGTATCGGGTCAAAGAAACAAATATTGCGCTACCGAACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  <  2:2927735/100‑1 (MQ=255)
                                   aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  >  2:934019/1‑100 (MQ=255)
                                   aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAAt                                                                  >  2:1662045/1‑100 (MQ=255)
                                      gggTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTa                                                               <  1:2606493/100‑1 (MQ=255)
                                       ggTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTaa                                                              >  1:2051793/1‑100 (MQ=255)
                                          aTCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCg                                                           <  1:2353538/100‑1 (MQ=255)
                                          aTCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCg                                                           <  1:1107383/100‑1 (MQ=255)
                                          aTCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCg                                                           >  1:3175451/1‑100 (MQ=255)
                                             gggTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCgg                                                          >  2:3555633/1‑98 (MQ=255)
                                               gTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTc                                                      >  1:1065232/1‑100 (MQ=255)
                                                 cAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtt                                                    >  2:3201874/1‑100 (MQ=255)
                                                 cAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtt                                                    >  1:227195/1‑100 (MQ=255)
                                                  aaaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCttt                                                   <  2:2201914/100‑1 (MQ=255)
                                                  aaaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCttt                                                   <  2:1228596/100‑1 (MQ=255)
                                                   aaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                  >  1:1319079/1‑100 (MQ=255)
                                                   aaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                  >  2:1340118/1‑100 (MQ=255)
                                                   aaGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCttt                                                   <  1:3886067/99‑1 (MQ=255)
                                                      aaaCAAATATTGCACTACCGCACACTTCGAAAAGATTGTTGATTATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATc                                               <  2:3379468/100‑1 (MQ=255)
                                                      aaaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATc                                               <  1:4122125/100‑1 (MQ=255)
                                                      aaaCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATc                                               <  2:2656251/100‑1 (MQ=255)
                                                       aaCAAATATTGCACTACCGCACACTGCGAAAAGCTTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGATCTTTTATCa                                              >  2:2734119/1‑100 (MQ=255)
                                                        aCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAg                                             >  1:931277/1‑100 (MQ=255)
                                                        aCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAg                                             <  1:433058/100‑1 (MQ=255)
                                                         cAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGa                                            >  2:288288/1‑100 (MQ=255)
                                                         cAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGAGATCTTTTATCAGa                                            >  1:3903796/1‑100 (MQ=255)
                                                          aaaTATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAt                                           <  1:1205841/100‑1 (MQ=255)
                                                          aaaTATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAt                                           >  1:3575545/1‑100 (MQ=255)
                                                            atatTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAt                                           >  2:715178/1‑98 (MQ=255)
                                                              atTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATc                                       >  1:2512532/1‑100 (MQ=255)
                                                              atTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATc                                       >  2:2308701/1‑100 (MQ=255)
                                                               ttGCACTACCGCACACTGCGAAAAGATTGTGGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCg                                      <  1:3497912/100‑1 (MQ=255)
                                                                tGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGACGAACTAAATTTAAGCGAGTTCTTTTATCAGATTATCGt                                     >  2:1504332/1‑100 (MQ=255)
                                                                 gCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCg                                      <  2:1562751/98‑1 (MQ=255)
                                                                 gCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGtt                                    >  1:1017813/1‑100 (MQ=255)
                                                                   aCTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  <  1:2991931/100‑1 (MQ=255)
                                                                   aCTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  >  1:1701705/1‑100 (MQ=255)
                                                                   aCTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  <  1:3079169/100‑1 (MQ=255)
                                                                   aCTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  <  2:2065044/100‑1 (MQ=255)
                                                                    cTACCGCACACTGCGAAAAGATTGTTGAAGATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                  <  2:3425078/99‑1 (MQ=255)
                                                                      aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTAt                               >  2:3166792/1‑100 (MQ=255)
                                                                         gCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATtct                            <  2:1404306/100‑1 (MQ=255)
                                                                            cacTTCGAAAAGATTGGTTAATATTATTCAACAAGAATGCGTGATACGGCGAATTAAATTTAAGCGGGTTCGTTTATCAGATTATCGTTGCTATTCTCa                          <  2:3255935/99‑1 (MQ=255)
                                                                            cacTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCaa                         <  1:2393164/100‑1 (MQ=255)
                                                                              cTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAAtt                       <  1:3835330/100‑1 (MQ=255)
                                                                               tGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGTTTCTTTTATCAGATTATCGTTGCTATTCTCAATTc                      <  1:3982927/100‑1 (MQ=255)
                                                                                gCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCa                     >  1:856592/1‑100 (MQ=255)
                                                                                gCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCa                     >  1:1776935/1‑100 (MQ=255)
                                                                                 cGAAAAGATTGTTGAATATTATTCATCAAAAAGCCGAGATACGGCGAATTAAATTTAAGCGGGTTCTGTGACCAGATGATCGTTGCTATTCTCAAt                        <  2:486538/96‑1 (MQ=255)
                                                                                  gAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGCTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATc                   >  2:3307476/1‑100 (MQ=255)
                                                                                  gAAAAGATTGTTGAATAAAAAAAAAAAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATc                   <  2:1745124/100‑1 (MQ=255)
                                                                                  cgcctttttgttaGATAATATTCAGCTAAACGGCGGGCTACGGGGAATTAAATTTAATCGGGTTCTTTAGTCAGATTATCGTTGCTATTCACAATTCagc                   >  2:3293284/8‑98 (MQ=255)
                                                                                   gcctttttgttgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  >  1:133604/7‑100 (MQ=255)
                                                                                   aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  <  1:3771128/100‑1 (MQ=255)
                                                                                   aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCg                  <  1:2865829/100‑1 (MQ=255)
                                                                                    aaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGa                 <  2:1725105/100‑1 (MQ=255)
                                                                                     aaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaa                >  1:84974/1‑100 (MQ=255)
                                                                                       gATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaaaa              >  2:3362712/1‑100 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  1:3246094/1‑100 (MQ=255)
                                                                                           gttgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGc          >  1:1619899/1‑100 (MQ=255)
                                                                                             tgAATAATATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCca        >  1:337507/1‑100 (MQ=255)
                                                                                               aaTATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCcac       <  2:2585282/99‑1 (MQ=255)
                                                                                                   ttattCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaaaa  >  1:686456/1‑100 (MQ=255)
                                                                                                   ttattCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaaaa  <  1:1813800/100‑1 (MQ=255)
                                                                                                   ttattCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACaaaaa  >  1:3756385/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: