Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,393 A→T 15.1% 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,3930AT15.1% 336.8 / 44.6 131intergenic (‑29/‑175)moeA/iaaAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (47/64);  new base T (10/10);  total (57/74)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 133 total aligned reads displayed.

GTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAA  >  NC_000913/866295‑866492
                                                                                                  |                                                                                                   
gTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAt                                                                                                    <  1:2485849/100‑1 (MQ=255)
  ggggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAttt                                                                                                  <  1:388967/100‑1 (MQ=255)
    ggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                >  2:1408794/1‑100 (MQ=255)
     gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGcc                                                                                               >  2:2710993/1‑100 (MQ=255)
     gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGcc                                                                                               <  2:1132239/100‑1 (MQ=255)
     gTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGcc                                                                                               <  1:1545672/100‑1 (MQ=255)
       gACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                             <  2:1400331/100‑1 (MQ=255)
       gACGCGAGAAAGCATCTCATTAAGAGCGGTGTCGAGCGACATCAATCCGGTGGTAAAAACAAAAAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtc                                                                                             <  2:2954754/100‑2 (MQ=255)
         cgcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccg                                                                                           <  1:3340520/100‑4 (MQ=255)
         cgcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccg                                                                                           <  2:278336/100‑4 (MQ=255)
         cgcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccg                                                                                           <  1:855432/100‑4 (MQ=255)
          gcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccgc                                                                                          <  2:3196329/100‑5 (MQ=255)
           cgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtccgca                                                                                         <  1:2131159/100‑6 (MQ=255)
               aaaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                     <  2:1078668/100‑1 (MQ=255)
               aaaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCa                                                                                     <  1:2876619/100‑1 (MQ=255)
                 aGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGa                                                                                   <  1:2775290/100‑1 (MQ=255)
                  gCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaa                                                                                  >  2:456332/1‑100 (MQ=255)
                   cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaa                                                                                 >  1:2630989/1‑100 (MQ=255)
                    aTCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaaa                                                                                >  1:122910/1‑100 (MQ=255)
                        cATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAAc                                                                              >  1:1442127/1‑98 (MQ=255)
                        cATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                            <  2:2870872/100‑1 (MQ=255)
                        cATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                            >  1:2776901/1‑100 (MQ=255)
                         aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGcc                                                                           <  1:2673557/100‑1 (MQ=255)
                         aTTAAGCGCAGTGTCGACCGACATCAATCCGGTCGTAAATTCAATGAAATCTCTCCTGCAGAGGCAACATCGCATTTGCCTATTATGTCAGAAAAACGcc                                                                           <  2:525419/100‑1 (MQ=255)
                            aaGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCcaca                                                                        <  2:1637989/100‑1 (MQ=255)
                            aaGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCcaca                                                                        >  1:3322509/1‑100 (MQ=255)
                              gcgcGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGa                                                                      >  1:2539567/1‑100 (MQ=255)
                                 nggTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGAc                                                                     >  2:2205854/2‑98 (MQ=255)
                                  ggTGTCGAGCGACATCAATGCGGTGGTAAATTCCATGAAAACACTCCGGCGTAGGCAAAATCGAATTGGCCTATTCTGTCAGAAAAACGCCACAGACTGt                                                                  <  2:351605/100‑1 (MQ=255)
                                  ggTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGt                                                                  <  2:2055346/100‑1 (MQ=255)
                                    tgtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTAt                                                                >  1:1502062/1‑100 (MQ=255)
                                        gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCa                                                            <  1:456334/100‑1 (MQ=255)
                                         aGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcc                                                             <  1:2520923/98‑1 (MQ=255)
                                          gCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCAcc                                                          >  2:2625081/1‑100 (MQ=255)
                                          gCGACATCAATCCGGGGGTAAATTCTATGAAAACACGCCGGCGGAGGCAAAATCGAATTTGCCTATTAGGTCAGAAAAACGCCACAGACTGTATGCCAcc                                                          <  2:710559/100‑1 (MQ=255)
                                            gACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTc                                                        >  1:1875125/1‑100 (MQ=255)
                                            gACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTc                                                        >  2:2111367/1‑100 (MQ=255)
                                             aCATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCg                                                       <  2:2637343/100‑1 (MQ=255)
                                               aTCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCggg                                                     <  2:495531/100‑1 (MQ=255)
                                                 cAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCg                                                   >  2:2505130/1‑100 (MQ=255)
                                                  aaTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGt                                                  >  1:2910210/1‑100 (MQ=255)
                                                   aTCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTa                                                 >  1:849244/1‑100 (MQ=255)
                                                    tCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAg                                                >  2:616413/1‑100 (MQ=255)
                                                     ccGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAgc                                               >  2:2806305/1‑100 (MQ=255)
                                                      cGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAgcg                                              <  2:3072807/100‑1 (MQ=255)
                                                      cGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAgcg                                              >  1:3072830/1‑100 (MQ=255)
                                                       ggtggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAgcgc                                             <  2:1897121/100‑1 (MQ=255)
                                                          ggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTgg                                          <  1:1669569/100‑1 (MQ=255)
                                                          ggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTgg                                          >  2:1660649/1‑100 (MQ=255)
                                                           gtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTggg                                         <  1:2111380/100‑1 (MQ=255)
                                                           gtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTggg                                         <  2:2951327/100‑1 (MQ=255)
                                                             acaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTc                                       >  1:2467256/3‑100 (MQ=255)
                                                             aaaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTc                                       <  1:1664386/100‑1 (MQ=255)
                                                             aaaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTc                                       >  2:22293/1‑100 (MQ=255)
                                                              aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                      >  2:2452317/1‑100 (MQ=255)
                                                              aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                      >  2:2803192/1‑100 (MQ=255)
                                                              aaTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTcc                                      >  1:1388373/1‑100 (MQ=255)
                                                               aTTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                     >  2:1773901/1‑100 (MQ=255)
                                                                   cATGAAAACACTCCTGCGGAGGCAAAATAGAATTTGCCTATTATGTCAGAAAAACGCCATAGACTGAATGTCAACTAGAGTGTAGCGATGGGTcc                                      >  2:3150738/1‑95 (MQ=255)
                                                                    aTGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCtt                                >  2:1167418/1‑100 (MQ=255)
                                                                    aTGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCtt                                <  2:530985/100‑1 (MQ=255)
                                                                    aTGAAAACACTCCTGCAGAGGCAAAATCGAATTTGCATATTAAATCAAAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGcc                                  <  2:3052369/98‑1 (MQ=255)
                                                                     tGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCATATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCttt                               >  2:1810022/1‑100 (MQ=255)
                                                                      gAAAACACTCCTGCGGAGGCAAAATCGATGTTGCCTATTAGGTCAGAAAAACGCCACAGACTGTTTGCCACCTCGGGCGTATCGCTGGTTCCTGcc                                  <  2:2511543/96‑1 (MQ=255)
                                                                      gAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTa                              >  2:2336067/1‑100 (MQ=255)
                                                                       aaaaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCCCAGACTGTATGCCAACTCGCGCGTAGCGCTGGGTCCTGCCTTTAc                             >  1:1253613/1‑100 (MQ=255)
                                                                       aaaaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGGCAGATAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTAc                             <  2:698984/100‑1 (MQ=255)
                                                                        aaaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACa                            >  2:1230531/1‑100 (MQ=255)
                                                                            acTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGcc                        <  1:1030198/100‑1 (MQ=255)
                                                                             tgacataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                       >  1:853394/10‑100 (MQ=38)
                                                                             cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                       <  2:1406128/100‑1 (MQ=255)
                                                                             cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                       <  1:527619/100‑1 (MQ=255)
                                                                             cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                       >  2:527617/1‑100 (MQ=255)
                                                                              tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCat                      <  1:3031161/100‑1 (MQ=255)
                                                                              gacataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCat                      >  1:3176305/9‑100 (MQ=38)
                                                                              gacataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCa                       >  1:2365687/9‑99 (MQ=38)
                                                                               ccTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCata                     <  1:1109894/100‑1 (MQ=255)
                                                                               ccTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCata                     >  1:1296793/1‑100 (MQ=255)
                                                                                cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                    >  2:795524/7‑100 (MQ=39)
                                                                                cataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                    >  1:1363025/7‑100 (MQ=39)
                                                                                cTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                    >  1:921156/1‑100 (MQ=255)
                                                                                cTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                    >  1:1745685/1‑100 (MQ=255)
                                                                                cTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatat                    >  1:3496596/1‑100 (MQ=255)
                                                                                 tGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATc                   <  1:134163/100‑1 (MQ=255)
                                                                                 tGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATc                   <  1:2379560/100‑1 (MQ=255)
                                                                                 ataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATc                   >  1:116701/6‑100 (MQ=255)
                                                                                  gCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc                  >  2:1201167/1‑100 (MQ=255)
                                                                                  gCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATcc                  <  1:439712/100‑1 (MQ=255)
                                                                                      aGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGAGCTGGGTCATGCCTTTACATGCCATATCCATc               >  2:1601034/1‑99 (MQ=255)
                                                                                       ggCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCtt             <  1:331788/100‑1 (MQ=255)
                                                                                        gCAAAATCGANTTTGCNTCTTATGTTAGAANAANGCCACAGNATGTNTGCCCTCTCGGGCGTAGCGCTGNGTCCTGCCTTTACATGCCCTCTACATCt              >  2:3220814/1‑98 (MQ=255)
                                                                                         cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc           <  1:3247502/100‑1 (MQ=255)
                                                                                         cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc           >  2:2028166/1‑100 (MQ=255)
                                                                                         cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc           <  1:37903/100‑1 (MQ=255)
                                                                                         cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc           >  1:1006483/1‑100 (MQ=255)
                                                                                         cAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCACTGGGTCCTGCCTTTACATGCCATATCCATCTTTc           <  1:2260569/100‑1 (MQ=255)
                                                                                           aaaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCta         <  1:906498/100‑1 (MQ=255)
                                                                                            aaTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGTTCCTGCCTTTACATGCCATATCCATCTTTCta         <  1:274315/99‑1 (MQ=255)
                                                                                              tCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTatat      >  1:883883/1‑100 (MQ=255)
                                                                                               cGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATAtt     <  2:1083953/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAA  >  NC_000913/866295‑866492

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: