Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,068 T→C 100% pseudogene (300/528 nt) nohQ ← pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,0680TC99.0% 7112.3 / 19.1 3537pseudogene (300/528 nt)nohQpseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base C (1699/1796);  total (1719/1818)
Fisher's exact test for biased strand distribution p-value = 8.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 3537 total aligned reads displayed.

TTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTT  >  NC_000913/1635935‑1636213
                                                                                                                                     |                                                                                                                                                 
   ccAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTc                                                                                                                                  <  5:208180/148‑1 (MQ=12)
    cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGCATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg                                                                                                                                >  2:32412/1‑149 (MQ=17)
    cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCg                                                                                                                                 <  7:282669/148‑1 (MQ=12)
          cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACCGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                                                                                                                          >  5:285965/1‑149 (MQ=255)
          cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                                                                                                                          >  2:25923/1‑149 (MQ=12)
           ggAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGTGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                                                                                                                          >  1:40892/1‑148 (MQ=255)
              aaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcg                                                                                                                      <  5:178019/149‑1 (MQ=12)
               aaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                                                                                                                     <  1:126354/149‑1 (MQ=12)
               aaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                                                                                                                     <  2:208002/149‑1 (MQ=12)
                  gccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGt                                                                                                                   >  2:177662/1‑148 (MQ=12)
                         cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGa                                                                                                            <  5:126021/148‑1 (MQ=12)
                         cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGAt                                                                                                           <  5:260256/149‑1 (MQ=12)
                            cGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGTGTCGGCCTGCGCACGCGTAAGTCGATGGc                                                                                                        <  7:137972/149‑1 (MQ=255)
                            cGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc                                                                                                        <  3:69482/149‑1 (MQ=12)
                            cGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc                                                                                                        <  6:159048/149‑1 (MQ=12)
                               cggggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCACTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGtt                                                                                                     <  2:348578/147‑1 (MQ=255)
                               cAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCCTTCTTCAGTTCCTGTGCGTCGGCCTGCGCCAGCGTAAGTCGCCGGCGtt                                                                                                     >  8:369760/1‑149 (MQ=11)
                                aGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTc                                                                                                    <  5:75824/149‑1 (MQ=12)
                                     gagcccGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACt                                                                                               <  4:319514/146‑1 (MQ=12)
                                     gagcccGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGAACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCTTAAGTCGATGGCGTTCGTACt                                                                                               <  1:279005/146‑1 (MQ=11)
                                        cccGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCg                                                                                                       >  4:252274/1‑138 (MQ=32)
                                        cccGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCaa                                                                                            >  8:344133/1‑149 (MQ=17)
                                        cccGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCa                                                                                             >  4:277550/1‑148 (MQ=17)
                                          cGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGTAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATa                                                                                          <  7:362701/149‑1 (MQ=12)
                                          cGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCTCACGCGTAAGTCGATGGCGTTCGTACTCaa                                                                                            >  8:229748/1‑147 (MQ=12)
                                               agaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTtcc                                                                                     <  3:98125/149‑1 (MQ=12)
                                               agaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTtcc                                                                                     >  5:101496/1‑149 (MQ=12)
                                               agaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTtcc                                                                                     <  3:143165/149‑1 (MQ=12)
                                               agaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCAGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTtcc                                                                                     >  2:148215/1‑149 (MQ=255)
                                                gaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCct                                                                                    >  7:189086/1‑149 (MQ=12)
                                                gaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCct                                                                                    <  6:208310/149‑1 (MQ=12)
                                                gaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACAACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCct                                                                                    >  2:116009/1‑149 (MQ=255)
                                                  aTACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCAGCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctgg                                                                                  <  6:232075/149‑1 (MQ=255)
                                                  aTACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACGACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctgg                                                                                  >  4:300487/1‑149 (MQ=255)
                                                   tACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggc                                                                                 >  2:47561/1‑149 (MQ=12)
                                                    aCTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggct                                                                                >  4:68781/1‑149 (MQ=12)
                                                      tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgc                                                                                                                           <  6:182407/104‑1 (MQ=12)
                                                      tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg                                                                                                          >  7:332522/1‑121 (MQ=17)
                                                      tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGttc                                                                                      >  3:26491/1‑139 (MQ=12)
                                                      tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctg                                                                               >  3:329179/1‑148 (MQ=12)
                                                      tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctg                                                                               >  3:110578/1‑148 (MQ=12)
                                                        gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGgag                                                                            <  8:197006/149‑1 (MQ=12)
                                                          aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGAt                                                                          <  8:101500/149‑1 (MQ=12)
                                                          aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGAt                                                                          <  5:369636/149‑1 (MQ=12)
                                                           aagtCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGTCATTCTTCAGTTCCTGTGCGTCTGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCTATAGTTCCTGGCTGGAGATc                                                                         <  4:42105/146‑1 (MQ=11)
                                                            tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcgc                                                                                                                         <  4:226612/100‑1 (MQ=16)
                                                            tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCt                                                                        <  6:317823/149‑1 (MQ=12)
                                                            tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCt                                                                        <  7:22115/149‑1 (MQ=12)
                                                             ttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTg                                                                       >  4:193147/1‑149 (MQ=12)
                                                               cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTgcc                                                                     >  6:133879/1‑147 (MQ=12)
                                                               cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTgc                                                                      <  3:282932/148‑2 (MQ=12)
                                                               cACCTGCGATCCGCGACAGCAAGAAAGTACAGAATGCGGTTTCCACCACTTGAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTGGGCCTGCGCACGCGTAAGTGGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTgc                                                                      <  7:104904/148‑2 (MQ=11)
                                                                aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGcct                                                                    <  2:379528/149‑1 (MQ=255)
                                                                aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCGCTGGCAGTCTGCAGTTCCTGTGCGGCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGcct                                                                    <  8:253442/149‑1 (MQ=11)
                                                                  cTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCg                                                                  >  1:333085/1‑149 (MQ=255)
                                                                    gCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCt                                                                >  3:351695/1‑149 (MQ=255)
                                                                     cGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGc                                                                 <  4:343685/147‑1 (MQ=32)
                                                                     cGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTg                                                               <  2:231617/149‑1 (MQ=255)
                                                                          cgcgACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAAGAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTg                                                          >  8:163848/1‑149 (MQ=11)
                                                                                  cACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAg                                                    <  1:199498/147‑1 (MQ=255)
                                                                                  cACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGtt                                                  >  7:308123/1‑149 (MQ=255)
                                                                                   aCGTAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCCCGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGAGCTGCCTCGCGGGCCTGCCGCAGttc                                                 >  8:288598/1‑149 (MQ=11)
                                                                                   aCGAAAGTACAGAATGCGGGTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                                                                                                                     <  3:44010/81‑1 (MQ=11)
                                                                                     gAAAGTACAGCATGCGGTTTCCACCCCTTCAGCGGAGTCTCTGGCATTCGTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTACTGGCTCGAGATCTGCCTCGCTGGCCTGCCGCAGttctt                                               >  1:69347/1‑149 (MQ=11)
                                                                                     gAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttctt                                               >  8:314208/1‑149 (MQ=255)
                                                                                     gAAAGTACAGAATGCGGTTTCCACCACGTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTGCTGAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttctt                                               >  6:80377/1‑149 (MQ=11)
                                                                                      aaaGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttcttc                                              <  7:16258/149‑1 (MQ=255)
                                                                                         gTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAAc                                           <  1:91389/149‑1 (MQ=255)
                                                                                          tACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTCCTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACt                                          >  6:107820/1‑149 (MQ=25)
                                                                                                aTGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttcttc                                              >  6:26034/1‑139 (MQ=255)
                                                                                                 tGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGa                                                                                                            <  3:228392/76‑1 (MQ=255)
                                                                                                     gTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAg                               >  5:230600/1‑149 (MQ=11)
                                                                                                      tttCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGc                              <  7:49788/149‑1 (MQ=11)
                                                                                                          caccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCtttt                          <  4:165031/149‑1 (MQ=11)
                                                                                                          caccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCtttt                          >  8:71410/1‑149 (MQ=11)
                                                                                                          caccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCtttt                          >  2:78325/1‑149 (MQ=11)
                                                                                                           accacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTc                         >  1:160407/1‑149 (MQ=11)
                                                                                                            ccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCg                        >  7:253902/1‑149 (MQ=11)
                                                                                                            ccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCg                        <  1:100552/149‑1 (MQ=11)
                                                                                                               cTTCCGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtc                     >  3:198499/1‑149 (MQ=11)
                                                                                                                ttCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtct                    >  6:45581/1‑149 (MQ=11)
                                                                                                                ttCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtct                    >  8:125149/1‑149 (MQ=11)
                                                                                                                  cAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAAGCGTAAGTCGATGGCGTTCGTACTCACTAGTTCATGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtct                    <  1:247340/147‑1 (MQ=11)
                                                                                                                      ggAGTCTCTGGCATTCTTCAGTTCCTGTGCGGCGGCCGGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTTGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAAttt              <  5:51773/149‑1 (MQ=11)
                                                                                                                        aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa            <  5:348463/149‑1 (MQ=11)
                                                                                                                        aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa            <  2:289372/149‑1 (MQ=11)
                                                                                                                        aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa            <  6:346454/149‑1 (MQ=11)
                                                                                                                        aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa            <  7:387573/149‑1 (MQ=11)
                                                                                                                          tctcTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGc          <  4:135081/149‑1 (MQ=11)
                                                                                                                          tctcTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGc          <  5:316990/149‑1 (MQ=11)
                                                                                                                            tctGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttctt                                               >  1:226562/1‑110 (MQ=11)
                                                                                                                            tctGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCAt        >  3:170710/1‑149 (MQ=11)
                                                                                                                            tctGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGAGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCAt        >  1:366027/1‑149 (MQ=11)
                                                                                                                             ctGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATc       >  6:370748/1‑149 (MQ=11)
                                                                                                                             ctGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATc       >  7:76821/1‑149 (MQ=11)
                                                                                                                              tGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttcttc                                              <  5:151361/109‑1 (MQ=11)
                                                                                                                                gCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATc       <  3:310797/146‑1 (MQ=11)
                                                                                                                                gCATTCTTCAGTTCCTGTGCGGCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATcc      <  6:286919/147‑1 (MQ=11)
                                                                                                                                  aTTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCttt  >  1:300057/1‑149 (MQ=11)
                                                                                                                                  aTTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCttt  <  8:46511/149‑1 (MQ=11)
                                                                                                                                     |                                                                                                                                                 
TTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTT  >  NC_000913/1635935‑1636213

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