Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,575,931 A→G pseudogene (2013/2937 nt) yjiV → pseudogene; conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,575,9310AG97.9% 519.5 / ‑2.5 142pseudogene (2013/2937 nt)yjiVpseudogene; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (66/73);  total (68/74)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 142 total aligned reads displayed.

AATGCCGGAAATCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAAGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATCGCGT  >  NC_000913/4575784‑4576077
                                                                                                                                                   |                                                                                                                                                  
aaTGCCGGAAATCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGggc                                                                                                                                                 >  2:137365/1‑151 (MQ=255)
  tGCCGGAAATCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGggcgg                                                                                                                                               >  1:59936/1‑151 (MQ=255)
      ggAAATCTGGAAGCTGGATCATTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTc                                                                                                                                           <  4:1622037/151‑1 (MQ=255)
       gAAATCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTct                                                                                                                                          <  3:365084/151‑1 (MQ=255)
       gAAATCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTct                                                                                                                                          >  3:82073/1‑151 (MQ=255)
        aaaTCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTct                                                                                                                                          >  4:937592/1‑150 (MQ=255)
         aaTCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTctac                                                                                                                                        >  3:251649/1‑151 (MQ=255)
         aaTCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTctac                                                                                                                                        <  3:336369/151‑1 (MQ=255)
         aaTCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTctac                                                                                                                                        >  2:355743/1‑151 (MQ=255)
           tCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTctact                                                                                                                                       >  1:137544/1‑150 (MQ=255)
           tCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTctact                                                                                                                                       >  2:323056/1‑150 (MQ=255)
              ggAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCg                                                                                                                                   >  2:856114/1‑151 (MQ=255)
               gAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCGCAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCGTCTACTATc                                                                                                                                    <  2:376158/149‑1 (MQ=255)
               gAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGt                                                                                                                                  >  2:353233/1‑151 (MQ=255)
               gAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGt                                                                                                                                  <  4:1372838/151‑1 (MQ=255)
               gAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGt                                                                                                                                  <  4:1441998/151‑1 (MQ=255)
                 aGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTAc                                                                                                                                >  4:1301770/1‑151 (MQ=255)
                 aGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTAc                                                                                                                                >  2:1439435/1‑151 (MQ=255)
                         cGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGc                                                                                                                        >  1:399506/1‑151 (MQ=255)
                           ttCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCAc                                                                                                                      >  4:1373447/1‑151 (MQ=255)
                             cacaGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCActc                                                                                                                    >  2:430230/1‑151 (MQ=255)
                              acaGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCActct                                                                                                                   <  1:399737/151‑1 (MQ=255)
                               gaggacgagaaCGCAGAGGGTGGGGAAAAAGGAGGCGACTGGGTCAGACATGGTGGGCGGACCGCACATTTAAAAAAGAAAAACGTCTGGACAGACATTAACGATCGAGTCGTGGAAGGCGGCTTCTCCTATCGTCCTGCCGCGCAct                                                                                                                     >  2:775093/11‑148 (MQ=255)
                                    aGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCcag                                                                                                             <  4:23515/151‑1 (MQ=255)
                                    aGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCcag                                                                                                             <  2:629087/151‑1 (MQ=255)
                                       cTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCcagc                                                                                                            >  1:89121/1‑149 (MQ=255)
                                        tACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAgt                                                                                                         >  4:1171233/1‑151 (MQ=255)
                                         aCGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAgtc                                                                                                        >  4:189141/1‑151 (MQ=255)
                                         aCGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAgtc                                                                                                        <  3:283850/151‑1 (MQ=255)
                                              gagGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAgtcgtc                                                                                                     >  2:69470/1‑149 (MQ=255)
                                                 ggTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAAc                                                                                                >  2:1365167/1‑151 (MQ=255)
                                                   tCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGt                                                                                              >  4:685028/1‑151 (MQ=255)
                                                         aaaaGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCaa                                                                                        <  2:1077464/151‑1 (MQ=255)
                                                          aaaGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCaaa                                                                                       <  4:734783/151‑1 (MQ=255)
                                                             gCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCt                                                                                    <  2:1080390/151‑1 (MQ=255)
                                                             gCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGCCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCt                                                                                    <  3:398724/151‑1 (MQ=255)
                                                              cACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTa                                                                                   <  4:794220/151‑1 (MQ=255)
                                                               aCGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTAt                                                                                  <  4:531503/151‑1 (MQ=255)
                                                               aCGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTAt                                                                                  <  2:1543692/151‑1 (MQ=255)
                                                                cgcgACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATg                                                                                 <  4:908299/151‑1 (MQ=255)
                                                                 gcgACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATg                                                                                 >  3:61526/1‑150 (MQ=255)
                                                                        ggTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATg                                                                         >  3:403597/1‑151 (MQ=255)
                                                                         gTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGt                                                                        <  3:92482/151‑1 (MQ=255)
                                                                         gTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGt                                                                        <  2:965286/151‑1 (MQ=255)
                                                                              ccATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAg                                                                   <  1:229453/151‑1 (MQ=255)
                                                                              ccATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAg                                                                   <  4:1277801/151‑1 (MQ=255)
                                                                               cATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGa                                                                  <  4:1683239/151‑1 (MQ=255)
                                                                               cATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGa                                                                  <  2:857618/151‑1 (MQ=255)
                                                                                  gATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAAt                                                                <  1:108274/150‑1 (MQ=255)
                                                                                   aTCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATg                                                               <  2:592755/150‑1 (MQ=255)
                                                                                    tCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATgg                                                              <  4:857234/150‑1 (MQ=255)
                                                                                      cGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCa                                                           <  4:137596/151‑1 (MQ=255)
                                                                                      cGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCa                                                           <  4:1129993/151‑1 (MQ=255)
                                                                                      cGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCa                                                           <  4:1587780/151‑1 (MQ=255)
                                                                                       gCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCaa                                                          <  4:1109157/151‑1 (MQ=255)
                                                                                       gCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCaa                                                          <  1:359896/151‑1 (MQ=255)
                                                                                        cTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAAc                                                         >  2:220414/1‑151 (MQ=255)
                                                                                           aTCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTa                                                       >  2:811994/1‑150 (MQ=255)
                                                                                              gCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATg                                                   >  1:73580/1‑151 (MQ=255)
                                                                                               cacaATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGt                                                  <  4:1659376/151‑1 (MQ=255)
                                                                                                 caATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTtc                                                >  2:536642/1‑151 (MQ=255)
                                                                                                   aTTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTtctc                                              <  2:1482414/151‑1 (MQ=255)
                                                                                                   aTTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTtctc                                              >  2:4123/1‑151 (MQ=255)
                                                                                                      gCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCaa                                            <  2:680164/150‑1 (MQ=255)
                                                                                                       cGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACt                                          >  2:1124683/1‑151 (MQ=255)
                                                                                                           cGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGtt                                       <  3:98973/150‑1 (MQ=255)
                                                                                                           cGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTc                                      <  3:51277/151‑1 (MQ=255)
                                                                                                           cGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTc                                      <  4:1427926/151‑1 (MQ=255)
                                                                                                           cGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTc                                      <  4:46586/151‑1 (MQ=255)
                                                                                                             gAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTcca                                    <  4:584938/151‑1 (MQ=255)
                                                                                                              aaaaCGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTccac                                   <  1:106357/151‑1 (MQ=255)
                                                                                                              aaaaCGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTccac                                   >  2:54638/1‑151 (MQ=255)
                                                                                                                aaCGTCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTccacca                                 <  4:430831/151‑1 (MQ=255)
                                                                                                                    tCTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCAt                             >  4:1027230/1‑151 (MQ=255)
                                                                                                                     cTGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATg                            >  2:1503885/1‑151 (MQ=255)
                                                                                                                      tGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATgg                           >  2:226606/1‑151 (MQ=255)
                                                                                                                      tGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATgg                           >  1:188990/1‑151 (MQ=255)
                                                                                                                      tGGACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATgg                           >  1:245201/1‑151 (MQ=255)
                                                                                                                       ggACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCAt                             <  2:948451/148‑1 (MQ=255)
                                                                                                                       ggACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCAt                             <  2:1485888/148‑1 (MQ=255)
                                                                                                                       ggACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCAt                             <  4:208286/148‑1 (MQ=255)
                                                                                                                       ggACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATgg                           <  2:1178582/150‑1 (MQ=255)
                                                                                                                       ggACAGACATTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATg                            <  4:634550/149‑1 (MQ=255)
                                                                                                                        gacagacaTTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGaa                         <  1:202499/151‑1 (MQ=255)
                                                                                                                          cagacaTTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGaaa                        <  4:1470839/150‑1 (MQ=255)
                                                                                                                           agacaTTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATg                      >  2:937982/1‑151 (MQ=255)
                                                                                                                              caTTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATgggg                   <  4:1437484/151‑1 (MQ=255)
                                                                                                                               aTTAACGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATggggg                  <  1:335056/151‑1 (MQ=255)
                                                                                                                                  aaCGATCGTGTCGTGGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGtt                <  2:1413311/150‑1 (MQ=255)
                                                                                                                                        cgtgtcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCg         <  4:1192698/151‑1 (MQ=255)
                                                                                                                                        cgtgtcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCg         <  4:1188268/151‑1 (MQ=255)
                                                                                                                                        cgtgtcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCg         >  1:157381/1‑151 (MQ=255)
                                                                                                                                          tgtcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGAt       >  4:1675464/1‑151 (MQ=255)
                                                                                                                                           gtcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGAt       <  2:188830/150‑1 (MQ=255)
                                                                                                                                           gtcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATc      <  4:10393/151‑1 (MQ=255)
                                                                                                                                            tcgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATcg     >  2:499143/1‑151 (MQ=255)
                                                                                                                                             cgtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATcg     <  3:25059/150‑1 (MQ=255)
                                                                                                                                              gtgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTTTGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATcgcg   >  2:129367/1‑151 (MQ=255)
                                                                                                                                              gtgGAGGGCGGCTTCTACCATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATcgcg   >  2:686506/1‑151 (MQ=255)
                                                                                                                                               tgGAGGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATCGCGt  >  2:154499/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATGCCGGAAATCTGGAAGCTGGATCGTTCACAGGAAGACTACGCAGAGGGTCTGGCAAAAGCACGCGACTGGGTCAGCCATGATCCGCTGATCGCACAATTGCGGTCGGAAAACGTCTGGACAGACATTAACGATCGTGTCGTGGAAGGCGGCTTCTACTATCGTACTGCCGAGCACTCTGCCCAGCAGTCGTCAGAACGTCTGCAAAGCTATGAAAAGATGTTCAAGAATGGGCAACTAAATGTTCTCAACTGTTCCACCACCATGGAAATGGGGGTTGATACCGATCGCGT  >  NC_000913/4575784‑4576077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: