Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,842,291 G→A C234Y (TGC→TAC)  yicP → putative adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,842,2910GA99.6% 948.0 / ‑5.3 250C234Y (TGC→TAC) yicPputative adenine deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (138/111);  total (138/112)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 252 total aligned reads displayed.

TGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATC  >  NC_000913/3842202‑3842380
                                                                                         |                                                                                         
 gCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCAc                                                                                       >  2:2637665/1‑93 (MQ=255)
  cATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACg                                                                                      <  2:3836588/93‑1 (MQ=255)
   aTTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGa                                                                                     >  2:228930/1‑93 (MQ=255)
   aTTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGa                                                                                     >  2:2843182/1‑93 (MQ=255)
       cGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGt                                                                                 <  2:1913078/93‑1 (MQ=255)
        gCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGtt                                                                                <  1:4136054/93‑1 (MQ=255)
          cACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTAt                                                                              >  2:2055584/1‑93 (MQ=255)
           aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATc                                                                             >  1:3786875/1‑93 (MQ=255)
           aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATc                                                                             <  1:2831006/93‑1 (MQ=255)
            ccTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCa                                                                            >  1:5985087/1‑93 (MQ=255)
               gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCt                                                                         <  1:4436419/93‑1 (MQ=255)
                 cGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTgg                                                                       <  1:4535304/93‑1 (MQ=255)
                 cGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTgg                                                                       <  1:4332384/93‑1 (MQ=255)
                 cGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTgg                                                                       <  2:4224987/93‑1 (MQ=255)
                  gCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGACAGTTATCAGCTgca                                                                      >  1:1600461/1‑91 (MQ=255)
                  gCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGga                                                                      >  2:4977637/1‑93 (MQ=255)
                   cTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGgaa                                                                     >  1:2230312/1‑93 (MQ=255)
                   cTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGgaa                                                                     >  1:1092139/1‑93 (MQ=255)
                   cTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGgaa                                                                     <  2:4613774/93‑1 (MQ=255)
                     ggACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGgaaga                                                                   <  2:1826837/93‑1 (MQ=255)
                      gACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGgaagaa                                                                  <  2:574855/93‑1 (MQ=255)
                      gACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGgaagaa                                                                  <  2:4162577/93‑1 (MQ=255)
                       aCGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGaagaag                                                                 <  2:3602924/93‑1 (MQ=255)
                        cGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAgg                                                                >  2:2569273/1‑93 (MQ=255)
                         ggTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGa                                                               <  2:1582232/93‑1 (MQ=255)
                          gTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGAc                                                              <  2:6358369/93‑1 (MQ=255)
                             aCTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGcc                                                           >  2:3328690/1‑93 (MQ=255)
                             aCTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGcc                                                           <  1:322001/93‑1 (MQ=255)
                               tGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCgg                                                         <  2:2741442/93‑1 (MQ=255)
                               tGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCgg                                                         >  2:341399/1‑93 (MQ=255)
                                  ccGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGaaa                                                      >  2:1164866/1‑93 (MQ=255)
                                  ccGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGaaa                                                      >  2:2131165/1‑93 (MQ=255)
                                   cGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAt                                                     <  1:2026741/93‑1 (MQ=255)
                                    gggTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAtt                                                    >  2:755930/1‑93 (MQ=255)
                                      gTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                  >  2:6148363/1‑93 (MQ=255)
                                      gTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                  >  2:1419169/1‑93 (MQ=255)
                                      gTTTGGGGGGTAAAGAACTTAACGCCTATATTACTGCGGGGATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                  >  1:5659485/1‑93 (MQ=255)
                                       tttGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGCATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACa                                                 >  1:5233791/1‑93 (MQ=255)
                                        ttGGGGGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACaa                                                >  1:1286271/1‑93 (MQ=255)
                                           ggtggtAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTc                                             <  1:3912350/93‑1 (MQ=255)
                                           ggtggtAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTc                                             <  1:5828311/93‑1 (MQ=255)
                                              ggtAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGc                                          >  1:60488/1‑93 (MQ=255)
                                              ggtAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGc                                          >  1:1017681/1‑93 (MQ=255)
                                                tAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCAt                                        >  1:6347760/1‑93 (MQ=255)
                                                tAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCAt                                        >  1:5866446/1‑93 (MQ=255)
                                                 aaaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATg                                       <  1:322879/93‑1 (MQ=255)
                                                 aaaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATg                                       >  1:5528672/1‑93 (MQ=255)
                                                  aaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGt                                      >  1:1863000/1‑93 (MQ=255)
                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAATAAGGACGCCGGAAATTACAACTCGGCATGTc                                     <  1:5148743/93‑1 (MQ=255)
                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTc                                     <  1:6277915/93‑1 (MQ=255)
                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTc                                     >  1:4621612/1‑93 (MQ=255)
                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTc                                     >  1:4861404/1‑93 (MQ=255)
                                                     aaCTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGt                                   >  2:508818/1‑93 (MQ=255)
                                                      aCTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGtt                                  >  1:3018773/1‑93 (MQ=255)
                                                      aCTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGtt                                  >  2:5988849/1‑93 (MQ=255)
                                                       cTTAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtg                                 >  1:1381434/1‑93 (MQ=255)
                                                        ttAACGGCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                <  2:1186217/93‑1 (MQ=255)
                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                <  2:4912952/93‑1 (MQ=255)
                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                <  2:3015021/93‑1 (MQ=255)
                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                >  2:6236085/1‑93 (MQ=255)
                                                            cGCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTgatgat                            <  1:5005362/93‑1 (MQ=255)
                                                             gCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATc                           <  2:4943004/93‑1 (MQ=255)
                                                             gCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATc                           <  1:219242/93‑1 (MQ=255)
                                                             gCCTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATc                           <  1:5296809/93‑1 (MQ=255)
                                                              ccTATATTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATcc                          <  2:5960123/93‑1 (MQ=255)
                                                                 atatTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCcgcg                       >  2:1370268/1‑93 (MQ=255)
                                                                 atatTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCcgcg                       <  2:268438/93‑1 (MQ=255)
                                                                   atTACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGaa                     <  2:4722637/93‑1 (MQ=255)
                                                                    ttACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAg                    <  2:1103789/93‑1 (MQ=255)
                                                                    ttACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAg                    <  1:5607842/93‑1 (MQ=255)
                                                                    ttACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAg                    <  2:3097504/93‑1 (MQ=255)
                                                                     tACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                   <  1:1571575/93‑1 (MQ=255)
                                                                     tACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                   <  2:2675967/93‑1 (MQ=255)
                                                                     tACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                   >  1:3127860/1‑93 (MQ=255)
                                                                     tACTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                   <  1:4774508/93‑1 (MQ=255)
                                                                       cTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                 <  2:3729999/93‑1 (MQ=255)
                                                                       cTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                 <  2:5141278/93‑1 (MQ=255)
                                                                       cTGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                 <  1:3642082/93‑1 (MQ=255)
                                                                        tGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTc                <  1:2680992/93‑1 (MQ=255)
                                                                        tGCGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTc                <  2:5556131/93‑1 (MQ=255)
                                                                          cGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCg              >  1:2851824/1‑93 (MQ=255)
                                                                          cGGGTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCg              <  2:611627/93‑1 (MQ=255)
                                                                           gggTATTGAAACCTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             <  2:230594/93‑1 (MQ=255)
                                                                           gggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             <  2:933774/93‑1 (MQ=255)
                                                                           gggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             <  2:4830459/93‑1 (MQ=255)
                                                                           gggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             >  2:4115549/1‑93 (MQ=255)
                                                                           gggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             <  1:819906/93‑1 (MQ=255)
                                                                           gggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             >  1:2488771/1‑93 (MQ=255)
                                                                           gggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGc             <  1:1191073/93‑1 (MQ=255)
                                                                            ggTATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCt            >  1:3846177/1‑93 (MQ=255)
                                                                              tATTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTTATGATCCGCGAAGGGTCCGCTGc          <  1:5370733/93‑1 (MQ=255)
                                                                               aTTGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGcc         >  2:4423642/1‑93 (MQ=255)
                                                                                ttGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc        >  2:3292325/1‑93 (MQ=255)
                                                                                ttGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc        <  2:1548079/93‑1 (MQ=255)
                                                                                ttGAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc        <  1:1036260/93‑1 (MQ=255)
                                                                                  gAAAACTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGGCCGCTGCCCGc      >  2:762917/1‑93 (MQ=255)
                                                                                   aaaaCTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCa     >  2:1023301/1‑93 (MQ=255)
                                                                                     aaCTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAt   >  1:4837085/1‑93 (MQ=255)
                                                                                     aaCTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAt   >  1:4305830/1‑93 (MQ=255)
                                                                                      aCTACCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAtc  >  2:3377989/1‑93 (MQ=255)
                                                                                         |                                                                                         
TGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATC  >  NC_000913/3842202‑3842380

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: