Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,393 A→T 16.6% intergenic (‑29/‑175) moeA ← / → iaaA molybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913866,3930AT16.6% 401.1 / 70.7 175intergenic (‑29/‑175)moeA/iaaAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (59/87);  new base T (18/11);  total (77/98)
Fisher's exact test for biased strand distribution p-value = 4.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-02

Only 100 of 177 total aligned reads displayed.

AGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATA  >  NC_000913/866294‑866487
                                                                                                   |                                                                                              
  tGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAtt                                                                                              <  2:3352312/100‑1 (MQ=255)
   ggggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAttt                                                                                             >  1:1316500/1‑100 (MQ=255)
    gggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTg                                                                                            <  2:2724647/100‑1 (MQ=255)
      gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGcc                                                                                          <  1:95414/100‑1 (MQ=255)
       tGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGccn                                                                                         <  1:2071056/100‑2 (MQ=255)
        gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtc                                                                                        <  1:3408944/100‑2 (MQ=255)
        gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                        >  1:4386569/1‑100 (MQ=255)
        gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                        >  1:4386459/1‑100 (MQ=255)
         aCGCGAGAACGCATCTCATTAAGCGCCGTCTCGTGCCCCCTCAATCCGGTGGTAACTTCCATGAAAACTTTCCTTTTTATGCAATCTCGAATTTGCCt                                                                                         <  1:1385985/98‑1 (MQ=255)
         aCGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtcc                                                                                       <  2:21110/100‑3 (MQ=255)
         aCGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtcc                                                                                       <  2:225389/100‑3 (MQ=255)
            cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccgca                                                                                    <  2:1400960/100‑6 (MQ=255)
            cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccgca                                                                                    <  2:3391186/100‑6 (MQ=255)
            cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtattat                                                                                    <  2:2432252/100‑1 (MQ=255)
            cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtattat                                                                                    >  1:1021514/1‑100 (MQ=255)
              agaAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGt                                                                                  >  1:4090033/1‑100 (MQ=255)
              agaAAGCATCTCATTAAGCGCAGTGTCTCGCGACATCTATCCGGTGGTAAATGCTATGACAGCCCTCCTGCGCACGCAAAATCGAAATTGCCta                                                                                        <  1:3662891/94‑1 (MQ=255)
               gaAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTc                                                                                 <  2:1346729/100‑1 (MQ=255)
                aaaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                <  2:3467411/100‑1 (MQ=255)
                aaaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                <  1:2184653/100‑1 (MQ=255)
                 aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               >  1:334927/1‑100 (MQ=255)
                 aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  2:1001661/100‑1 (MQ=255)
                 aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  2:2904971/100‑1 (MQ=255)
                 aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  2:1429482/100‑1 (MQ=255)
                  aGCAACTCATTAAGCGCTGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  1:3720584/99‑1 (MQ=255)
                   gCATCTCATTAAGCGCGGTGTCGGGCGACATCAATCCGGTGGTAAATTCCATGAAAACTCTCCTGCGGCGGCAAATTCGATTTTGCCtccg                                                                                      <  1:178147/91‑4 (MQ=255)
                    cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                            <  1:2197458/100‑1 (MQ=255)
                        tcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACg                                                                        >  2:56942/1‑100 (MQ=255)
                        tcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACg                                                                        >  2:3802527/1‑100 (MQ=255)
                            tAAGCGCGGTGTCGAGCGACATCAATCCGGTAGTAAATTCCATGAAAACCCTCCTGCGTAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCcac                                                                    <  1:2402138/100‑1 (MQ=255)
                              aGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAg                                                                  <  1:2942445/100‑1 (MQ=255)
                               gcgcGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGa                                                                 >  2:2786621/1‑100 (MQ=255)
                                cgcgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGAc                                                                <  1:10274/100‑1 (MQ=255)
                                  cgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTg                                                              <  1:3481238/100‑1 (MQ=255)
                                  cgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTg                                                              <  2:3750182/100‑1 (MQ=255)
                                    gtgtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACCGACTGTa                                                            <  2:3286775/100‑1 (MQ=255)
                                    gtgtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTa                                                            <  1:404081/100‑1 (MQ=255)
                                    gtgtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAAACGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTa                                                            >  2:677510/1‑100 (MQ=255)
                                    gtgtCGAGCAACAGAAAACCAGTGGTAAATTCCATGAAAACACTCATGCGGAGGCAAAATCGAATTTGCATAATATATAAAAAAAACGCCACAGACTGTa                                                            <  1:3525746/100‑1 (MQ=255)
                                      gtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATg                                                          >  2:3603626/1‑100 (MQ=255)
                                        cGATCGACATCAAATCGGTGGTAAATTCCATGCAACCCCTCCTCCGGAGGCAAAATCGAATTTTCCTATTATGTCAGAAAAACGCCACAGACTGTATg                                                          <  2:3982960/98‑1 (MQ=255)
                                        cGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcc                                                        >  1:3057604/1‑100 (MQ=255)
                                         gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGTCa                                                       >  1:2515701/1‑100 (MQ=255)
                                         gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                       >  2:1319580/1‑100 (MQ=255)
                                         gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                       >  2:766375/1‑100 (MQ=255)
                                            cGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCt                                                    >  1:41970/1‑100 (MQ=255)
                                             gACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTc                                                   >  2:87260/1‑100 (MQ=255)
                                             gACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTc                                                   <  2:1716106/100‑1 (MQ=255)
                                               cATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCgg                                                 >  2:4175587/1‑100 (MQ=255)
                                                aTCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCggg                                                <  1:61979/100‑1 (MQ=255)
                                                       cGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAgcg                                         >  2:367630/1‑100 (MQ=255)
                                                          tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTg                                      <  1:796114/100‑1 (MQ=255)
                                                          tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTg                                      <  1:1597554/100‑1 (MQ=255)
                                                           ggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTgg                                     <  2:2625140/100‑1 (MQ=255)
                                                             tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGt                                   <  1:3897120/100‑1 (MQ=255)
                                                               aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                 <  2:963705/100‑1 (MQ=255)
                                                                aTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                <  1:1372181/100‑1 (MQ=255)
                                                                aTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                >  2:4267779/1‑100 (MQ=255)
                                                                aTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACACCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                >  2:702940/1‑100 (MQ=255)
                                                                aTTCCATGAAAACACGCCTGCGGTGGCAAAATCGGATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACTTCCGGCGTAGCGCTGGGTCCt                                >  1:2525496/1‑100 (MQ=255)
                                                                 ttCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTg                               <  1:688021/100‑1 (MQ=255)
                                                                 ttCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTg                               <  2:1891/100‑1 (MQ=255)
                                                                  tCCAGGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACTGACTGTATGCCACCTCGGTCGTAGCGCTGGGTCCTGc                              <  1:1209817/100‑1 (MQ=255)
                                                                   ccATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGcc                             >  2:2351373/1‑100 (MQ=255)
                                                                    cATGAAAACACTCCTGCGGAGGCAAACTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGACGTAGCGCTGGGTCCTGc                              >  1:2695424/1‑98 (MQ=255)
                                                                    cATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCt                            >  2:1587132/1‑100 (MQ=255)
                                                                    cATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCt                            >  2:2381147/1‑100 (MQ=255)
                                                                    cATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCt                            >  2:3513572/1‑100 (MQ=255)
                                                                     aTGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCtt                           >  2:1508480/1‑100 (MQ=255)
                                                                       gAAAACGCTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                         <  2:935740/100‑1 (MQ=255)
                                                                       gAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                         <  2:1761199/100‑1 (MQ=255)
                                                                        aaaaCTCTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                         <  1:3350450/99‑1 (MQ=255)
                                                                        aaaaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTAc                        >  1:3108200/1‑100 (MQ=255)
                                                                         aaaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                       <  1:1955528/100‑1 (MQ=255)
                                                                         aaaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAAACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                       >  1:2814981/1‑100 (MQ=255)
                                                                           acacTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                     <  2:2393663/100‑1 (MQ=255)
                                                                           acacTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                     <  2:1910575/100‑1 (MQ=255)
                                                                            cacTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGc                    <  2:2650051/100‑1 (MQ=255)
                                                                             acTCCTGCGGAGGCACAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGcc                   <  2:4337589/100‑1 (MQ=255)
                                                                             acTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGcc                   <  1:1391511/100‑1 (MQ=255)
                                                                              tgacataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                  >  1:3757968/10‑100 (MQ=38)
                                                                              tgacataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                  >  1:1854420/10‑100 (MQ=38)
                                                                              cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                  >  2:1617968/1‑100 (MQ=255)
                                                                                 cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTTGGTCCTGCCTTTACATGCCatat               >  1:339273/7‑100 (MQ=38)
                                                                                 cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  1:2608293/7‑100 (MQ=39)
                                                                                 cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  2:1632980/7‑100 (MQ=39)
                                                                                 cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  1:2608368/7‑100 (MQ=39)
                                                                                 cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACATCGGGCGTAGCACTGGGTCCTGCCTTTACATGCCaca                >  1:2413053/7‑97 (MQ=38)
                                                                                 cataataGGCAAACTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  1:2268376/7‑100 (MQ=255)
                                                                                   taataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc             >  2:3403867/5‑100 (MQ=255)
                                                                                   taataGGCAAATACGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc             >  2:2893197/5‑100 (MQ=255)
                                                                                   gCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc             <  2:2134268/100‑1 (MQ=255)
                                                                                    cGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCa            >  2:754706/1‑100 (MQ=255)
                                                                                      gcggCAAAATCGAATTTGCCTATTATGGCAGAAAAACGCCACAGACTGTAGGCCACCTCTGGCGTAGCGCTGGGTCCTGCCTTTACATGCCAGATCCa            <  1:1787627/96‑1 (MQ=255)
                                                                                       aGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCt         >  2:607357/1‑100 (MQ=255)
                                                                                          cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc      >  2:3100092/1‑100 (MQ=255)
                                                                                          cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc      >  1:4031690/1‑100 (MQ=255)
                                                                                           aaaaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCt     >  2:40104/1‑100 (MQ=255)
                                                                                              ntCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCtata  >  1:3992403/2‑100 (MQ=255)
                                                                                              aTCGAATTTGCCTATTATGTCAGAATGTCGCCACTGACTGTTTGCCCCCTCGGGCGTATCGATGGGTCCGGCCTTTACATGCCATAGCCATCTTTCtat   <  1:2362149/99‑1 (MQ=255)
                                                                                                   |                                                                                              
AGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATA  >  NC_000913/866294‑866487

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: