Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,770,087 C→A 24.6% 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,0870CA24.6% 353.6 / 120.3 183pseudogene (55/114 nt)yibUpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (65/73);  new base A (20/25);  total (85/98)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 191 total aligned reads displayed.

ATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTAATA  >  NC_000913/3769988‑3770186
                                                                                                   |                                                                                                   
aTGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGc                                                                                                     <  1:1675877/100‑1 (MQ=255)
 tGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCa                                                                                                    >  1:3298525/1‑100 (MQ=255)
  ggAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAg                                                                                                   <  2:644453/100‑1 (MQ=255)
   gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGa                                                                                                  <  1:1132791/100‑1 (MQ=255)
   gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGa                                                                                                  >  1:2582935/1‑100 (MQ=255)
   gAAAAACACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGa                                                                                                  >  1:3830496/1‑100 (MQ=255)
    aaaaGCACTCCTGCTCTTTTTTTTTTTTTTTTTTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGCTGTATAGTCTCATCCTGCGgaagaa                                                                                                 >  2:1033484/1‑100 (MQ=255)
     aaaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAAc                                                                                                >  1:1241641/1‑100 (MQ=255)
        gCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAg                                                                                             >  1:1202096/1‑100 (MQ=255)
        gCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAg                                                                                             >  1:2933754/1‑100 (MQ=255)
        gCACTCCTGCTCTTTTTTGATTAGAGATGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTGTCATCCTGCGGCAGAACAAg                                                                                             >  2:3561137/1‑100 (MQ=255)
         cACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGa                                                                                            <  2:1137748/100‑1 (MQ=255)
         cACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGa                                                                                            >  2:2977870/1‑100 (MQ=255)
           cTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACg                                                                                          <  2:2110922/100‑1 (MQ=255)
            tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACgg                                                                                         <  1:1945317/100‑1 (MQ=255)
             ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGCTCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGa                                                                                        >  1:3933034/1‑100 (MQ=255)
              cTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACACCTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACCAGACgg                                                                                         >  2:2592507/1‑98 (MQ=255)
                   ctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATaaaaa                                                                                  >  1:1821990/1‑100 (MQ=255)
                    ttttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATaaaca                                                                                  >  2:72567/1‑97 (MQ=255)
                     tttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGTGGCAGAACAAGACGGATAAAAAAt                                                                                <  2:3040200/100‑1 (MQ=255)
                       tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATcc                                                                              <  1:2304635/100‑1 (MQ=255)
                        ttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCa                                                                             <  2:2340649/100‑1 (MQ=255)
                         tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCaa                                                                            <  1:1919321/100‑1 (MQ=255)
                         tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCaa                                                                            >  1:3720774/1‑100 (MQ=255)
                         tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCaa                                                                            <  2:1239245/100‑1 (MQ=255)
                         tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGACTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCaa                                                                            >  1:4291/1‑100 (MQ=255)
                          gATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAAt                                                                           >  1:3407543/1‑100 (MQ=255)
                          gATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                            >  1:214950/1‑99 (MQ=255)
                           aTTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCAATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATg                                                                          <  2:1325771/100‑1 (MQ=255)
                            ttAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATga                                                                         >  2:2427866/1‑100 (MQ=255)
                             tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                        <  2:150813/100‑1 (MQ=255)
                              ngagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATGCTGCGACAGAACAAGACGGATAAAAAATCCAATgaga                                                                       >  2:2786260/2‑100 (MQ=255)
                               gagaCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGagag                                                                      <  2:1475320/100‑1 (MQ=255)
                               gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGagag                                                                      <  1:3549849/100‑1 (MQ=255)
                                agCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGagagn                                                                     <  2:3637995/100‑2 (MQ=255)
                                 gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGca                                                                    <  1:973700/100‑1 (MQ=255)
                                   cGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCAACCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGcac                                                                   <  2:423441/99‑1 (MQ=255)
                                    gTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGGGAGCacac                                                                 <  2:1721236/100‑1 (MQ=255)
                                    gTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCacac                                                                 <  1:242672/100‑1 (MQ=255)
                                     tCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACt                                                                <  1:245146/100‑1 (MQ=255)
                                     tCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACt                                                                <  1:3503056/100‑1 (MQ=255)
                                       gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                              >  1:1450927/1‑100 (MQ=255)
                                       gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                              >  2:2496065/1‑100 (MQ=255)
                                        cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAt                                                             >  1:2372313/1‑100 (MQ=255)
                                         aggaggAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATc                                                            <  1:2294804/100‑1 (MQ=255)
                                          ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                           <  1:1887548/100‑1 (MQ=255)
                                           gaggaACAACTCATATCAACATCTCGAACTCGACGAAGTATAGTCTCATCCTGCGGCAGAAAAAGAAAAATAAAAAATCCAATGAGAGca                                                                    <  2:3698590/90‑1 (MQ=255)
                                            aggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAg                                                         >  2:3482226/1‑100 (MQ=255)
                                               aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTTATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                      <  2:2900831/100‑1 (MQ=255)
                                               aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                      <  2:2943632/100‑1 (MQ=255)
                                               aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                      >  2:1807561/1‑100 (MQ=255)
                                                 caacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtt                                                    <  1:1511199/100‑1 (MQ=255)
                                                  aacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                   <  1:3299145/100‑1 (MQ=255)
                                                   acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                   >  2:3802737/1‑99 (MQ=255)
                                                   acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                   >  2:2612919/1‑99 (MQ=255)
                                                    cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                                 >  1:3068414/1‑100 (MQ=255)
                                                    cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                                 <  2:1969765/100‑1 (MQ=255)
                                                     tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                                <  2:117309/100‑1 (MQ=255)
                                                      cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                               >  1:2102008/1‑100 (MQ=255)
                                                      cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                               <  1:840573/100‑1 (MQ=255)
                                                       atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                              <  2:924384/100‑1 (MQ=255)
                                                       atatCATCAATCCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGAGAAAAAATCCAATGTGAGCACCCTCATCCAGAACTTATTaa                                               <  1:3807985/99‑1 (MQ=255)
                                                        tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                             >  1:362755/1‑100 (MQ=255)
                                                        tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                             <  1:3989459/100‑1 (MQ=255)
                                                         atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                            >  1:3860744/1‑100 (MQ=255)
                                                            atcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAAt                                         <  1:1555029/100‑1 (MQ=255)
                                                            atcagTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAGCCAGAACTTATTAACATAAt                                         <  2:2747733/100‑1 (MQ=255)
                                                             tcatTTCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAAt                                         <  2:588205/99‑1 (MQ=255)
                                                               atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                      <  2:3123620/100‑1 (MQ=255)
                                                               atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                      >  1:1354736/1‑100 (MQ=255)
                                                                 tGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAata                                    >  2:624085/1‑100 (MQ=255)
                                                                  gCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                                   <  1:1807574/100‑1 (MQ=255)
                                                                   ctaTCTCGACGATGTATAGTCTCATCCTGCGACAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATAtt                                  >  1:2782402/3‑100 (MQ=255)
                                                                   cGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATAtt                                  <  1:2704389/100‑1 (MQ=255)
                                                                   cGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATAtt                                  <  2:458516/100‑1 (MQ=255)
                                                                     aTCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaa                                >  2:3770286/1‑100 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                              >  1:1426410/1‑100 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                              <  2:1115982/100‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                              <  2:3222038/100‑1 (MQ=255)
                                                                        tagacgaTGTATAGTATCATCCTGCAGCAGAACAAGAAAGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                                 <  2:1565581/94‑1 (MQ=255)
                                                                         cgacgaTCTATAGCCTCATCCTGCGGCAGAACAAGACGGATACAAAATCCAATGAGAGCACACTCATCGAAAACTAATTAACAAAATGAATAGAAAAGaa                            >  2:2183018/1‑100 (MQ=255)
                                                                          gacgaTGTATAGTCTCAACCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                              <  2:1000960/97‑1 (MQ=255)
                                                                             gaTGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaaa                         <  2:3860719/99‑1 (MQ=255)
                                                                              atatataGACTCATCCTGCGGCAGAAAAAGACGGATAAAAAATCCAATGAGAGCACACACATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatc                       <  2:2831129/97‑1 (MQ=255)
                                                                                  ataGTCTCATCCTGCGGCAGAAAAAAAAAAATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcat                     <  2:2196271/98‑1 (MQ=255)
                                                                                  ataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCg                   <  1:2200550/100‑1 (MQ=255)
                                                                                     gTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg                >  1:2417497/1‑100 (MQ=255)
                                                                                     gTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg                >  2:1430065/1‑100 (MQ=255)
                                                                                       ctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCTATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCt              <  1:2908075/100‑1 (MQ=255)
                                                                                         cATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTaa            >  1:2824543/1‑100 (MQ=255)
                                                                                         cATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTaa            >  1:3567411/1‑100 (MQ=255)
                                                                                           tCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATg          >  2:2126766/1‑100 (MQ=255)
                                                                                             cTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa        >  1:3146822/1‑100 (MQ=255)
                                                                                             cTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa        >  2:4068611/1‑100 (MQ=255)
                                                                                              tGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAg       >  2:3951634/1‑100 (MQ=255)
                                                                                               gCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGt      <  1:206412/100‑1 (MQ=255)
                                                                                                  gCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTAat   <  1:2299719/100‑1 (MQ=255)
                                                                                                   cAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTAata  >  2:1358185/1‑100 (MQ=255)
                                                                                                   cAGAACAAGACGGATAAAAAATACAATGAGAGCACACTCAAACAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTAata  <  2:833796/100‑1 (MQ=255)
                                                                                                   cAGAACAAGACGGAAAAAAAAACCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTa     <  2:541679/97‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTAATA  >  NC_000913/3769988‑3770186

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: