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.

AGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTG  >  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)
                      aTTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑gg                                                                                                                            <  1:1552479/151‑1 (MQ=255)
                        ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGt                                                                                                                           >  1:681021/1‑150 (MQ=255)
                        ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGAATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGt                                                                                                                           <  1:1798538/150‑1 (MQ=255)
                          ccATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTaca                                                                                                                        <  1:1357275/151‑1 (MQ=255)
                           cATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTacac                                                                                                                       >  1:195087/1‑151 (MQ=255)
                               gtgCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACacaaca                                                                                                                   <  1:794498/151‑1 (MQ=255)
                                 gCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAacac                                                                                                                  <  1:1301232/150‑1 (MQ=255)
                                 gCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACt                                                                                                                 <  2:1012550/151‑1 (MQ=255)
                                    tCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACttt                                                                                                               >  1:332451/1‑150 (MQ=255)
                                     ccACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTgg                                                                                                             <  1:1010897/151‑1 (MQ=255)
                                     ccACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTg                                                                                                              <  1:457477/150‑1 (MQ=255)
                                      cACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACttt                                                                                                               >  1:949277/1‑148 (MQ=255)
                                       aCCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTgg                                                                                                             >  2:173071/1‑149 (MQ=255)
                                        cccTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTAt                                                                                                          <  2:1452292/151‑1 (MQ=255)
                                         ccTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTAt                                                                                                          <  1:173801/150‑1 (MQ=255)
                                           tGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGc                                                                                                       >  2:1905768/1‑151 (MQ=255)
                                              ttACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCa                                                                                                     >  1:1765692/1‑150 (MQ=255)
                                               tACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCa                                                                                                     <  2:673650/149‑1 (MQ=255)
                                               tACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCa                                                                                                     <  2:1485801/149‑1 (MQ=255)
                                                aCCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCaaa                                                                                                   <  2:121992/150‑1 (MQ=255)
                                                 ccTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAG‑CCTTGCCCTGAACTTTTCGGAAGCGCCGCCCGCCTTTTTCAATAAAAGCGGGTAACCAAACCTTTGGATTCGCCAAAt                                                                                                  >  2:1974337/1‑150 (MQ=255)
                                                   tttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTc                                                                                                >  2:1252460/1‑150 (MQ=255)
                                                   tttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCg                                                                                               >  2:689167/1‑151 (MQ=255)
                                                   tttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCg                                                                                               >  2:1613795/1‑151 (MQ=255)
                                                    ttCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCg                                                                                               >  1:1781899/1‑150 (MQ=255)
                                                       gccgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGAcc                                                                                            >  1:238435/1‑150 (MQ=255)
                                                        ccgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGAcc                                                                                            >  1:1117128/1‑149 (MQ=255)
                                                        ccgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAg                                                                                          <  2:1531321/151‑1 (MQ=255)
                                                         cgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         <  2:863377/151‑1 (MQ=255)
                                                         cgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         <  2:2063818/151‑1 (MQ=255)
                                                         cgcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         <  2:647595/151‑1 (MQ=255)
                                                          gcTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGc                                                                                         >  2:412437/1‑150 (MQ=255)
                                                            ttCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCAcc                                                                                      >  1:1551297/1‑151 (MQ=255)
                                                             tCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCAc                                                                                       <  1:249911/149‑1 (MQ=255)
                                                              cACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCg                                                                                     >  1:875464/1‑150 (MQ=255)
                                                              cACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGa                                                                                    <  2:445235/151‑1 (MQ=255)
                                                                   gCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCtt                                                                               >  1:900577/1‑151 (MQ=255)
                                                                     gTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCt                                                                             <  1:22113/151‑1 (MQ=255)
                                                                     gTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCt                                                                             <  2:1317460/151‑1 (MQ=255)
                                                                         cGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCAt                                                                         >  1:856401/1‑151 (MQ=255)
                                                                            gtgtGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                      >  2:1420324/1‑151 (MQ=255)
                                                                            gtgtGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                      <  2:1924299/151‑1 (MQ=255)
                                                                             tgtgGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGttt                                                                     >  2:1780606/1‑151 (MQ=255)
                                                                                    aGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAg                                                              <  2:1512994/151‑1 (MQ=255)
                                                                                      aaCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAg                                                              >  1:1310403/1‑149 (MQ=255)
                                                                                      aaCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGtt                                                            <  2:367821/151‑1 (MQ=255)
                                                                                      aaCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGtt                                                            >  2:1730362/1‑151 (MQ=255)
                                                                                       aCGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGtt                                                            >  2:1438725/1‑150 (MQ=255)
                                                                                           ggCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGa                                                        <  2:2000899/150‑1 (MQ=255)
                                                                                            gcagcaTGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGt                                                      <  2:976898/151‑1 (MQ=255)
                                                                                                 aTGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAggg                                                 >  1:515460/1‑151 (MQ=255)
                                                                                                   gCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTGTTGTATAGTTACGAGGCa                                                    >  1:193544/1‑146 (MQ=255)
                                                                                                    ccGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGt                                                >  1:838906/1‑149 (MQ=255)
                                                                                                    ccGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTaa                                              <  2:348198/151‑1 (MQ=255)
                                                                                                     cGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTaa                                              >  2:648423/1‑150 (MQ=255)
                                                                                                     cGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAAc                                             >  2:1078036/1‑151 (MQ=255)
                                                                                                      gTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                            <  2:1080275/151‑1 (MQ=255)
                                                                                                      gTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                            <  2:2036708/151‑1 (MQ=255)
                                                                                                       tACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACgg                                           >  1:1398053/1‑151 (MQ=255)
                                                                                                       tACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACgg                                           >  2:1842200/1‑151 (MQ=255)
                                                                                                           ctGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGa                                         <  2:2044515/149‑1 (MQ=255)
                                                                                                           ctGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACg                                       >  1:1782132/1‑151 (MQ=255)
                                                                                                             gCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCa                                     <  2:2012401/151‑1 (MQ=255)
                                                                                                                 tGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTg                                 >  1:1810704/1‑151 (MQ=255)
                                                                                                                   aTCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                               <  2:333854/151‑1 (MQ=255)
                                                                                                                     ccAGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGatt                              <  1:606897/150‑1 (MQ=255)
                                                                                                                       aGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGattatt                           >  1:1647046/1‑151 (MQ=255)
                                                                                                                       aGCCCTGGCCCTGACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGattatt                           <  2:1361544/151‑1 (MQ=255)
                                                                                                                         ccctggccctgACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGc                         <  2:1678225/151‑1 (MQ=255)
                                                                                                                             ggccctgACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACg                     <  1:1121471/151‑1 (MQ=255)
                                                                                                                                cctgACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCTACGGcc                  >  1:2015608/1‑151 (MQ=255)
                                                                                                                                   gACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctgc               <  1:1842234/151‑1 (MQ=255)
                                                                                                                                   gACCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctg                <  2:103954/150‑1 (MQ=255)
                                                                                                                                    aCCTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctgct              >  1:1166380/1‑151 (MQ=255)
                                                                                                                                     ccTTTCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCctgctg             <  1:615799/151‑1 (MQ=255)
                                                                                                                                       tttCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCTACGGCCTGCTgg            >  1:37221/1‑150 (MQ=255)
                                                                                                                                        ttCAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGAc          >  2:1215737/1‑151 (MQ=255)
                                                                                                                                          cAGGAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACg         <  2:2075835/150‑1 (MQ=255)
                                                                                                                                            ggAAGCGCAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCat      <  2:1280320/151‑1 (MQ=255)
                                                                                                                                                gcgcAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTg  <  1:1921029/151‑1 (MQ=255)
                                                                                                                                                gcgcAGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTg  <  1:736928/151‑1 (MQ=255)
                                                                                                                                                 cgcaGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTg  >  2:1942862/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGC‑GGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTG  >  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: ‑