Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,838,233 T→C T433T (ACA→ACG yicM ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,2330TC99.2% 810.4 / ‑3.2 252T433T (ACA→ACGyicMputative transport protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (104/146);  total (106/146)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 268 total aligned reads displayed.

AACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAAC  >  NC_000913/3838143‑3838325
                                                                                          |                                                                                            
aacaacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAcgt                                                                                            >  1:1519429/1‑93 (MQ=255)
aacaacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAcgt                                                                                            <  2:4580829/93‑1 (MQ=255)
aacaacTACTGACCAAAGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAcgt                                                                                            <  1:4195425/93‑1 (MQ=255)
  caacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGc                                                                                          <  1:4327515/93‑1 (MQ=255)
  caacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGc                                                                                          <  2:3708442/93‑1 (MQ=255)
   aacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGcc                                                                                         <  2:2773841/93‑1 (MQ=255)
   aacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGcc                                                                                         >  2:2868649/1‑93 (MQ=255)
    actactGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCg                                                                                        <  2:4258335/93‑1 (MQ=255)
     ctactGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCgg                                                                                       >  1:4041072/1‑93 (MQ=255)
         tGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGacaa                                                                                   >  1:1262627/1‑93 (MQ=255)
         tGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGacaa                                                                                   <  2:691876/93‑1 (MQ=255)
          gACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGacaac                                                                                  >  1:4663143/1‑93 (MQ=255)
           aCCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGacaaca                                                                                 >  1:184997/1‑93 (MQ=255)
            ccAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACAt                                                                                >  1:2036310/1‑93 (MQ=255)
            ccAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACAt                                                                                <  2:3638879/93‑1 (MQ=255)
             cAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATc                                                                               <  1:3209981/93‑1 (MQ=255)
              aaCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCa                                                                              >  2:2801200/1‑93 (MQ=255)
              aaCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCa                                                                              >  1:3438505/1‑93 (MQ=255)
               aCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCaa                                                                             <  1:4466860/93‑1 (MQ=255)
                cGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAAc                                                                            <  1:2453492/93‑1 (MQ=255)
                 gCTTCTTTACCTGTTCTCAGGATTTCTTCATTTGCGCCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACg                                                                           <  1:4213774/93‑1 (MQ=255)
                 gCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACg                                                                           <  2:4087535/93‑1 (MQ=255)
                 gCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACg                                                                           >  1:2305374/1‑93 (MQ=255)
                 gCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACg                                                                           >  2:2986703/1‑93 (MQ=255)
                 gCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACg                                                                           <  2:4699150/93‑1 (MQ=255)
                 gCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACg                                                                           >  1:655687/1‑93 (MQ=255)
                  cttcttTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACgg                                                                          >  2:1406357/1‑93 (MQ=255)
                    tcttTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGACAGCAACATCAACGTGGCGGACAACATCAACGGCg                                                                        >  1:1897743/1‑93 (MQ=255)
                     cttTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGa                                                                       <  1:1731462/93‑1 (MQ=255)
                      tttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGaa                                                                      <  1:3563469/93‑1 (MQ=255)
                      tttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGaa                                                                      <  2:717516/93‑1 (MQ=255)
                       ttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAg                                                                     <  1:1070835/93‑1 (MQ=255)
                           cTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAg                                                                 >  1:4115292/1‑93 (MQ=255)
                            tGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGa                                                                >  1:1637868/1‑93 (MQ=255)
                              atatCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGAcc                                                              <  1:1017349/93‑1 (MQ=255)
                                atCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCaa                                                            >  1:1043551/1‑93 (MQ=255)
                                 tCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGATGTCAGACCAat                                                           >  1:4313604/1‑93 (MQ=255)
                                 tCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAat                                                           <  1:931114/93‑1 (MQ=255)
                                    ggATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAAtatt                                                        <  1:2114800/93‑1 (MQ=255)
                                    ggATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAAtatt                                                        >  2:4691517/1‑93 (MQ=255)
                                     gATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAAtatta                                                       >  2:2463405/1‑93 (MQ=255)
                                     gATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAAtatta                                                       <  2:1023303/93‑1 (MQ=255)
                                        ttcttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATcc                                                    <  2:1282907/93‑1 (MQ=255)
                                         tcttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCa                                                   <  2:2323497/93‑1 (MQ=255)
                                          cttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAg                                                  <  1:4733563/93‑1 (MQ=255)
                                           ttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAgc                                                 >  1:1762646/1‑93 (MQ=255)
                                             cATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAgcgc                                               >  2:2397119/1‑93 (MQ=255)
                                                tttCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATa                                            >  1:2800248/1‑93 (MQ=255)
                                                tttCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATa                                            <  2:4315246/93‑1 (MQ=255)
                                                tttCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATa                                            <  1:2709702/93‑1 (MQ=255)
                                                 ttCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATaa                                           <  1:834862/93‑1 (MQ=255)
                                                  tCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAc                                          <  2:3471445/93‑1 (MQ=255)
                                                  tCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAc                                          <  2:2262064/93‑1 (MQ=255)
                                                  tCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAc                                          <  2:2076968/93‑1 (MQ=255)
                                                  tCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAc                                          <  1:1323784/93‑1 (MQ=255)
                                                  tCACCTTTGCAGTAACCAACAAGGCAGGCAGCAACATCAACGTGCCGGACAACATCACCGGCGAAGTCAGACCAATATTATCCAGCGCATAAc                                          <  1:4288957/93‑1 (MQ=255)
                                                   cACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAACGCATAAcc                                         <  2:3145126/93‑1 (MQ=255)
                                                     ccTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAccgc                                       <  2:3407736/93‑1 (MQ=255)
                                                     ccTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAccgc                                       >  1:2667095/1‑93 (MQ=255)
                                                      cTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAccgcc                                      >  2:4150187/1‑93 (MQ=255)
                                                      cTTTGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAccgcc                                      <  1:2000870/93‑1 (MQ=255)
                                                       tttGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAccgccg                                     <  1:2403458/93‑1 (MQ=255)
                                                       tttGCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAAccgccg                                     <  1:451386/93‑1 (MQ=255)
                                                          gCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGAtt                                  >  1:915147/1‑93 (MQ=255)
                                                          gCAGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGAtt                                  <  2:2649277/93‑1 (MQ=255)
                                                            aGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGc                                >  2:1943450/1‑93 (MQ=255)
                                                            aGTAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGc                                >  1:60390/1‑93 (MQ=255)
                                                              tAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCcg                              <  1:531986/93‑1 (MQ=255)
                                                              tAACCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCcg                              <  2:94914/93‑1 (MQ=255)
                                                               aaCCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCcgc                             <  1:724240/93‑1 (MQ=255)
                                                               aaCCAACAATGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCcgc                             >  2:8425/1‑93 (MQ=255)
                                                                  caacaaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGcc                          >  1:2834705/1‑93 (MQ=255)
                                                                  caacaaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGcc                          <  2:2545042/93‑1 (MQ=255)
                                                                    acaaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCcac                        <  2:437966/93‑1 (MQ=255)
                                                                    acaaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCcac                        <  1:3385979/93‑1 (MQ=255)
                                                                     caaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCcaca                       <  2:308108/93‑1 (MQ=255)
                                                                     caaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCcaca                       <  2:2634591/93‑1 (MQ=255)
                                                                     caaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCcaca                       <  1:1586919/93‑1 (MQ=255)
                                                                      aaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGGCCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAg                      <  1:142717/93‑1 (MQ=255)
                                                                      aaTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAg                      >  2:3428832/1‑93 (MQ=255)
                                                                       aTGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAgg                     >  1:2944778/1‑93 (MQ=255)
                                                                        tGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGt                    >  2:2287867/1‑93 (MQ=255)
                                                                        tGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGt                    >  2:3558793/1‑93 (MQ=255)
                                                                        tGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGt                    >  2:3591970/1‑93 (MQ=255)
                                                                        tGCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGt                    <  2:644342/93‑1 (MQ=255)
                                                                         gCAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTa                   <  1:1895273/93‑1 (MQ=255)
                                                                          cAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTAt                  >  1:3351929/1‑93 (MQ=255)
                                                                          cAGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTAt                  <  1:2957314/93‑1 (MQ=255)
                                                                           aGTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTAtt                 <  2:291748/93‑1 (MQ=255)
                                                                            gTCAGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTa                <  1:4710307/93‑1 (MQ=255)
                                                                               aGCAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCa             <  2:3947878/93‑1 (MQ=255)
                                                                                 cAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAg           <  1:1969579/93‑1 (MQ=255)
                                                                                 cAACATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAg           <  2:2015527/93‑1 (MQ=255)
                                                                                   aCATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCt         >  1:3130564/1‑93 (MQ=255)
                                                                                   aCATCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCt         >  1:4241630/1‑93 (MQ=255)
                                                                                     aTCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGa       >  1:4530424/1‑93 (MQ=255)
                                                                                      tCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGaa      <  2:1622798/93‑1 (MQ=255)
                                                                                      tCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGaa      >  1:3850518/1‑93 (MQ=255)
                                                                                      tCAACGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGAGCCACAGGTATTAGCAAGCTGaa      >  2:3381089/1‑93 (MQ=255)
                                                                                          tgtgCCGGACAATATCATCGGCTAAGTAAGACAAATATTATCTAGGGCAAATCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAAc  <  2:4478962/93‑1 (MQ=255)
                                                                                          |                                                                                            
AACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAAC  >  NC_000913/3838143‑3838325

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