Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,770,087 C→A 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,0870CA94.7% 419.9 / 9.5 132pseudogene (55/114 nt)yibUpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (60/65);  total (65/67)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 133 total aligned reads displayed.

GAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTA  >  NC_000913/3769991‑3770183
                                                                                                |                                                                                                
gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                               <  2:126220/100‑1 (MQ=255)
gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                               >  2:2870942/1‑100 (MQ=255)
 aaaaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGgaagaa                                                                                              <  2:2831765/100‑1 (MQ=255)
   aaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCGTCCTGCGGAAGAACa                                                                                            <  2:733865/100‑1 (MQ=255)
   aaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTCCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACa                                                                                            >  1:3262951/1‑100 (MQ=255)
    aGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACaa                                                                                           <  2:3464698/100‑1 (MQ=255)
      cACTCCTGCTCTTTTTTTTTTATATTCGTCGCAGGAGGAACAACTCATAACATCATTGCGATCTCGACGATGTATAGTCTCATTCTGCGGAAGCAcca                                                                                           >  1:1988563/1‑96 (MQ=255)
       aCTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGAc                                                                                        <  1:3220842/100‑1 (MQ=255)
         tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACgg                                                                                      <  1:1853129/100‑1 (MQ=255)
         tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACgg                                                                                      <  2:1452761/100‑1 (MQ=255)
            tGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATa                                                                                   <  1:1946631/100‑1 (MQ=255)
             gCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaa                                                                                  >  2:227882/1‑100 (MQ=255)
               tctTTTTTGATTATATCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaaaa                                                                                >  1:156055/1‑100 (MQ=255)
                ctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaaaaa                                                                               >  1:716367/1‑100 (MQ=255)
                  tttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAAt                                                                             <  2:1900795/100‑1 (MQ=255)
                    tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATcc                                                                           >  2:616733/1‑100 (MQ=255)
                       gATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaan                                                                        <  2:3008028/100‑2 (MQ=255)
                       gATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAAt                                                                        <  2:3187500/100‑1 (MQ=255)
                       ccgatctgCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAAt                                                                        <  2:1493512/93‑1 (MQ=255)
                       ccgatctgCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAAt                                                                        <  2:1493395/93‑1 (MQ=255)
                        cgatctgCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATg                                                                       <  1:3138222/94‑1 (MQ=255)
                        aatAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATATCGACGATGTATAGTCTCATCCTGCGGAAGAAAAAGAAAAAAAAAAAATcc                                                                           <  1:1091593/94‑1 (MQ=255)
                          tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                     <  2:513711/100‑1 (MQ=255)
                          tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                     >  2:3094423/1‑100 (MQ=255)
                           agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGTCGATGTATAGTCTCATCCTACAGAACAACAAGACGGCTAAAAAATCCAATgaga                                                                    >  1:2661473/1‑100 (MQ=255)
                            tagacGTCGCAGGAGTAACAACTCATATCATCGTTGCGATCGCGACGATGTATAGTCTCATCCTGCGGAATAACAACTCGGATAAAAAATCCAATGagag                                                                   <  1:2602609/96‑1 (MQ=255)
                            gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                   <  1:2085880/100‑1 (MQ=255)
                              gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAAtagatcgg                                                                 >  1:1493490/1‑93 (MQ=255)
                              gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgagatcg                                                                 >  2:3138272/1‑97 (MQ=255)
                              gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGca                                                                 >  2:2642928/1‑100 (MQ=255)
                              gCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGca                                                                 >  1:2996400/1‑100 (MQ=255)
                                 gTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCacac                                                              <  2:3351383/100‑1 (MQ=255)
                                   cGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCacacgt                                                            >  2:213886/1‑98 (MQ=255)
                                   cGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTc                                                            >  2:2935631/1‑100 (MQ=255)
                                    gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                           >  1:1055138/1‑100 (MQ=255)
                                     cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAt                                                          >  2:1270957/1‑100 (MQ=255)
                                     cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAt                                                          >  2:1363992/1‑100 (MQ=255)
                                       ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                        >  1:3382507/1‑100 (MQ=255)
                                       ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCacacgtctga                                                        >  2:1130594/1‑94 (MQ=255)
                                          ggaACAACTCATATCATCATTGCGATCTCGACGATGTGTAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                     >  1:21454/1‑100 (MQ=255)
                                            aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                   >  2:1910774/1‑100 (MQ=255)
                                            aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                   <  1:436855/100‑1 (MQ=255)
                                             acaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACt                                                  <  1:379999/100‑1 (MQ=255)
                                              caacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtt                                                 <  1:926559/100‑1 (MQ=255)
                                              caacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtt                                                 <  1:3444656/100‑1 (MQ=255)
                                              caacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtt                                                 >  1:2818934/1‑100 (MQ=255)
                                               aacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                >  1:581609/1‑100 (MQ=255)
                                               aacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGACAAAAAATCCAATGAGAGCACACTCATCCAGAACtta                                                >  2:2436451/1‑100 (MQ=255)
                                                acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttat                                               >  2:573590/1‑100 (MQ=255)
                                                acTCATATCATCATTGCGATCACAACGATGTATAGTCTCATCCTGCGGAAGAAAAAAAAAAAAAAAAAATCCAATGAGAGCACACTCATCCAGAACttat                                               <  1:2374730/100‑1 (MQ=255)
                                                 cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                              <  2:312547/100‑1 (MQ=255)
                                                  tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                             <  2:878918/100‑1 (MQ=255)
                                                    atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                           <  1:1529125/100‑1 (MQ=255)
                                                    atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAAc                                           >  2:2126752/1‑100 (MQ=255)
                                                     tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                          <  1:1910745/100‑1 (MQ=255)
                                                     tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                          >  2:1667886/1‑100 (MQ=255)
                                                       tcaacaTTGCGATCTCGACTATGTATAGTCTCATCCTTCGGCAGAAGAAAACCGATAAAAAATCCAAGGAGAGCACACTCATCCAGAACTTATTAACaga                                        >  1:1156231/1‑98 (MQ=255)
                                                        catcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATaa                                       >  1:46164/1‑100 (MQ=255)
                                                        catcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATaa                                       >  1:1460231/1‑100 (MQ=255)
                                                          tcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACTCTCATCCAGAACTTATTAACATAATg                                     <  2:2113435/100‑1 (MQ=255)
                                                           catTGCAATCTCAACAATGTATAGTCTCATCCTGCGGAACAAAAAAAAAAAAAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATg                                     <  1:710856/99‑1 (MQ=255)
                                                            atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                   >  2:1084933/1‑100 (MQ=255)
                                                             ttGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAat                                  <  2:1741477/100‑1 (MQ=255)
                                                              tGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAata                                 >  1:1445221/1‑100 (MQ=255)
                                                              tGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAata                                 <  2:1319455/100‑1 (MQ=255)
                                                               gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                                >  1:354230/1‑100 (MQ=255)
                                                               gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                                >  1:3328191/1‑100 (MQ=255)
                                                               gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                                >  1:3073856/1‑100 (MQ=255)
                                                                cGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATAtt                               >  2:1736378/1‑100 (MQ=255)
                                                                cGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATAtt                               <  2:1014132/100‑1 (MQ=255)
                                                                  aTCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaa                             >  1:1981532/1‑100 (MQ=255)
                                                                  aTCTAGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACAGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaa                             <  1:1659476/100‑1 (MQ=255)
                                                                    ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                           <  1:2625528/100‑1 (MQ=255)
                                                                    ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                           <  2:448608/100‑1 (MQ=255)
                                                                    ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                           <  1:83309/100‑1 (MQ=255)
                                                                     tcGACGATGTATAGTCTCATCCTGCGGAAGCACAAGACGGATAAAAAATCCAATGAGCGAACACTCATCCAGAACTTATTAACATAATGAATATTAAAGa                          >  1:1193194/1‑100 (MQ=255)
                                                                     tcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGa                          >  1:1377099/1‑100 (MQ=255)
                                                                      cgacgaTGTATAGTCTCATCCTGAGGAAGAACAAGACGGATAAAAAATCCAATGAGAGAACACTAATCAAGCACTTATTAACATAATGAATATAAAAGa                          >  1:2842945/1‑99 (MQ=255)
                                                                       gacgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaa                        >  2:421925/1‑100 (MQ=255)
                                                                         cgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaaa                      <  2:1019974/100‑1 (MQ=255)
                                                                           aggTATAGACTCATACTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatc                    <  1:3015598/98‑1 (MQ=255)
                                                                           aTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatc                    <  1:5251/100‑1 (MQ=255)
                                                                           aTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatc                    <  1:1084920/100‑1 (MQ=255)
                                                                            tGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                   >  1:664263/1‑100 (MQ=255)
                                                                              tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc                 <  2:210523/100‑1 (MQ=255)
                                                                               ataGTCGCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCAGACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatag                >  1:749136/1‑98 (MQ=255)
                                                                                  gTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg             <  2:2053054/100‑1 (MQ=255)
                                                                                  gTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg             >  1:1365798/1‑100 (MQ=255)
                                                                                    ctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCt           <  1:1406082/100‑1 (MQ=255)
                                                                                      cATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTaa         >  2:2951382/1‑100 (MQ=255)
                                                                                       aaccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATACAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAAt        <  1:1112779/98‑1 (MQ=255)
                                                                                       aTCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAAt        <  1:2780816/100‑1 (MQ=255)
                                                                                        tCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATg       <  2:2254757/100‑1 (MQ=255)
                                                                                         ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATgn      <  1:1028567/100‑2 (MQ=255)
                                                                                         ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa      <  2:351787/100‑1 (MQ=255)
                                                                                         ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa      <  2:191352/100‑1 (MQ=255)
                                                                                         ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa      >  2:1353924/1‑100 (MQ=255)
                                                                                          cTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa     >  2:3113006/1‑100 (MQ=255)
                                                                                          cTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa     >  1:1467344/1‑100 (MQ=255)
                                                                                           tGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAg    <  2:20304/100‑1 (MQ=255)
                                                                                                |                                                                                                
GAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGTA  >  NC_000913/3769991‑3770183

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