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,0870CA96.7% 323.2 / 0.9 100pseudogene (55/114 nt)yibUpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (50/46);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTA  >  NC_000913/3769990‑3770183
                                                                                                 |                                                                                                
ggAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAg                                                                                                <  1:1005488/100‑1 (MQ=255)
ggAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAg                                                                                                <  1:1328842/100‑1 (MQ=255)
ggAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAg                                                                                                <  2:386894/100‑1 (MQ=255)
  aaaaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGgaagaa                                                                                              <  1:684350/100‑1 (MQ=255)
   aaaGCACTCCTGCTCTTTTTTGAGTAGAGCCGTCGCAGTAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAAc                                                                                             <  1:1492863/100‑1 (MQ=255)
    aaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACa                                                                                            <  1:2365736/100‑1 (MQ=255)
     aGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACaa                                                                                           <  1:2550239/100‑1 (MQ=255)
     aGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACaa                                                                                           >  1:1623031/1‑100 (MQ=255)
       cACTCCTGCTCTTTTTTAATTAGAGCCGTCGCAGACGGAACCACTAATATCATCATTGCGATCTCGACGATGTATAATCTCATCCTGCGGAAGAACAAGa                                                                                         >  2:491200/1‑100 (MQ=255)
         cTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGGCTCATCCTGCGGAAGAACAAGACg                                                                                       >  1:1151239/1‑100 (MQ=255)
           ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGa                                                                                     <  1:1128028/100‑1 (MQ=255)
           ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGa                                                                                     >  1:730205/1‑100 (MQ=255)
            cTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGAt                                                                                    <  2:1693098/100‑1 (MQ=255)
            cTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGAt                                                                                    <  2:1914889/100‑1 (MQ=255)
             tGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATa                                                                                   <  1:2077886/100‑1 (MQ=255)
               ctctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaaa                                                                                 >  1:692390/1‑100 (MQ=255)
               ctctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaaa                                                                                 >  1:2331666/1‑100 (MQ=255)
                 ctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTAATCCTGCGAAAGAACAAGAAGGATaaaaa                                                                               >  2:1238959/1‑100 (MQ=255)
                 ctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAATCTCATCCTGCGGAAGAACAAGACGGATaaaaa                                                                               >  2:633005/1‑100 (MQ=255)
                 ctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGAAGCTGTATAGTCTCATACTGCGTAAGAACAAGACGGATaaaaa                                                                               >  2:1386833/1‑100 (MQ=255)
                   tttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAAt                                                                             <  2:2487426/100‑1 (MQ=255)
                   tttttGATTAGAGCAGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGAAGATGTATGTTCTCATCCTGCGGAAGAACAAGACGGATAAAAAAt                                                                             >  2:606149/1‑100 (MQ=255)
                     tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                           >  1:1381207/1‑100 (MQ=255)
                      ttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCAACCTGCGGAAGAACAAGACGGATAAAAAATCCa                                                                          <  1:1668109/100‑1 (MQ=255)
                       tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                         <  1:1872820/100‑1 (MQ=255)
                        gATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAAt                                                                        >  2:1595320/1‑100 (MQ=255)
                        gATTAGAGCCGTCGCAGGAGGAACAACACATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                         <  2:518410/99‑1 (MQ=255)
                         aTTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATg                                                                       <  2:1975113/100‑1 (MQ=255)
                            agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                    >  1:123623/1‑100 (MQ=255)
                             gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                    >  2:861596/1‑99 (MQ=255)
                             gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                   <  1:1489680/100‑1 (MQ=255)
                              agCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGc                                                                  <  1:758392/100‑1 (MQ=255)
                              agCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGc                                                                  <  1:1521147/100‑1 (MQ=255)
                              agCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGc                                                                  <  1:393138/100‑1 (MQ=255)
                               gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGca                                                                 >  1:2370126/1‑100 (MQ=255)
                               gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGca                                                                 >  2:2138410/1‑100 (MQ=255)
                                ccGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGcac                                                                >  2:58367/1‑100 (MQ=255)
                                  gTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGca                                                                 >  2:2332713/1‑97 (MQ=255)
                                  gTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCacac                                                              >  1:1654746/1‑100 (MQ=255)
                                  gTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCacac                                                              <  2:237915/100‑1 (MQ=255)
                                   tCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACt                                                             <  2:1437248/100‑1 (MQ=255)
                                        ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGTATAAAAAATCCAATGAGAGCACACTCATcc                                                        >  2:920943/1‑100 (MQ=255)
                                        ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                        >  2:1165727/1‑100 (MQ=255)
                                           ggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                     <  1:312110/100‑1 (MQ=255)
                                              acaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACt                                                  <  2:1399571/100‑1 (MQ=255)
                                                aacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                <  2:2155617/100‑1 (MQ=255)
                                                aacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                >  2:653602/1‑100 (MQ=255)
                                                 acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                <  2:239622/99‑1 (MQ=255)
                                                  cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                              <  1:1994696/100‑1 (MQ=255)
                                                  cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                              <  1:2208368/100‑1 (MQ=255)
                                                   tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                             <  1:824811/100‑1 (MQ=255)
                                                    cATATCATCATTGCGATCTCGACGATGTCTAGTCTCCTCCTGCGGAAGAACAAGACGGATAAAACATCCGACGAGAGGAAACTCATCCCGAACttatta                                             >  2:748982/1‑99 (MQ=255)
                                                    cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCCATGAGAGCACACTCATCCAGAACCTATTaa                                            >  2:245843/1‑100 (MQ=255)
                                                    cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                            <  1:37407/100‑1 (MQ=255)
                                                     atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                           >  1:259851/1‑100 (MQ=255)
                                                      tataATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCTGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                           <  2:1409175/99‑1 (MQ=255)
                                                      tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                          <  1:1620683/100‑1 (MQ=255)
                                                      tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                          >  2:1620671/1‑100 (MQ=255)
                                                       atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                         <  2:628581/100‑1 (MQ=255)
                                                       atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                         >  2:1380928/1‑100 (MQ=255)
                                                       atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                         >  2:1288840/1‑100 (MQ=255)
                                                        tcatcaTTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATa                                        >  2:1901288/1‑100 (MQ=255)
                                                        tcatcaTTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                         <  2:2270186/99‑1 (MQ=255)
                                                         catcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATaa                                       >  2:1616258/1‑100 (MQ=255)
                                                          atcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAAt                                      <  1:1901302/100‑1 (MQ=255)
                                                          atcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAAt                                      <  2:2248714/100‑1 (MQ=255)
                                                            catTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACACGACGGGTAAAAAATCCACTGAGAGCACACACATCCAGAACTTATTAACATAATGa                                    >  2:2355492/1‑100 (MQ=255)
                                                             atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                   >  2:815333/1‑100 (MQ=255)
                                                              ttGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAat                                  >  2:504189/1‑100 (MQ=255)
                                                              ttGCGAACCCGACGATGTAGAGTCTCATCCTGCGGAAGAACAAGACTGATAAAAGAGCCAAGGAGAGCCCACTCGGCCAGAACCTAGTAACATAATGagt                                  <  2:448542/100‑3 (MQ=255)
                                                                gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                                <  1:1196080/100‑1 (MQ=255)
                                                                  gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                              >  2:1969181/1‑100 (MQ=255)
                                                                  gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                              >  2:1671061/1‑100 (MQ=255)
                                                                  gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                              >  2:386046/1‑100 (MQ=255)
                                                                   aTCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaa                             <  1:653609/100‑1 (MQ=255)
                                                                   aTCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaa                             >  2:276899/1‑100 (MQ=255)
                                                                    tctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaaa                            <  2:658448/100‑1 (MQ=255)
                                                                    tctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaaa                            <  1:2473662/100‑1 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaaa                            <  2:2059439/99‑1 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                           >  1:1323559/1‑100 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                           >  2:1800281/1‑100 (MQ=255)
                                                                     ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                           >  1:939927/1‑100 (MQ=255)
                                                                      tcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCCAAACTTATTAACATAATGACTATTAAAGa                          >  2:355554/1‑100 (MQ=255)
                                                                      tcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGa                          >  1:1181927/1‑100 (MQ=255)
                                                                          cgaTGTATAGTCTCCTCCTGCGGAACAACAAGACGGATCAAAAATCCAATGAGAGCACACTCCTACAGAACTTACTAACATAATGACTATTACAGcaaca                      >  2:2356268/1‑98 (MQ=255)
                                                                          cgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaaa                      >  2:447348/1‑100 (MQ=255)
                                                                               tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc                 >  2:756754/1‑100 (MQ=255)
                                                                               tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc                 >  1:1061855/1‑100 (MQ=255)
                                                                               tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTAGTAACATAATGAATATTAAAGAAAAatcatc                 >  2:1894150/1‑100 (MQ=255)
                                                                                ataatgtcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCg                <  2:303161/94‑1 (MQ=255)
                                                                                ataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCg                <  1:1873771/100‑1 (MQ=255)
                                                                                ataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCg                <  1:1605512/100‑1 (MQ=255)
                                                                                  aGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGaa              <  2:1061847/100‑1 (MQ=255)
                                                                                   gTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg             >  2:1315783/1‑100 (MQ=255)
                                                                                     ctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCt           >  1:2089744/1‑100 (MQ=255)
                                                                                        aTCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAAt        >  1:2113616/1‑100 (MQ=255)
                                                                                        aTCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAAt        <  1:756761/100‑1 (MQ=255)
                                                                                           cTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa     >  2:2226947/1‑100 (MQ=255)
                                                                                            tGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAg    <  1:1892654/100‑1 (MQ=255)
                                                                                              cGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTa  <  1:2415538/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTA  >  NC_000913/3769990‑3770183

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: