Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,388 A→T 15.8% intergenic (‑24/‑180) 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,3880AT15.8% 413.0 / 66.7 175intergenic (‑24/‑180)moeA/iaaAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (58/86);  new base T (17/10);  total (77/98)
Fisher's exact test for biased strand distribution p-value = 3.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-02

Only 100 of 181 total aligned reads displayed.

GTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATA  >  NC_000913/866291‑866487
                                                                                                 |                                                                                                   
gTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATc                                                                                                   >  2:2479101/1‑100 (MQ=255)
 tCAGTGGGGTGACGCGGGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGACACAATCg                                                                                                  >  1:2116996/1‑100 (MQ=255)
 tCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCg                                                                                                  >  2:1456831/1‑100 (MQ=255)
  cAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCg                                                                                                  <  1:562372/99‑1 (MQ=255)
   aGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                <  1:2335975/100‑1 (MQ=255)
   aGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                <  2:2402084/100‑1 (MQ=255)
      ggggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAttt                                                                                             >  1:1316500/1‑100 (MQ=255)
       gggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTg                                                                                            <  2:436626/100‑1 (MQ=255)
       gggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTg                                                                                            <  2:2724647/100‑1 (MQ=255)
         gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGcc                                                                                          <  2:705431/100‑1 (MQ=255)
         gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGcc                                                                                          <  1:95414/100‑1 (MQ=255)
          tGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGccn                                                                                         <  1:2071056/100‑2 (MQ=255)
           gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                        >  1:4386459/1‑100 (MQ=255)
           gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                        <  2:2070792/100‑1 (MQ=255)
           gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                        >  1:4386569/1‑100 (MQ=255)
            aCGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtcc                                                                                       <  2:225389/100‑3 (MQ=255)
            aCGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtcc                                                                                       <  2:3773270/100‑3 (MQ=255)
               cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccgca                                                                                    <  2:3391186/100‑6 (MQ=255)
               cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtattat                                                                                    >  1:1021514/1‑100 (MQ=255)
                 agaAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGt                                                                                  <  2:3355355/100‑1 (MQ=255)
                 agaAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGt                                                                                  >  1:4090033/1‑100 (MQ=255)
                 agaAAGCATCTCATTAAGCGCAGTGTCTCGCGACATCTATCCGGTGGTAAATGCTATGACAGCCCTCCTGCGCACGCAAAATCGAAATTGCCta                                                                                        <  1:3662891/94‑1 (MQ=255)
                   aaaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccgcaggag                                                                                <  2:2123449/100‑10 (MQ=38)
                    aaGCATCTCATTAAGCGCGGTGTCGAGCGCCATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  1:3720665/100‑1 (MQ=255)
                    aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  2:2904971/100‑1 (MQ=255)
                    aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               <  2:1429482/100‑1 (MQ=255)
                    aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAg                                                                               >  1:334927/1‑100 (MQ=255)
                      gCATCTCATTAAGCGCGGTGTCGGGCGACATCAATCCGGTGGTAAATTCCATGAAAACTCTCCTGCGGCGGCAAATTCGATTTTGCCtccg                                                                                      <  1:178147/91‑4 (MQ=255)
                       cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                            >  2:4284607/1‑100 (MQ=255)
                       cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                            <  1:2197458/100‑1 (MQ=255)
                          ctcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAAc                                                                         >  1:3428489/1‑100 (MQ=255)
                            cATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                       >  2:3481256/1‑100 (MQ=255)
                              ttAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCca                                                                     <  2:1344992/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                                                              <  2:3750182/100‑1 (MQ=255)
                                     cgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTg                                                              <  1:3481238/100‑1 (MQ=255)
                                     cgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTg                                                              <  2:3827033/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)
                                           cGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcc                                                        >  2:3423536/1‑100 (MQ=255)
                                           cGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcc                                                        >  1:3057604/1‑100 (MQ=255)
                                            gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGTCa                                                       >  1:2515701/1‑100 (MQ=255)
                                            gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                       >  2:766375/1‑100 (MQ=255)
                                            gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                       <  2:3174633/100‑1 (MQ=255)
                                            gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                       >  2:1380997/1‑100 (MQ=255)
                                            gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                       >  2:3309498/1‑100 (MQ=255)
                                               cGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCt                                                    >  1:41970/1‑100 (MQ=255)
                                                gACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTc                                                   <  2:1716106/100‑1 (MQ=255)
                                                gACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTc                                                   >  2:87260/1‑100 (MQ=255)
                                                  cATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCgg                                                 >  2:4175587/1‑100 (MQ=255)
                                                   aTCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCggg                                                <  1:1866430/100‑1 (MQ=255)
                                                          cGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAgcg                                         >  2:367630/1‑100 (MQ=255)
                                                             tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTg                                      <  1:1597554/100‑1 (MQ=255)
                                                             tggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTg                                      <  2:3559884/100‑1 (MQ=255)
                                                                tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGt                                   >  2:1237209/1‑100 (MQ=255)
                                                                tAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGt                                   <  1:3897120/100‑1 (MQ=255)
                                                                  aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                 <  2:963705/100‑1 (MQ=255)
                                                                  aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                 <  1:958158/100‑1 (MQ=255)
                                                                   aTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                <  1:1372181/100‑1 (MQ=255)
                                                                   aTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                >  1:1391116/1‑100 (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:3806247/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)
                                                                        aTGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCtt                           >  2:1508480/1‑100 (MQ=255)
                                                                          gAAAACGCTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                         <  2:935740/100‑1 (MQ=255)
                                                                          gAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                         <  2:1761199/100‑1 (MQ=255)
                                                                          gAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATGCTGTCAGAACAACGCCACACCCTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                         <  1:2396125/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)
                                                                             aaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                       <  1:688020/99‑1 (MQ=255)
                                                                               cacTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGc                    <  2:2650051/100‑1 (MQ=255)
                                                                                acTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGcc                   <  2:2247762/100‑1 (MQ=255)
                                                                                acTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGcc                   <  1:1391511/100‑1 (MQ=255)
                                                                                 cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                  >  2:1617968/1‑100 (MQ=255)
                                                                                 cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                  <  2:2602668/100‑1 (MQ=255)
                                                                                    cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTTGGTCCTGCCTTTACATGCCatat               >  1:339273/7‑100 (MQ=38)
                                                                                    cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  1:2608368/7‑100 (MQ=39)
                                                                                    cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  1:2608293/7‑100 (MQ=39)
                                                                                    cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACATCGGGCGTAGCACTGGGTCCTGCCTTTACATGCCaca                >  1:2413053/7‑97 (MQ=38)
                                                                                    cataataGGCAAACTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat               >  1:2268376/7‑100 (MQ=255)
                                                                                      taataGGCAAATACGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc             >  2:2893197/5‑100 (MQ=255)
                                                                                       aataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCa            >  1:714040/4‑100 (MQ=255)
                                                                                         gAGGCAAAATCGAATTTGCCTATTATGTCAGAAAACCGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCa            <  1:1787692/98‑1 (MQ=255)
                                                                                         gAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATc          >  2:1640454/1‑100 (MQ=255)
                                                                                          aGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCt         >  2:607357/1‑100 (MQ=255)
                                                                                             cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCtt        >  1:4029055/1‑98 (MQ=255)
                                                                                             cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc      >  2:3100092/1‑100 (MQ=255)
                                                                                              aaaaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCt     >  2:40104/1‑100 (MQ=255)
                                                                                               aaaTCGAACTTGCCTATTATGTCAGAAAAACGCCACACACTGTATGCCAACTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCta    <  2:2392707/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATA  >  NC_000913/866291‑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: