Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,635,974 T→G 100% pseudogene (394/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,635,9740TG99.4% 1693.5 / ‑1.4 883pseudogene (394/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 (2/5);  new base G (335/541);  total (337/546)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 883 total aligned reads displayed.

GCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTAC  >  NC_000913/1635832‑1636119
                                                                                                                                              |                                                                                                                                                 
gcgcGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTc                                                                                                                                             <  1:223135/149‑1 (MQ=35)
                    cacccccGGTATCAGTTCATCCAGCGCGGCGGCTTTGTTCATGGCGTTGATGATATCCCGTTTCAGGAAATCACCATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCACTTTCAc                                                                                                                         <  6:62314/146‑1 (MQ=1)
                                                        gtgcATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCg                                                                                      <  1:370292/145‑1 (MQ=1)
                                                                                                         tCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGGGTCTCTGGCATTCTTCAGTTCCTGTGCGt                                     <  4:165347/148‑1 (MQ=9)
                                                                                                         tCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGt                                     <  4:251481/148‑1 (MQ=9)
                                                                                                          ccAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTc                                    <  8:54106/148‑1 (MQ=12)
                                                                                                          ccAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTc                                    <  2:232339/148‑1 (MQ=12)
                                                                                                          ccAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTc                                    <  5:337823/148‑1 (MQ=12)
                                                                                                          ccAGTTCCGGAAAACGCCGCTGCACCGACAGGGGCAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTc                                    <  1:80632/148‑1 (MQ=255)
                                                                                                           cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGTAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGGTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCg                                   <  1:268007/148‑1 (MQ=11)
                                                                                                           cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCCCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCgg                                  >  5:147/1‑149 (MQ=255)
                                                                                                           cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCgg                                  >  6:27571/1‑149 (MQ=12)
                                                                                                           cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCg                                   <  8:65726/148‑1 (MQ=12)
                                                                                                           cAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCg                                   <  7:48777/148‑1 (MQ=12)
                                                                                                            aGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGc                                 >  5:35899/1‑149 (MQ=12)
                                                                                                            aGTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGc                                 >  2:144754/1‑149 (MQ=12)
                                                                                                             gTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGcc                                >  2:166170/1‑149 (MQ=12)
                                                                                                             gTTCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGcc                                >  8:242203/1‑149 (MQ=12)
                                                                                                              ttCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCGGTCGAGAATACTGGCAATTTCACCGGCGAGCCGCGACAGCACGAAAGTACAGACTGCGGTTTCCACCACGTCAGCGGCGTCTCTGGCAGTCTTCGGGTCCGGTGCGTCGGCCt                               >  4:237432/1‑149 (MQ=11)
                                                                                                              ttCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCt                               >  5:81410/1‑149 (MQ=12)
                                                                                                              ttCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCt                               <  6:350490/149‑1 (MQ=12)
                                                                                                              ttCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCt                               <  5:301581/149‑1 (MQ=12)
                                                                                                              ttCCGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACAGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCt                               >  1:137170/1‑149 (MQ=255)
                                                                                                                ccGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGGTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgc                             >  4:66282/1‑149 (MQ=255)
                                                                                                                 cGGAAAACGCCGCTGCACCGTCAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGACAGTACAGAATGCGGTTTCCACCACTTCAGCAGAGGCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                            <  2:375449/149‑1 (MQ=11)
                                                                                                                 cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                            >  8:289865/1‑149 (MQ=12)
                                                                                                                 cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                            >  4:310188/1‑149 (MQ=12)
                                                                                                                 cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                            <  7:373090/149‑1 (MQ=12)
                                                                                                                 cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgcg                            >  8:110283/1‑149 (MQ=12)
                                                                                                                 cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgc                             <  2:108484/148‑1 (MQ=12)
                                                                                                                 cGGAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgc                             <  5:63586/148‑1 (MQ=12)
                                                                                                                  ggAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCg                                   >  4:189128/1‑141 (MQ=12)
                                                                                                                   gAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCa                          >  5:38062/1‑149 (MQ=12)
                                                                                                                   gAAAACGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCa                          >  3:75408/1‑149 (MQ=12)
                                                                                                                    aaaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgc                             <  8:358962/145‑1 (MQ=12)
                                                                                                                    aaaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAc                         <  7:265652/149‑1 (MQ=12)
                                                                                                                    aaaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAc                         <  3:12218/149‑1 (MQ=12)
                                                                                                                    aaaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTCCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAc                         <  5:209803/149‑1 (MQ=255)
                                                                                                                     aaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcg                        >  7:64713/1‑149 (MQ=12)
                                                                                                                     aaaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcg                        >  7:330530/1‑149 (MQ=12)
                                                                                                                      aaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGACTACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                       <  3:188016/149‑1 (MQ=255)
                                                                                                                      aaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGACTGCGGGTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                       <  3:33992/149‑1 (MQ=11)
                                                                                                                      aaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                       <  1:304235/149‑1 (MQ=12)
                                                                                                                      aaCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAAAGCGGTTTCCACCCCTTCAGCGGAGTCTCTGGCATGCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                       >  4:334926/1‑149 (MQ=11)
                                                                                                                       aCGCCGCTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgcg                      <  3:169109/149‑1 (MQ=12)
                                                                                                                        cgccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGt                     <  2:340006/149‑1 (MQ=12)
                                                                                                                         gccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCCCCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGt                     >  3:150530/1‑148 (MQ=255)
                                                                                                                         gccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTgc                             >  1:5504/1‑140 (MQ=12)
                                                                                                                         gccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGt                     >  1:248836/1‑148 (MQ=12)
                                                                                                                          ccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTaa                   <  6:204041/149‑1 (MQ=12)
                                                                                                                          ccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTaa                   <  8:134359/149‑1 (MQ=12)
                                                                                                                          ccgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTaa                   <  1:149439/149‑1 (MQ=12)
                                                                                                                           cgcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGt                     >  4:275584/1‑146 (MQ=12)
                                                                                                                            gcTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCCATCAGCGCCAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTGTGGCATTCTTCAGATCCTGTGCGTCGGCCTGCGCACGCGAAAGt                 >  2:197091/1‑149 (MQ=11)
                                                                                                                             cTGCACCGACAGGGGGAGCCCGTCGAGAATACTTGAAATTTCACCTGCGATCCGCGACAGCACGAACGTAAAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTTGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTc                <  3:298790/149‑1 (MQ=11)
                                                                                                                             cTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTc                <  3:38519/149‑1 (MQ=12)
                                                                                                                             cTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTc                <  4:145855/149‑1 (MQ=12)
                                                                                                                             cTGCACCGACAGGGGGAGCCCGTCGAGAATACTGGCAAGTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTc                <  8:379888/149‑1 (MQ=255)
                                                                                                                             cTGCACCGACAGGGGGAGCCCGGCGAGACTACTGGCAATTTCACCTGCGATCCGCGACCGCCCGAAAGTACAGAAGGCGGGGTCCACCACTTCAGCGGAGTCGCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTc                <  7:178361/149‑1 (MQ=11)
                                                                                                                                cACCGACAGGGGGTGCGCGTCGAGAATACTGCCAATTTCACCTGCGATCCGATACAGCACGATAGTACAGAATGCGGTTTCCACCAATTAAGAGGAGTCTCTGGCTTTGTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGa              <  2:290444/148‑1 (MQ=11)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCg                                   <  3:264243/127‑1 (MQ=12)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGa              <  3:240168/148‑1 (MQ=12)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGAt             <  5:190058/149‑1 (MQ=12)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGAt             >  2:42885/1‑149 (MQ=12)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGAt             >  5:101560/1‑149 (MQ=12)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGAt             <  8:253064/149‑1 (MQ=12)
                                                                                                                                cACCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGAt             <  5:260256/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCGCGTCGAGAATACTGGCAATTTCACCTTCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTTGCATTCTTCAGTTCCTGTTCGTCGGCCTGTGCACGCGTAAGTCGATg            <  1:344977/149‑1 (MQ=11)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGCTCCGCGCCAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            >  3:216474/1‑149 (MQ=11)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  3:65177/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            >  6:16126/1‑149 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  4:367715/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  6:2314/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  1:397465/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  7:305912/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  5:16126/149‑1 (MQ=12)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTGCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  1:105996/149‑1 (MQ=255)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAACTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  1:193928/149‑1 (MQ=255)
                                                                                                                                 aCCGACAGGGGGAGCCCGTCGAGAATACGGGCAATTTCACCTGCGATCCGCGAGAGGACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCTTTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            <  7:339108/149‑1 (MQ=11)
                                                                                                                                  ccGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATgg           <  4:12634/149‑1 (MQ=12)
                                                                                                                                  ccGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGGCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATgg           <  7:143346/149‑1 (MQ=255)
                                                                                                                                   cGACAGGGGGATCCCGTCTAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCAACCACTTCAGCGGAGTCTCTGGCATTCTTGAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc          <  6:263513/149‑1 (MQ=11)
                                                                                                                                   cGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCCCCCCTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc          >  6:118426/1‑149 (MQ=11)
                                                                                                                                   cGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc          <  7:118279/149‑1 (MQ=12)
                                                                                                                                   cGACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc          <  6:82914/149‑1 (MQ=12)
                                                                                                                                    gACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCg         >  2:60313/1‑149 (MQ=12)
                                                                                                                                    gACAGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCg         <  2:254492/149‑1 (MQ=12)
                                                                                                                                       aGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTc      <  8:309480/149‑1 (MQ=12)
                                                                                                                                       aGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTc      <  2:86719/149‑1 (MQ=12)
                                                                                                                                       aGGGGGAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTc      <  7:230424/149‑1 (MQ=12)
                                                                                                                                        gggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCg     <  5:129796/149‑1 (MQ=12)
                                                                                                                                        gggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCg     <  3:280115/149‑1 (MQ=12)
                                                                                                                                         ggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGCCAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCAcgc                       >  6:36712/1‑130 (MQ=255)
                                                                                                                                         ggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCGCGCGTAAGTCGATGGCGTTCg     >  6:33281/1‑148 (MQ=255)
                                                                                                                                         ggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCCCGCGTAAGTCGATGGCGTTCg     >  8:179074/1‑148 (MQ=255)
                                                                                                                                         ggggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGt    <  6:189039/149‑1 (MQ=12)
                                                                                                                                         ggggAGCCCATCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGt        <  7:288703/145‑1 (MQ=255)
                                                                                                                                          gggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATg            >  6:305372/1‑140 (MQ=12)
                                                                                                                                           ggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGc          >  6:139945/1‑141 (MQ=12)
                                                                                                                                           ggAGCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTCTTCAGTTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTAc  >  8:22115/1‑149 (MQ=12)
                                                                                                                                              |                                                                                                                                                 
GCGCGGTTATATTCACTCAGCAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTGCGCACGCGTAAGTCGATGGCGTTCGTAC  >  NC_000913/1635832‑1636119

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: