Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,407 G→A pseudogene (143/669 nt) rph ← ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,4070GA99.2% 496.0 / ‑4.2 122pseudogene (143/669 nt)rphribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (50/70);  total (51/71)
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.38e-01

Only 100 of 123 total aligned reads displayed.

AGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTG  >  NC_000913/3816259‑3816553
                                                                                                                                                    |                                                                                                                                                   
     gggCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGcc                                                                                                                                              <  1:652355/151‑1 (MQ=255)
        cGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCcttc                                                                                                                                           <  1:1139898/151‑1 (MQ=255)
        cGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCcttc                                                                                                                                           <  2:1412788/151‑1 (MQ=255)
         gATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCcttct                                                                                                                                          <  1:790698/151‑1 (MQ=255)
         gATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCcttct                                                                                                                                          <  2:2103779/151‑1 (MQ=255)
          aTCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCATCACGCcttctt                                                                                                                                         <  1:1363337/151‑1 (MQ=255)
           tCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCttcttc                                                                                                                                        >  1:2012444/1‑151 (MQ=255)
           tCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCttcttc                                                                                                                                        >  1:825890/1‑151 (MQ=255)
             aGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCaa                                                                                                                                      <  1:379787/151‑1 (MQ=255)
             aGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCa                                                                                                                                       >  1:1937692/1‑150 (MQ=255)
             aGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGAATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCa                                                                                                                                       >  2:1798504/1‑150 (MQ=255)
                cGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATag                                                                                                                                   <  2:1133251/151‑1 (MQ=255)
                  cTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATaga                                                                                                                                  <  2:1664929/150‑1 (MQ=255)
                  cTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATag                                                                                                                                   >  2:1757304/1‑149 (MQ=255)
                    ggATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATagag                                                                                                                                 <  1:722527/149‑1 (MQ=255)
                    ggATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAgg                                                                                                                                <  2:875279/150‑1 (MQ=255)
                    ggATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGc                                                                                                                               <  1:1616749/151‑1 (MQ=255)
                    ggATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGc                                                                                                                               <  2:937095/151‑1 (MQ=255)
                      aTTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCgg                                                                                                                            <  1:1552479/151‑1 (MQ=255)
                        ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGt                                                                                                                           >  1:681021/1‑150 (MQ=255)
                        ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGAATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGt                                                                                                                           <  1:1798538/150‑1 (MQ=255)
                          ccATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTaca                                                                                                                        <  1:1357275/151‑1 (MQ=255)
                           cATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTacac                                                                                                                       >  1:195087/1‑151 (MQ=255)
                               gtgCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACacaaca                                                                                                                   <  1:794498/151‑1 (MQ=255)
                                 gCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAacac                                                                                                                  <  1:1301232/150‑1 (MQ=255)
                                 gCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACt                                                                                                                 <  2:1012550/151‑1 (MQ=255)
                                    tCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACttt                                                                                                               >  1:332451/1‑150 (MQ=255)
                                     ccACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTgg                                                                                                             <  1:1010897/151‑1 (MQ=255)
                                     ccACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTg                                                                                                              <  1:457477/150‑1 (MQ=255)
                                      cACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACttt                                                                                                               >  1:949277/1‑148 (MQ=255)
                                       aCCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTgg                                                                                                             >  2:173071/1‑149 (MQ=255)
                                        cccTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTAt                                                                                                          <  2:1452292/151‑1 (MQ=255)
                                         ccTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTAt                                                                                                          <  1:173801/150‑1 (MQ=255)
                                           tGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGc                                                                                                       >  2:1905768/1‑151 (MQ=255)
                                              ttACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCa                                                                                                     >  1:1765692/1‑150 (MQ=255)
                                               tACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCa                                                                                                     <  2:673650/149‑1 (MQ=255)
                                               tACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCa                                                                                                     <  2:1485801/149‑1 (MQ=255)
                                                aCCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCaaa                                                                                                   <  2:121992/150‑1 (MQ=255)
                                                 ccTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCTTGCCCTGAACTTTTCGGAAGCGCCGCCCGCCTTTTTCAATAAAAGCGGGTAACCAAACCTTTGGATTCGCCAAAt                                                                                                  >  2:1974337/1‑150 (MQ=255)
                                                   tttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTc                                                                                                >  2:1252460/1‑150 (MQ=255)
                                                   tttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCg                                                                                               >  2:689167/1‑151 (MQ=255)
                                                   tttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCg                                                                                               >  2:1613795/1‑151 (MQ=255)
                                                    ttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCg                                                                                               >  1:1781899/1‑150 (MQ=255)
                                                       gccgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGAcc                                                                                            >  1:238435/1‑150 (MQ=255)
                                                        ccgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGAcc                                                                                            >  1:1117128/1‑149 (MQ=255)
                                                        ccgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAg                                                                                          <  2:1531321/151‑1 (MQ=255)
                                                         cgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         <  2:863377/151‑1 (MQ=255)
                                                         cgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         <  2:2063818/151‑1 (MQ=255)
                                                         cgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         <  2:647595/151‑1 (MQ=255)
                                                          gcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         >  2:412437/1‑150 (MQ=255)
                                                            ttCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCAcc                                                                                      >  1:1551297/1‑151 (MQ=255)
                                                             tCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCAc                                                                                       <  1:249911/149‑1 (MQ=255)
                                                              cACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCg                                                                                     >  1:875464/1‑150 (MQ=255)
                                                              cACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGa                                                                                    <  2:445235/151‑1 (MQ=255)
                                                                   gCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCtt                                                                               >  1:900577/1‑151 (MQ=255)
                                                                     gTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCt                                                                             <  1:22113/151‑1 (MQ=255)
                                                                     gTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCt                                                                             <  2:1317460/151‑1 (MQ=255)
                                                                         cGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCAt                                                                         >  1:856401/1‑151 (MQ=255)
                                                                            gtgtGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                      >  2:1420324/1‑151 (MQ=255)
                                                                            gtgtGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                      <  2:1924299/151‑1 (MQ=255)
                                                                             tgtgGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGttt                                                                     >  2:1780606/1‑151 (MQ=255)
                                                                                    aGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAg                                                              <  2:1512994/151‑1 (MQ=255)
                                                                                      aaCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAg                                                              >  1:1310403/1‑149 (MQ=255)
                                                                                      aaCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGtt                                                            <  2:367821/151‑1 (MQ=255)
                                                                                      aaCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGtt                                                            >  2:1730362/1‑151 (MQ=255)
                                                                                       aCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGtt                                                            >  2:1438725/1‑150 (MQ=255)
                                                                                           ggCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGa                                                        <  2:2000899/150‑1 (MQ=255)
                                                                                            gcagcaTGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGt                                                      <  2:976898/151‑1 (MQ=255)
                                                                                                 aTGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAggg                                                 >  1:515460/1‑151 (MQ=255)
                                                                                                   gCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTGTTGTATAGTTACGAGGCa                                                    >  1:193544/1‑146 (MQ=255)
                                                                                                    ccGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGt                                                >  1:838906/1‑149 (MQ=255)
                                                                                                    ccGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTaa                                              <  2:348198/151‑1 (MQ=255)
                                                                                                     cGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTaa                                              >  2:648423/1‑150 (MQ=255)
                                                                                                     cGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAAc                                             >  2:1078036/1‑151 (MQ=255)
                                                                                                      gTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                            <  2:1080275/151‑1 (MQ=255)
                                                                                                      gTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                            <  2:2036708/151‑1 (MQ=255)
                                                                                                       tACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACgg                                           >  1:1398053/1‑151 (MQ=255)
                                                                                                       tACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACgg                                           >  2:1842200/1‑151 (MQ=255)
                                                                                                           ctGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGa                                         <  2:2044515/149‑1 (MQ=255)
                                                                                                           ctGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACg                                       >  1:1782132/1‑151 (MQ=255)
                                                                                                             gCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCa                                     <  2:2012401/151‑1 (MQ=255)
                                                                                                                 tGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTg                                 >  1:1810704/1‑151 (MQ=255)
                                                                                                                   aTCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                               <  2:333854/151‑1 (MQ=255)
                                                                                                                     ccAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGatt                              <  1:606897/150‑1 (MQ=255)
                                                                                                                       aGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGattatt                           >  1:1647046/1‑151 (MQ=255)
                                                                                                                       aGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGattatt                           <  2:1361544/151‑1 (MQ=255)
                                                                                                                         ccctggccctgACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGc                         <  2:1678225/151‑1 (MQ=255)
                                                                                                                             ggccctgACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACg                     <  1:1121471/151‑1 (MQ=255)
                                                                                                                                cctgACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCTACGGcc                  >  1:2015608/1‑151 (MQ=255)
                                                                                                                                   gACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctgc               <  1:1842234/151‑1 (MQ=255)
                                                                                                                                   gACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctg                <  2:103954/150‑1 (MQ=255)
                                                                                                                                    aCCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctgct              >  1:1166380/1‑151 (MQ=255)
                                                                                                                                     ccTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctgctg             <  1:615799/151‑1 (MQ=255)
                                                                                                                                       tttCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCTACGGCCTGCTgg            >  1:37221/1‑150 (MQ=255)
                                                                                                                                        ttCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGAc          >  2:1215737/1‑151 (MQ=255)
                                                                                                                                          cAGGAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACg         <  2:2075835/150‑1 (MQ=255)
                                                                                                                                            ggAAGCGCAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCat      <  2:1280320/151‑1 (MQ=255)
                                                                                                                                                gcgcAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTg  <  1:1921029/151‑1 (MQ=255)
                                                                                                                                                gcgcAGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTg  <  1:736928/151‑1 (MQ=255)
                                                                                                                                                 cgcaGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTg  >  2:1942862/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTG  >  NC_000913/3816259‑3816553

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