Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,370 C→G 6.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 6.9% 353.8 / 11.2 101intergenic (+303/+119)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (54/40);  new base G (3/4);  total (57/44)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 105 total aligned reads displayed.

CATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTT  >  NC_000913/1642222‑1642517
                                                                                                                                                    |                                                                                                                                                   
cATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtc                                                                                                                                                     <  1:648698/149‑1 (MQ=255)
    gtgtTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                                                                  <  4:152476/148‑1 (MQ=255)
       tttAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                                                              >  2:190629/1‑149 (MQ=255)
       tttAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                                                              <  1:539288/149‑1 (MQ=255)
        ttAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                                                              <  3:129568/148‑1 (MQ=255)
         tAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCggg                                                                                                                                            >  2:402115/1‑149 (MQ=255)
           gggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtt                                                                                                                                          <  1:692456/149‑1 (MQ=255)
           gggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtt                                                                                                                                          <  1:96381/149‑1 (MQ=255)
            ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                                                                         <  2:114316/149‑1 (MQ=255)
            ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtt                                                                                                                                          <  4:759377/148‑1 (MQ=255)
              aCATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                                                                       >  1:266393/1‑149 (MQ=255)
               cATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAAACAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                                                                       >  4:169680/1‑148 (MQ=255)
                   gTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                                                                  <  2:771807/149‑1 (MQ=255)
                    tACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                                                                  <  3:209091/148‑1 (MQ=255)
                         gaagaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                              >  3:39998/1‑147 (MQ=255)
                         gaagaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                              <  4:39998/147‑1 (MQ=255)
                         gaagaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCtt                                                                                                                            >  2:76775/1‑149 (MQ=255)
                         gaagaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                                                             >  4:262880/1‑148 (MQ=255)
                            gaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                                                                         >  2:593585/1‑149 (MQ=255)
                             aaaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                                                                         >  4:294860/1‑148 (MQ=255)
                               aaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                                                                       >  3:884840/1‑148 (MQ=255)
                               aaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTc                                                                                                                      >  1:606755/1‑149 (MQ=255)
                                 cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                                                                    <  1:76775/149‑1 (MQ=255)
                                 cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                                                                    <  1:994633/149‑1 (MQ=255)
                                 cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCg                                                                                                                     <  3:294860/148‑1 (MQ=255)
                                    tttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                                                                    >  3:24031/1‑146 (MQ=255)
                                    tttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                                                                    >  3:405581/1‑146 (MQ=255)
                                    tttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                                                                    <  4:405581/146‑1 (MQ=255)
                                     ttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                                                                >  2:388781/1‑149 (MQ=255)
                                         aCATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACa                                                                                                             >  3:967263/1‑148 (MQ=255)
                                           aTCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAg                                                                                                          >  1:723428/1‑149 (MQ=255)
                                           aTCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAg                                                                                                          >  1:535800/1‑149 (MQ=255)
                                            tCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAg                                                                                                          >  4:168522/1‑148 (MQ=255)
                                               ggCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTATTAACGTCGGGTACACAAAGccc                                                                                                       >  3:661901/1‑148 (MQ=255)
                                                gCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCa                                                                                                      >  4:185384/1‑148 (MQ=255)
                                                  aaaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCg                                                                                                   >  1:98855/1‑149 (MQ=255)
                                                      aaCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                                                                >  4:617097/1‑148 (MQ=255)
                                                      aaCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                                                                               >  2:160174/1‑149 (MQ=255)
                                                       aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaa                                                                                              >  1:646223/1‑149 (MQ=255)
                                                        ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGt                                                                                                                                           >  4:971051/1‑103 (MQ=255)
                                                        ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGt                                                                                                                                           <  3:971051/103‑1 (MQ=255)
                                                         cAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                                                                               <  4:405275/130‑1 (MQ=255)
                                                         cAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                                                                               >  3:405275/1‑130 (MQ=255)
                                                                ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                                                                      >  4:690949/1‑148 (MQ=255)
                                                                ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                                                                      >  3:658850/1‑148 (MQ=255)
                                                                  aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGTCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                                                                    <  3:559432/148‑1 (MQ=255)
                                                                    cTAAATAATAAGTTAACAGACATGAGTCCCGGGATGTGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCa                                                                                 >  2:1077414/1‑149 (MQ=255)
                                                                      aaaTAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt                                                                                         >  4:49432/1‑139 (MQ=255)
                                                                      aaaTAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt                                                                                         <  3:49432/139‑1 (MQ=255)
                                                                          aatnaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCANTGCCCCATTTAAAGCANNAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAGGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                                                                            >  3:680548/1‑148 (MQ=255)
                                                                          aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                                                                           >  2:1049924/1‑149 (MQ=255)
                                                                            taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                                                                         <  1:190629/149‑1 (MQ=255)
                                                                            taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATGTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                                                                          <  4:647606/148‑1 (MQ=255)
                                                                                  aaCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                                                                   >  2:672613/1‑149 (MQ=255)
                                                                                   aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                                                                  >  3:188589/1‑85 (MQ=255)
                                                                                   aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                                                                  <  4:188589/85‑1 (MQ=255)
                                                                                   aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                                                                   <  4:884840/148‑1 (MQ=255)
                                                                                   aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                                                                  <  1:28026/149‑1 (MQ=255)
                                                                                    cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAAGCCGCTCATCAGCGGGTTTTCTACTTTTCCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                                                                  >  4:1057333/1‑148 (MQ=255)
                                                                                    cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                                                                 <  1:672613/149‑1 (MQ=255)
                                                                                     aGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                                                                   <  4:695841/146‑1 (MQ=255)
                                                                                     aGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                                                                   >  3:695841/1‑146 (MQ=255)
                                                                                        cATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                                                                           <  4:162294/135‑1 (MQ=255)
                                                                                        cATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                                                                           >  3:162294/1‑135 (MQ=255)
                                                                                        cATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa                                                              <  3:971791/148‑1 (MQ=255)
                                                                                          tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca                                                           >  2:1009099/1‑149 (MQ=255)
                                                                                            aGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatc                                                         >  2:85007/1‑149 (MQ=255)
                                                                                                ccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg                                                      <  4:1029506/148‑1 (MQ=255)
                                                                                                 cGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                                                                           >  4:404425/1‑78 (MQ=255)
                                                                                                 cGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                                                                           <  3:404425/78‑1 (MQ=255)
                                                                                                       gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAt                                               >  3:171825/1‑148 (MQ=255)
                                                                                                            tCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCa                                          >  4:242440/1‑148 (MQ=255)
                                                                                                            tCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAg                                         >  1:699677/1‑149 (MQ=255)
                                                                                                            tCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAg                                         <  2:699677/149‑1 (MQ=255)
                                                                                                               aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCa                                       <  3:185384/148‑1 (MQ=255)
                                                                                                                cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCaaa                                     >  1:844673/1‑149 (MQ=255)
                                                                                                                    aCCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAAt                                  <  3:381655/148‑1 (MQ=255)
                                                                                                                       aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAtt                              <  2:606755/149‑1 (MQ=255)
                                                                                                                       aTTGCCCCATATAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTCTTTGAAAAATGCAAACATCATGCCGCCATCTTCAGCAAAAATCAt                               >  3:599684/1‑148 (MQ=255)
                                                                                                                          gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAt                           <  1:497417/149‑1 (MQ=255)
                                                                                                                           ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                          >  1:954357/1‑149 (MQ=255)
                                                                                                                           ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAt                           <  4:941480/148‑1 (MQ=255)
                                                                                                                             ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtctc                        <  1:1053612/149‑1 (MQ=255)
                                                                                                                             ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtctc                        >  1:510113/1‑149 (MQ=255)
                                                                                                                               aTTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTGTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGt                      >  1:186825/1‑149 (MQ=255)
                                                                                                                                tttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGt                      >  3:67629/1‑148 (MQ=255)
                                                                                                                                tttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTc                     >  1:361914/1‑149 (MQ=255)
                                                                                                                                    aaGCACAAGACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATGTTCAGCAAAAATCATTTATCTCGTCAct                  >  4:910467/1‑147 (MQ=255)
                                                                                                                                    aaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCAcc                  >  3:262727/1‑148 (MQ=255)
                                                                                                                                     agtagaaaaCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCt                 >  4:773431/6‑148 (MQ=39)
                                                                                                                                     aGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGACGCCATCTTCAGCAAAAATCATTTATCTCGTCACCt                 >  3:705859/1‑148 (MQ=255)
                                                                                                                                        agaaaaCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTcc              >  4:220798/3‑148 (MQ=255)
                                                                                                                                        acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCt             <  1:1077414/149‑1 (MQ=255)
                                                                                                                                           aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCaa          <  1:388781/149‑1 (MQ=255)
                                                                                                                                           aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCa           <  3:19396/148‑1 (MQ=255)
                                                                                                                                            aaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCaa          <  4:67629/148‑1 (MQ=255)
                                                                                                                                            aaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAt         <  2:535800/149‑1 (MQ=255)
                                                                                                                                             aCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAtt        <  2:954357/149‑1 (MQ=255)
                                                                                                                                                    catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCGGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCtt  <  3:690949/148‑1 (MQ=255)
                                                                                                                                                    catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTgcgcat  <  4:658850/148‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTT  >  NC_000913/1642222‑1642517

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