Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,749 C→A 100% D578Y (GAT→TAT)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,7490CA99.4% 510.7 / ‑4.9 152D578Y (GAT→TAT) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (78/73);  total (79/73)
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.47e-01

Only 100 of 152 total aligned reads displayed.

GAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATG  >  NC_000913/843653‑843844
                                                                                                |                                                                                               
gAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATcc                                                                                              <  1:3478117/100‑1 (MQ=255)
 aaTGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCa                                                                                             <  1:3883608/100‑1 (MQ=255)
  atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                            <  2:3447791/100‑1 (MQ=255)
  atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                            <  2:1306273/100‑1 (MQ=255)
  atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                            <  1:1841317/100‑1 (MQ=255)
      tgtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGtt                                                                                        >  2:3725897/1‑100 (MQ=255)
      tgggagaggCCCCGTTATGCTGTCGTACGCCCGCGGAGCGAAACGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtg                                                                                          <  2:3660722/90‑1 (MQ=255)
         gATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAt                                                                                     >  1:1315428/1‑100 (MQ=255)
          aTACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATg                                                                                    <  2:3088231/100‑1 (MQ=255)
             cGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCg                                                                                 >  2:685907/1‑100 (MQ=255)
               ccccGGTATCCTGTCGTACGCCCGCGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtt                                                                               >  2:3288628/1‑100 (MQ=255)
               ccccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCCGTGTTCATGCCGtt                                                                               >  1:4106593/1‑100 (MQ=255)
               ccccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtt                                                                               <  1:143236/100‑1 (MQ=255)
                cccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCg                                                                                 <  2:1326007/97‑1 (MQ=255)
                cccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGttt                                                                              >  1:3988326/1‑100 (MQ=255)
                  cGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTc                                                                            <  1:1369113/100‑1 (MQ=255)
                  cGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTc                                                                            <  1:1193868/100‑1 (MQ=255)
                     tATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCaaa                                                                         <  2:333532/100‑1 (MQ=255)
                      aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                        <  2:2704988/100‑1 (MQ=255)
                      aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAAc                                                                        <  1:1605450/100‑1 (MQ=255)
                        ccTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTg                                                                      <  2:1988756/100‑1 (MQ=255)
                          tGTCGTACGCCCACGGAGCGAATCGACATCCGTTGCACTGTGCCGGTCAACGGCCCGATAGTAACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaa                                                                    >  2:216947/1‑100 (MQ=255)
                            tCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaata                                                                  <  1:2302194/100‑1 (MQ=255)
                              gTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaa                                                                >  2:4148569/1‑100 (MQ=255)
                              gTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaa                                                                <  1:3127523/100‑1 (MQ=255)
                              gTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaa                                                                >  1:3202319/1‑100 (MQ=255)
                               tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                               <  1:2815480/100‑1 (MQ=255)
                               tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                               <  2:2356164/100‑1 (MQ=255)
                                 cGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaata                                                              >  2:2886472/1‑99 (MQ=255)
                                   cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAccc                                                           >  1:704756/1‑100 (MQ=255)
                                   cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAccc                                                           >  1:1957408/1‑100 (MQ=255)
                                    ccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAcccc                                                          >  2:1019884/1‑100 (MQ=255)
                                      aCGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCgg                                                        <  1:1843184/100‑1 (MQ=255)
                                       cGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTGCAAACTGAATAAATACCCCGGt                                                       >  2:3924556/1‑100 (MQ=255)
                                        ggAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTg                                                      <  1:3771130/100‑1 (MQ=255)
                                           gCGAATCCACATCCGTTCCACTGTGCTGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCTGTTTTCAAACTGAATAACTACCCCGGt                                                       <  2:2980940/96‑1 (MQ=255)
                                              aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                >  1:653919/1‑100 (MQ=255)
                                              aaTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCAGGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                <  2:2239185/100‑1 (MQ=255)
                                               aTCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAatat                                               <  1:1286940/100‑1 (MQ=255)
                                                tCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATc                                              <  1:2078358/100‑1 (MQ=255)
                                                tCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATc                                              <  1:685908/100‑1 (MQ=255)
                                                 cGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCt                                             >  2:3925102/1‑100 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                            >  1:2322967/1‑100 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                            >  1:317034/1‑100 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                            >  1:3150879/1‑100 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                            <  1:3236680/100‑1 (MQ=255)
                                                  gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATAGCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCtt                                            <  2:1614510/100‑1 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                          <  1:1964568/100‑1 (MQ=255)
                                                    cATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTc                                          >  1:1191118/1‑100 (MQ=255)
                                                     aTCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCa                                         <  2:2577693/100‑1 (MQ=255)
                                                     aTCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCa                                         <  2:4060265/100‑1 (MQ=255)
                                                     aTCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCACTGTTCATGCCGTTTTCAAACTGAATAAATCCCCCGGTGATAATACCttt                                           >  2:547760/1‑98 (MQ=255)
                                                      tCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                        <  1:3632156/100‑1 (MQ=255)
                                                       ccGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAcc                                       >  2:3134628/1‑100 (MQ=255)
                                                        cGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCa                                      >  1:685820/1‑100 (MQ=255)
                                                        cGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTAAAACTGAATAAATAGCCCGGTGataatggatttacg                                        >  1:3782375/1‑89 (MQ=255)
                                                         gTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAg                                     >  1:1973126/1‑100 (MQ=255)
                                                         gTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAg                                     <  2:2007708/100‑1 (MQ=255)
                                                              aCTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGttt                                <  2:3192109/100‑1 (MQ=255)
                                                                tgtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGc                              <  2:3176359/100‑1 (MQ=255)
                                                                 gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                             >  2:1363812/1‑100 (MQ=255)
                                                                 gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                             >  1:2884936/1‑100 (MQ=255)
                                                                 gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCg                             >  2:265418/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                           <  2:1360401/100‑1 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                           >  1:703471/1‑100 (MQ=255)
                                                                   gCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                           >  2:2642162/1‑100 (MQ=255)
                                                                    ccGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAc                          >  1:928597/1‑100 (MQ=255)
                                                                    ccGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAc                          >  2:2323129/1‑100 (MQ=255)
                                                                    ccGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAc                          <  1:3725918/100‑1 (MQ=255)
                                                                     cGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAcc                         >  2:1959493/1‑100 (MQ=255)
                                                                     cGGTCAACGGCCCGATAATCACCAAAGATCCAGTGTTCATGCCGTTTTCAAACTGAATAAACACCCCGGTTATAATATCTTTCACCAGCGTTTGCGAAcc                         <  2:1418734/100‑1 (MQ=255)
                                                                      ggTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAAcc                         <  1:1070451/99‑1 (MQ=255)
                                                                       gTCAACGGCCCGATAGTCACCAAATATCCCGTGTTGATGCCGTTTTCAAACTGAATAAAGACCCCGGTGATAAGATCTGTCACCAGCGTTTGCGa                            <  2:3479314/95‑1 (MQ=255)
                                                                        tCAACGGCCCGATAGTCACCAAATATCCAGGGTACATGCCGTTTTCAAACTGCATAAATACCCCGGTGATGATATCTTTCACCAGCGTTTGCGAACCaaa                      >  2:2512098/1‑100 (MQ=255)
                                                                         cAACGGCCCGATCGTCACCAAATATCCAGTGTTCATGCCGCTTTCAAACTGAATAAATACCCCAGTGATAATATCTTTCACCAGCGTTTGCGAACCaaa                      >  2:932023/1‑99 (MQ=255)
                                                                         cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                     >  1:233826/1‑100 (MQ=255)
                                                                         cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                     <  1:845567/100‑1 (MQ=255)
                                                                         cAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACACCGGTGATAATATATTTCACCAGAGTTTGAGAACCaaa                      >  2:3396446/1‑99 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAc                     >  2:1912006/1‑99 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                    >  1:1692368/1‑100 (MQ=255)
                                                                          aaCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACg                    <  1:621672/100‑1 (MQ=255)
                                                                           aCGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGGATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACga                   >  2:3158033/1‑100 (MQ=255)
                                                                           aCGGCCCGATAGTCACCAAATATCCAGTGTACACGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCAAACCAAACga                   >  2:2751108/1‑100 (MQ=255)
                                                                               cccGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATa               >  2:1214482/1‑100 (MQ=255)
                                                                                 cGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGc             <  1:3090512/100‑1 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc            <  1:2026143/100‑1 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc            <  1:2276248/100‑1 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc            <  2:3723797/100‑1 (MQ=255)
                                                                                  gATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc            <  2:2665098/100‑1 (MQ=255)
                                                                                     aGTCACCAAATATCCAGTTTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGa         <  2:1207376/100‑1 (MQ=255)
                                                                                     aGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGa         <  2:80337/100‑1 (MQ=255)
                                                                                      gTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAc        >  2:534399/1‑100 (MQ=255)
                                                                                       tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc       <  1:1903724/100‑1 (MQ=255)
                                                                                       tCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc       >  2:1808487/1‑100 (MQ=255)
                                                                                        cACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAAAGAGATAGCCAGACCt      >  1:3908766/1‑100 (MQ=255)
                                                                                        cACCAAATATCCAGTGTTCATGCCGTGTTCAAACTGAATAAATACCCAGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt      >  2:3018260/1‑100 (MQ=255)
                                                                                         aCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTa     >  2:1848464/1‑100 (MQ=255)
                                                                                           cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt   <  1:3008757/100‑1 (MQ=255)
                                                                                           cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt   <  1:3452194/100‑1 (MQ=255)
                                                                                           cAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAAt   >  1:1035193/1‑100 (MQ=255)
                                                                                                |                                                                                               
GAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATG  >  NC_000913/843653‑843844

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: