Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,128 C→T 100% pseudogene (240/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,1280CT98.7% 5379.3 / 22.3 2653pseudogene (240/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 C (22/16);  new base T (1382/1224);  total (1411/1242)
Fisher's exact test for biased strand distribution p-value = 6.25e-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 2680 total aligned reads displayed.

AGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCAC  >  NC_000913/1635982‑1636270
                                                                                                                                                  |                                                                                                                                              
agaATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCCGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTtcc                                                                                                                                              <  2:268921/149‑1 (MQ=255)
  aaTACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctg                                                                                                                                            <  4:155690/149‑1 (MQ=12)
   aTACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctgg                                                                                                                                           <  8:188173/149‑1 (MQ=12)
   aTACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGGCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctgg                                                                                                                                           <  3:138053/149‑1 (MQ=255)
       tgccAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCCCCACTTCAGCGGAGTATCTGGCAGTCTTCAGGTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGGTCGTACTCAATAGTGActgcctg                                                                                                                                        >  2:219197/4‑148 (MQ=11)
       tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTAATCAATAGTTCctggctg                                                                                                                                        >  3:246190/1‑148 (MQ=255)
       tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAAAGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctg                                                                                                                                        >  1:156657/1‑148 (MQ=255)
         gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctgg                                                                                                                                       <  7:343342/147‑1 (MQ=12)
         gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctg                                                                                                                                        >  4:68265/1‑146 (MQ=12)
         gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGga                                                                                                                                      <  4:18965/148‑1 (MQ=12)
         gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGga                                                                                                                                      <  6:202575/148‑1 (MQ=12)
           aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctg                                                                                                                                        >  1:122339/1‑144 (MQ=12)
           aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGAt                                                                                                                                   <  3:59733/149‑1 (MQ=12)
            aagtCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGTCATTCTTCAGTTCCTGTGCGTCTGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCTATAGTTCCTGGCTGGAGATc                                                                                                                                  <  4:42105/146‑1 (MQ=11)
            aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATc                                                                                                                                  <  6:154464/149‑1 (MQ=12)
                cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTgc                                                                                                                               <  6:35669/148‑2 (MQ=12)
                 aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCGCTGGCAGTCTGCAGTTCCTGTGCGGCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGcct                                                                                                                             <  8:253442/149‑1 (MQ=11)
                  ccTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGcctc                                                                                                                            <  3:317252/149‑1 (MQ=255)
                   cTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctg                                                                                                                                        >  2:374640/1‑136 (MQ=12)
                    tGCGATCCGCGACAGCACCAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGc                                                                                                                          >  7:153642/1‑149 (MQ=11)
                     gCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCt                                                                                                                         >  5:333402/1‑149 (MQ=255)
                      cGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGc                                                                                                                          <  5:2747/147‑1 (MQ=255)
                      cGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTg                                                                                                                        >  1:79142/1‑149 (MQ=255)
                      cGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTg                                                                                                                        <  4:119092/149‑1 (MQ=255)
                      cGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTg                                                                                                                        <  6:277075/149‑1 (MQ=255)
                       gATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTgcc                                                                                                                              <  8:154175/142‑3 (MQ=12)
                       gATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGc                                                                                                                          <  2:104256/146‑1 (MQ=255)
                       gATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGGTCGTACTCAATAGTTCCTGGCGGGAGATCTGCCTCGCTgg                                                                                                                       >  2:131589/1‑149 (MQ=11)
                         tCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGTGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGcc                                                                                                                     <  6:224719/149‑1 (MQ=11)
                         tCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGcc                                                                                                                     >  5:360344/1‑149 (MQ=255)
                         tCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGACGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGcc                                                                                                                     >  2:43249/1‑149 (MQ=11)
                            gcgACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCt                                                                                                                                 >  8:339405/1‑134 (MQ=12)
                            gcgACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGc                                                                                                                  >  3:369541/1‑149 (MQ=255)
                             cgACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTgg                                                                                                                       <  3:342030/143‑1 (MQ=255)
                             cgACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTgg                                                                                                                       >  4:342030/1‑143 (MQ=255)
                                   cACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGtt                                                                                                           <  6:302848/149‑1 (MQ=255)
                                    aCGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttc                                                                                                          >  7:91659/1‑149 (MQ=255)
                                    aCGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttc                                                                                                          >  6:106561/1‑149 (MQ=255)
                                    aCGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGGTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttc                                                                                                          <  5:147789/149‑1 (MQ=11)
                                     cGAAAGTACAGAATGCGGTTTGCTCGACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGGGTCGGCGTGCGCTCGCGTTAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttct                                                                                                         <  3:168464/149‑1 (MQ=11)
                                        aaGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCctggctgg                                                                                                                                       >  8:71344/1‑116 (MQ=255)
                                        aaGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttcttc                                                                                                       >  6:381188/1‑148 (MQ=255)
                                          gTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAAc                                                                                                    <  7:260301/149‑1 (MQ=255)
                                          gTACAGAATGCGGTTTCCACCACTTCAGCGGAGGCGCTGGCAGTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAAc                                                                                                    <  4:297951/149‑1 (MQ=11)
                                           tACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAAcc                                                                                                   >  6:248074/1‑148 (MQ=255)
                                              aGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTcc                                                                                                >  1:215901/1‑149 (MQ=11)
                                              aGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTcc                                                                                                >  3:292980/1‑149 (MQ=11)
                                              aGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTcc                                                                                                >  3:203084/1‑149 (MQ=11)
                                                      ggttCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAg                                                                                        <  1:293446/147‑1 (MQ=11)
                                                         tCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCtt                                                                                     >  5:224561/1‑149 (MQ=11)
                                                         tCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCATGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCt                                                                                      <  3:90146/148‑1 (MQ=11)
                                                           caccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCtttt                                                                                   >  8:11959/1‑149 (MQ=11)
                                                            accacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTc                                                                                  <  6:179376/149‑1 (MQ=11)
                                                            accacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTc                                                                                  <  8:222566/149‑1 (MQ=11)
                                                            accacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTc                                                                                  >  6:13113/1‑149 (MQ=11)
                                                             ccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCg                                                                                 <  1:100552/149‑1 (MQ=11)
                                                             ccacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCg                                                                                 >  5:297476/1‑149 (MQ=11)
                                                              cacTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGt                                                                                <  6:67403/149‑1 (MQ=11)
                                                               acTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGtt                                                                               <  2:345323/149‑1 (MQ=11)
                                                                cTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtc                                                                              <  6:92040/149‑1 (MQ=11)
                                                                 ttCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttcttc                                                                                                       <  5:380988/123‑1 (MQ=255)
                                                                 ttCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTc                                                                                  >  2:166207/1‑144 (MQ=11)
                                                                 ttCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTtct                                                                             >  7:212138/1‑149 (MQ=11)
                                                                    aGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCaa                                                                          >  1:375593/1‑149 (MQ=37)
                                                                     gCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAAt                                                                         <  8:314132/149‑1 (MQ=11)
                                                                     gCGGAGTCTCTGGCATTCTTCAGTTACTGTGAGTCGGCCTGCGCACGCGTAAGTCGATCGCGTTCGTACTCAATAGTTCCTTGCTGTAGATCTGCCTAGCTGGCCTGCCGCAGTTCTACAACCTCCCGGCGAAGCTTTTCGTTAACAAt                                                                         >  1:49691/1‑149 (MQ=11)
                                                                       ggAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCt                                                                                                                    >  4:95684/1‑104 (MQ=11)
                                                                        gAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCACCCTCCCGGCGCAGCTTTTCGTTCTCAAttt                                                                       >  5:245885/1‑148 (MQ=11)
                                                                         aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAGCCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa                                                                     <  3:10439/149‑1 (MQ=11)
                                                                         aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa                                                                     <  2:185963/149‑1 (MQ=11)
                                                                         aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa                                                                     <  8:279921/149‑1 (MQ=11)
                                                                         aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa                                                                     <  7:65523/149‑1 (MQ=11)
                                                                         aGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa                                                                     <  5:70713/149‑1 (MQ=11)
                                                                           gctctTGCATTCTTCTGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCTATAGTTCCTGGCTGGAGTTCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGc                                                                   <  1:318274/148‑1 (MQ=11)
                                                                             tctGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTTAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAAc                                                                                                    <  7:120233/114‑1 (MQ=11)
                                                                               tGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATcc                                                               >  8:399623/1‑149 (MQ=11)
                                                                               tGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATcc                                                               >  5:201033/1‑149 (MQ=11)
                                                                               tGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATcc                                                               >  7:397901/1‑149 (MQ=11)
                                                                                ggCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATccc                                                              >  1:350108/1‑149 (MQ=11)
                                                                                 gCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCt                                                             >  1:264349/1‑149 (MQ=11)
                                                                                 gCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACCCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCaa                                                                          >  2:235917/1‑136 (MQ=11)
                                                                                 gCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACCCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCaa                                                                          <  1:235917/136‑1 (MQ=11)
                                                                                   aTTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCttt                                                           >  7:241477/1‑149 (MQ=11)
                                                                                    ttcttCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTTCTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTc                                                          >  5:120337/4‑149 (MQ=11)
                                                                                    ttcttCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATATGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTc                                                          >  3:269889/4‑149 (MQ=11)
                                                                                    ttcttCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGAGCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTc                                                          <  4:216474/146‑1 (MQ=11)
                                                                                    ttcttCAGTTACTGTGCGTCGGACTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTc                                                          >  2:327774/4‑149 (MQ=11)
                                                                                       ttCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttcttc                                                                                                       <  2:297990/101‑1 (MQ=255)
                                                                                       ttCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCgg                                                        >  6:334413/1‑148 (MQ=255)
                                                                                         cAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCgg                                                        <  5:174277/146‑1 (MQ=255)
                                                                                         cAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGc                                                       <  1:185143/147‑1 (MQ=255)
                                                                                         cAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGACGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCgg                                                        <  7:369760/146‑1 (MQ=11)
                                                                                          agttCCTGTGCGGCGGCCTGCGCACGCGTAAGTCGATGGCGTGCGTACTCAATAGTTCCTGGCTGGAGAGCTGCCTCGCTGGCCGGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATa                                                    <  4:276263/146‑1 (MQ=11)
                                                                                           gttCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATAc                                                   >  2:77600/3‑149 (MQ=12)
                                                                                            ttCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATAc                                                   >  7:318284/2‑148 (MQ=12)
                                                                                             tCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCa                                                 >  8:366031/1‑149 (MQ=12)
                                                                                              ccTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCAt                                                >  7:179425/1‑149 (MQ=12)
                                                                                                tgtgCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCa                                                                     <  7:296450/126‑1 (MQ=12)
                                                                                                tgtgCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCGGGCCTGCCGCAGTTCTTCAACCTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTc                                                          <  4:24552/137‑1 (MQ=255)
                                                                                                    cGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTTCCTGGCTGGAGATCTGCCTCGCTGGCCTGCCGCAGttc                                                                                                          <  4:317884/85‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTACTCAATAGTCCCTGGCTGGAGATCTGTCTCGCTGGCCTGCCGCAGTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCAC  >  NC_000913/1635982‑1636270

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