Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,770,087 C→A 100% pseudogene (55/114 nt) yibU → pseudogene, H repeat‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,0870CA97.6% 402.8 / 0.6 127pseudogene (55/114 nt)yibUpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (59/64);  total (61/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 129 total aligned reads displayed.

GGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAA  >  NC_000913/3769990‑3770180
                                                                                                 |                                                                                             
ggAAAAGCACTCCTGCTGTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAg                                                                                             <  1:565094/100‑1 (MQ=255)
 gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                            >  1:372019/1‑100 (MQ=255)
 gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                            >  2:2783791/1‑100 (MQ=255)
 gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                            >  1:1873700/1‑100 (MQ=255)
  aaaaGCACTCCTGCTCTTTTTTGATTAGATCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGgaagaa                                                                                           <  2:1819651/100‑1 (MQ=255)
  aaaaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGATGATGTATAGTCTCATCCTGCGgaagaa                                                                                           <  2:563828/100‑1 (MQ=255)
  aaaaGCACTCCTCCTCTTTTTTGATTAGAGCCGTCGCAGGAGTAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                            <  2:2778788/99‑1 (MQ=255)
   aaaGCACTCCTGCTCTTTTTTTTTTTTTTTCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAAc                                                                                          >  2:2006698/1‑100 (MQ=255)
   aaaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAAc                                                                                          <  1:1089105/100‑1 (MQ=255)
     aGCACTCCTGCTCTTTTTTTTTTTTTTTTGTCGCAGGAGGAACAACTCATATAATCATTGCGATTTCGAAGATGTATAGTCTCATTCTGCGGAAGAAc                                                                                          >  2:2905634/1‑98 (MQ=255)
      gCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAg                                                                                       <  2:929056/100‑1 (MQ=255)
          tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGTAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCTTCCTGCGGAAGAACAAGACgg                                                                                   <  2:2575275/100‑1 (MQ=255)
           ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGAGATCTCGACGATGTATATTCTCATCCTGCAGAAGAACAAGACGGa                                                                                  >  2:2771658/1‑100 (MQ=255)
             tGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTAATCCTGCGGAAGAACAAGACGGATa                                                                                >  1:1889908/1‑100 (MQ=255)
              gCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaa                                                                               >  1:1356092/1‑100 (MQ=255)
                tctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaaaa                                                                             >  1:2245190/1‑100 (MQ=255)
                    ttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATc                                                                         >  1:269734/1‑100 (MQ=255)
                     tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                        <  2:2690713/100‑1 (MQ=255)
                     tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                        <  2:1308332/100‑1 (MQ=255)
                     tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                        <  2:1308999/100‑1 (MQ=255)
                     tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                        >  1:1613951/1‑100 (MQ=255)
                      ttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                        <  1:834898/99‑1 (MQ=255)
                       tGATTAGAGTCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCGCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                      <  1:961420/100‑1 (MQ=255)
                       tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                      >  2:1340289/1‑100 (MQ=255)
                       tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                      <  2:613872/100‑1 (MQ=255)
                         aTTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTTCGTTCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATg                                                                    <  1:1705894/100‑1 (MQ=255)
                           tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                  <  2:2525391/100‑1 (MQ=255)
                           tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                  <  2:2525346/100‑1 (MQ=255)
                            agagCCGTCGCTGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                 <  2:1779561/100‑1 (MQ=255)
                            agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                 >  1:812124/1‑100 (MQ=255)
                             gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaggg                                                                <  2:2086181/100‑3 (MQ=255)
                             gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                <  2:2835223/100‑1 (MQ=255)
                             gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                >  1:2322312/1‑100 (MQ=255)
                             gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                <  2:2014363/100‑1 (MQ=255)
                             aaaaaagagCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                <  2:334296/92‑1 (MQ=255)
                                 cGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGcaca                                                            <  2:447630/100‑1 (MQ=255)
                                  gtagCAGGAGGAACAAATAATATCATCATCGCGATCTCGACGATGTATAGTCTCATCCCGCGGAAGAACAACACGGCTAAAAAATCCACTGCGATCacac                                                           >  2:758741/4‑100 (MQ=255)
                                  gTCGCAGGAGGAACAACTCATATCATTATTTCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCacac                                                           <  1:2714283/100‑1 (MQ=255)
                                   tCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACt                                                          <  1:2006721/100‑1 (MQ=255)
                                     gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                        >  1:1328817/1‑100 (MQ=255)
                                     gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                        <  1:960651/100‑1 (MQ=255)
                                        ggaggaGCAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                     <  2:305371/100‑1 (MQ=255)
                                        ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                     <  2:2925434/100‑1 (MQ=255)
                                        ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                     <  2:2925485/100‑1 (MQ=255)
                                        ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAAAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                     >  1:158171/1‑100 (MQ=255)
                                         gaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCa                                                    <  1:320532/100‑1 (MQ=255)
                                          aggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAg                                                   <  2:2930022/100‑1 (MQ=255)
                                             aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGTACAACAAGACGGATAAAAAATCCAATGAAAGCACACTCATCCAGaa                                                 >  1:2099073/1‑99 (MQ=255)
                                             aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                <  2:1578054/100‑1 (MQ=255)
                                             aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCGCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                <  1:2088904/100‑1 (MQ=255)
                                              acaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACt                                               <  1:2621695/100‑1 (MQ=255)
                                               caacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtt                                              >  2:3029175/1‑100 (MQ=255)
                                                 acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttat                                            >  1:2026987/1‑100 (MQ=255)
                                                 acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttat                                            <  2:1438831/100‑1 (MQ=255)
                                                 acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttat                                            <  2:225449/100‑1 (MQ=255)
                                                  cTCATATCATTATTGTTTTTTCGACGTTGTATAGTCTCTTCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                           <  1:293519/100‑1 (MQ=255)
                                                   tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                          >  2:866799/1‑100 (MQ=255)
                                                    cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                         >  1:179547/1‑100 (MQ=255)
                                                    cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                         >  2:1484631/1‑100 (MQ=255)
                                                     atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                        <  2:1067993/100‑1 (MQ=255)
                                                     atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                        <  1:628865/100‑1 (MQ=255)
                                                     atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                        <  1:195750/100‑1 (MQ=255)
                                                     atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                        >  1:222558/1‑100 (MQ=255)
                                                      tatCATCATTGCGATCTCGACGATGTATAGTCTCATGCAGCGGAGGAACAAGACGGATAATAAATCCAATGAGAGCACACTCATCTAGAACTTATTAACa                                       >  2:1595712/1‑100 (MQ=255)
                                                      tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                       <  2:2792392/100‑1 (MQ=255)
                                                       atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttaccaac                                        >  2:2693720/1‑93 (MQ=255)
                                                       atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                      >  1:1746236/1‑100 (MQ=255)
                                                       atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                      >  1:429427/1‑100 (MQ=255)
                                                         catcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATaa                                    >  1:1270866/1‑100 (MQ=255)
                                                            catTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGa                                 <  1:2727226/100‑1 (MQ=255)
                                                             atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                >  2:2971693/1‑100 (MQ=255)
                                                             atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGa                                 <  1:2682036/99‑1 (MQ=255)
                                                              ttGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                >  2:665952/1‑99 (MQ=255)
                                                               tGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAata                              <  2:888221/100‑1 (MQ=255)
                                                                gCGATCTCGACGATGTATAGTCTCATCGTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                             <  2:631542/100‑1 (MQ=255)
                                                                gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                             >  2:214201/1‑100 (MQ=255)
                                                                 cGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                             <  1:518529/99‑1 (MQ=255)
                                                                 cGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATAAACATAATGAatat                             >  1:123872/1‑99 (MQ=255)
                                                                  gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                           >  1:828435/1‑100 (MQ=255)
                                                                  gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                           <  1:214203/100‑1 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGTACAAGTCGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                        <  1:2599598/100‑1 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCAGACTCATCCAGAACTTATTAACATAATGAATATTAAAg                        <  2:1613932/100‑1 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                        <  2:214457/100‑1 (MQ=255)
                                                                      tcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGa                       <  1:1524877/100‑1 (MQ=255)
                                                                      tcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGa                       >  2:108960/1‑100 (MQ=255)
                                                                      tcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAAAGAATATaaaan                        <  2:515900/99‑2 (MQ=255)
                                                                       cgacgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaa                      >  1:1580312/1‑100 (MQ=255)
                                                                       cgacgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaa                      >  2:1438770/1‑100 (MQ=255)
                                                                         acgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaa                    <  2:2350635/100‑1 (MQ=255)
                                                                            aTGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatc                 <  2:179545/100‑1 (MQ=255)
                                                                             tGTATAGTCTCATCCTGCGGAATAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                <  1:234971/100‑1 (MQ=255)
                                                                              gTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcat               >  1:2286823/1‑100 (MQ=255)
                                                                               tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc              >  1:2722880/1‑100 (MQ=255)
                                                                               tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc              >  1:2252202/1‑100 (MQ=255)
                                                                                 taGTCTCATCCTGCGGAAAAACAAAACTGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGa            <  1:1142654/100‑1 (MQ=255)
                                                                                   gTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg          >  2:1911895/1‑100 (MQ=255)
                                                                                         tCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATg    >  2:1512045/1‑100 (MQ=255)
                                                                                          ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTGTTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa   >  2:2129813/1‑100 (MQ=255)
                                                                                          ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa   >  2:2819424/1‑100 (MQ=255)
                                                                                           cTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa  >  1:1132504/1‑100 (MQ=255)
                                                                                                 |                                                                                             
GGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAA  >  NC_000913/3769990‑3770180

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