Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,370 C→G 7.9% intergenic (+303/+119) cspF → / ← quuQ Qin prophage; cold shock protein/Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,642,3700CG 7.9% 249.4 / 9.9 76intergenic (+303/+119)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base G (3/3);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

TTTTTATGAATAACCCCGTCTGTCTTGATGACTGGTTGATTGGCTTTAAAAGCTTATGCTGTACTTTGGCCGTAATAGCTCTGCTAATAATATAATAAGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCATACCCCTTATACCACTGATAATCGGTCAG  >  NC_000913/1642117‑1642602
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
tttttATGAATAACCCCGTCTGTCTTGATGACTGGTTGATTGGCTTTAAAAGCTTATGCTGTACTTTGGCCGTAATAGCTCTGCTAATAATATAATAAGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCAGCTCATCAGCGGGTTTTCTACtttt                                                                                                                                                                                                                      >  1:9280/1‑275 (MQ=255)
          aaaCCCCGGCTGCTGTGATGACTGGTTGATTGTCTTTAAAATCTTATGCTTAACTTTTGCCGACATAGCTCTGCTAAGTAAATAATAAGGCGACCCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGGTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACATTTTCTTAACGt                                                                                                                                                                                                             <  4:455219/273‑1 (MQ=255)
                cGTCTGTCTTGATGACTGGTTGATTGGCTTTAAAAGCTTATGCTGTACTTTGGCCGTAATAGCTCTGCTAATAATATAATAAGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTATTAAAGTCGGGta                                                                                                                                                                                                       >  1:143929/1‑274 (MQ=255)
                                                      tATGCTGTACTTTGGCCGTAATAGCTCTTCTAATAATATAATAAGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTGTAAACATCAGGCAAATAACTAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACAATACCATTGCCCCAGTTAAAGCACAAAACCCGCCCATCAGCGGGTTTTCTCCTTTTTCTTAACGTCGGGTATACAAAGCCCCTAGTTGAAAAAATGTTATCAATAttttt                                                                                                                                                                >  4:295258/1‑275 (MQ=255)
                                                                                             aataaGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACTTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCTTCGGGTCTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCACACATCATGTCGACCTCTTCAGCaaaa                                                                                                                          >  2:83577/1‑274 (MQ=255)
                                                                                             aataaGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCAATCTTCAGCaaaaa                                                                                                                         >  2:396915/1‑275 (MQ=255)
                                                                                                           ttGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcat                                                                                                                                                                                                                                         >  6:80438/1‑149 (MQ=255)
                                                                                                            tgtgTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGGCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCa                                                                                                          <  2:295262/275‑1 (MQ=255)
                                                                                                                tttAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                                                                                                                                                    <  7:199911/149‑1 (MQ=255)
                                                                                                                 ttAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCgg                                                                                                                                                                                                                                   <  5:85938/149‑1 (MQ=255)
                                                                                                                  tAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCggg                                                                                                                                                                                                                                  <  5:386645/149‑1 (MQ=255)
                                                                                                                  tAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTTAAAAATGCAAACATCATGTCGCCATATTAAGCAAAAATCATTTATCTCGTCACCTTc                                                                                                     >  4:88970/1‑274 (MQ=255)
                                                                                                                     ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtt                                                                                                                                                                                                                                <  6:309185/148‑1 (MQ=255)
                                                                                                                     ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTGAAAGCACAAAACCCGCTGATGAGCGGGtt                                                                                                                                                                                                                                >  5:309185/1‑148 (MQ=255)
                                                                                                                       aCATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttg                                                                                                                                                                                                                             <  6:166042/149‑2 (MQ=255)
                                                                                                                          ttGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                                                                                                                                                          <  6:137345/149‑1 (MQ=255)
                                                                                                                          ttGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAACAATTATGTCGCCATCTTCAGCAAAAATCATTTATCCCTTCACCTTCCTCAAtta                                                                                             >  4:30665/1‑273 (MQ=255)
                                                                                                                                tGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                                                                                                                    <  8:298365/149‑1 (MQ=255)
                                                                                                                                  gaagaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCcg                                                                                     >  4:260549/1‑274 (MQ=255)
                                                                                                                                       aaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                                                                                                                                                                 >  2:507024/1‑145 (MQ=255)
                                                                                                                                       aaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTc                                                                                                                                                                                                            >  4:423386/1‑150 (MQ=255)
                                                                                                                                          cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                                                                                                                                                          >  6:67031/1‑149 (MQ=255)
                                                                                                                                           aTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                                                                                                                                                        <  4:223012/118‑1 (MQ=255)
                                                                                                                                           aTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                                                                                                                                                             <  4:507019/145‑1 (MQ=255)
                                                                                                                                           aTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGt                                                                                                                                                                                                        <  2:423391/150‑1 (MQ=255)
                                                                                                                                            ttttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGt                                                                                                                                                                                                        >  7:339317/1‑149 (MQ=255)
                                                                                                                                            ttttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCTATTTTTTTTGAACATTGCAACCCTCATGTCGCCATCTTCAGCAAAATTCATTTATTTCGTAACCTTCCTCAATTGCGCTGCCGCGTAGCcttct                                                                         >  4:338917/1‑275 (MQ=255)
                                                                                                                                              ttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                                                                                                                                                      <  5:80438/149‑1 (MQ=255)
                                                                                                                                                    aTCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGCTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                                                                                                                                                                 >  4:545571/1‑132 (MQ=255)
                                                                                                                                                     tCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAACGCCCATCGTTGAAAAACTTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCTATCTTCAGCAAAATACATTTATTCTGTCCACTTACACAATTGCCATTCCGCTTCAGTTTCTTCCTGCCaa                                                                >  4:188985/1‑274 (MQ=255)
                                                                                                                                                      cAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGNATACAAAGCCCATCGTTTAAAAATTTTTCTCCATATTTTTTGAAAAATGCAAACATCATTTCCCCATCTTACACAAAAATNAATTATCTCCTAAACTCCCTCAATTNATCTACCCCGTCAGCTTCTTccacctacc                                                               >  3:94136/1‑268 (MQ=255)
                                                                                                                                                       aGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAACCATCATGTCGCCATCTTCAGCAAAAATCCTTTATCTCGTCACCTTCCTCAATTTCCCTTCCGCTTTTTCTTCTTCCTGCCatcc                                                              >  4:467474/1‑272 (MQ=255)
                                                                                                                                                        ggCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                                                                                                                                                             <  2:545576/132‑1 (MQ=255)
                                                                                                                                                         gCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCAt                                                                                                                                                                                           >  5:434299/1‑149 (MQ=255)
                                                                                                                                                          cAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATc                                                                                                                                                                                          >  8:14018/1‑149 (MQ=255)
                                                                                                                                                             aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                                                                                                                                                       <  7:14018/149‑1 (MQ=255)
                                                                                                                                                                aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTTGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTAc                                                     <  2:338921/275‑1 (MQ=255)
                                                                                                                                                                      aatagagacctataataATTGAACAGACCTGTGTCCCGTGATTCGATTAAAAAGTACCATTGCCCCATTTAAAGGACAATACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGAGTATACAAAGCCCATCGTTGACAAAATTTTATCCATATTGTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCATAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGtt                                                <  4:260131/263‑1 (MQ=255)
                                                                                                                                                                         ccAGCTAAATAATAATTTAACAGATATGAGTCCCGAGATGAGATTCAACATTACCATTGCCCCATTTAGAGCACGAAACCCGCTCTTCATCGGGTTTTCTACTTTTTCTTAAGGTCGGGTATGCAAAGCCCATCGTTGAAAAAATTTTa                                                                                                                                                                           >  8:295513/1‑149 (MQ=255)
                                                                                                                                                                              tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                                                                                                                                                               >  8:272102/1‑108 (MQ=255)
                                                                                                                                                                              tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                                                                                                                                                               <  7:272102/108‑1 (MQ=255)
                                                                                                                                                                                aataatAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTTACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACg                                     <  2:88971/275‑1 (MQ=255)
                                                                                                                                                                                       aggtAACAGACATGAGTCCCGGGAAGAGGTTCAACATTACCAATGCCCCATGTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATATGTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCAt                              <  4:321003/272‑1 (MQ=255)
                                                                                                                                                                                            acagacagGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCATAcc                           <  1:94136/273‑1 (MQ=255)
                                                                                                                                                                                            aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTAcc                                                    >  2:260135/1‑248 (MQ=255)
                                                                                                                                                                                            aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTcc                                  >  2:321007/1‑266 (MQ=255)
                                                                                                                                                                                               gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCATACCCCTta                      <  2:188987/275‑1 (MQ=255)
                                                                                                                                                                                                aCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTa                                                                                                                                                                           <  5:70497/126‑1 (MQ=255)
                                                                                                                                                                                                aCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTa                                                                                                                                                                           >  6:70497/1‑126 (MQ=255)
                                                                                                                                                                                                   tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAACAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTAACCTACCTCANATTCATTTTCGCGTATGTTTCTTCTTGCCAGCTCTTTGTTACAATTTTACCAATACCGCCCACATACACCTcattc                   >  3:119205/1‑272 (MQ=255)
                                                                                                                                                                                                       tCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCATAc                            <  1:119205/261‑1 (MQ=255)
                                                                                                                                                                                                         ccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGt                                                                                                                                           >  6:126710/1‑149 (MQ=255)
                                                                                                                                                                                                              aTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCa                                                                                                                                      <  6:434299/149‑1 (MQ=255)
                                                                                                                                                                                                                 agaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCCTCTTCAGCAAAAATCATTTATCTCGTCAACTTCCTCAATTGCGTCTCCGACNTGTCTTCTTCCTGCCAGCACTTTTTTTCCATTTTACCAATGACTCCAGCATCACCCTTCTAACACttaaaaacactt    >  3:124828/1‑264 (MQ=255)
                                                                                                                                                                                                                   aTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCTGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTTTCCATATTTTTTGAAAAATGCAACCATCATGTTGCCATCTTCAGCAATAATCATTTATCTCGTCACCTTCCTCACTTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCATACCTCTTATACCACTGATCATCGGTCAg  >  1:99480/1‑275 (MQ=255)
                                                                                                                                                                                                                   aTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTc                                                                                                                                 >  7:464028/1‑149 (MQ=255)
                                                                                                                                                                                                                    ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                                                                                                                    >  6:348563/1‑65 (MQ=255)
                                                                                                                                                                                                                    ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                                                                                                                    <  5:348563/65‑1 (MQ=255)
                                                                                                                                                                                                                    ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCa                                                                                                                                >  6:363058/1‑149 (MQ=255)
                                                                                                                                                                                                                          aTTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCaaaa                                                                                                                          >  5:308583/1‑149 (MQ=255)
                                                                                                                                                                                                                             aCCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTg                                                                    >  2:132078/1‑199 (MQ=255)
                                                                                                                                                                                                                               cATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAt                                                                                                                     <  6:308583/149‑1 (MQ=255)
                                                                                                                                                                                                                                 ttGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAg                                                                <  4:132077/199‑1 (MQ=255)
                                                                                                                                                                                                                                    ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAtt                                                                                                                    <  5:224928/145‑1 (MQ=255)
                                                                                                                                                                                                                                    ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAtt                                                                                                                    >  6:224928/1‑145 (MQ=255)
                                                                                                                                                                                                                                    ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                                                                                                                >  7:22130/1‑149 (MQ=255)
                                                                                                                                                                                                                                    ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                                                                                                                <  5:363058/149‑1 (MQ=255)
                                                                                                                                                                                                                                      ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtctc                                                                                                              <  5:67031/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAGAATGCAAAAATCATGTCGCCATCTTCAGCAAAAATCGTTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGGCGTCCGCATa                             <  2:295905/228‑1 (MQ=255)
                                                                                                                                                                                                                                         tttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTc                                                                                                           >  7:373515/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTc                                                                                                           <  8:339317/149‑1 (MQ=255)
                                                                                                                                                                                                                                                cacaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTcc                                                                                                    <  8:464028/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTc                                                                                                  <  5:126710/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                                                                                                                                                                     >  3:32285/1‑81 (MQ=255)
                                                                                                                                                                                                                                                     aaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAt                                                                                               >  6:226783/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                                                                                                                                                               >  7:259997/1‑84 (MQ=38)
                                                                                                                                                                                                                                                      aCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                                                                                                                                                               <  8:259997/84‑1 (MQ=37)
                                                                                                                                                                                                                                                       cccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                                                                                                                                                                 <  1:32285/81‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
TTTTTATGAATAACCCCGTCTGTCTTGATGACTGGTTGATTGGCTTTAAAAGCTTATGCTGTACTTTGGCCGTAATAGCTCTGCTAATAATATAATAAGCAGACTCATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCGCGTATGCTTCTTCCTGCCAGCACTTTGTTACCAGTTTACCAATGACGTCCGCATACCCCTTATACCACTGATAATCGGTCAG  >  NC_000913/1642117‑1642602

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