Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,086,965 G→T 5.0% intergenic (+116/‑37) fdhD → / → yiiG formate dehydrogenase formation protein/DUF3829 family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,086,9650GT 5.0% 2568.4 / 15.4 775intergenic (+116/‑37)fdhD/yiiGformate dehydrogenase formation protein/DUF3829 family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (243/465);  new base T (65/2);  total (308/467)
Fisher's exact test for biased strand distribution p-value = 9.45e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 784 total aligned reads displayed.

TTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCA  >  NC_000913/4086818‑4087114
                                                                                                                                                   |                                                                                                                                                     
tttATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGtt                                                                                                                                                     <  2:8033873/150‑1 (MQ=255)
  tataCCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCa                                                                                                                                                   <  2:11011963/150‑1 (MQ=255)
  tataCCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCa                                                                                                                                                   <  2:10538377/150‑1 (MQ=255)
  tataCCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCGTATTTCATTTATTATGTTcc                                                                                                                                                   >  2:77089/1‑149 (MQ=255)
  tataCCCATCCGCAGAGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCa                                                                                                                                                   <  2:15111585/150‑1 (MQ=255)
   ataCCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCAc                                                                                                                                                  <  1:8836515/150‑1 (MQ=255)
       ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATTTTCActcta                                                                                                                                              >  1:4743730/1‑149 (MQ=255)
       ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCActctc                                                                                                                                              >  2:10432789/1‑150 (MQ=255)
       ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATTTTCActgaa                                                                                                                                              >  1:7911353/1‑147 (MQ=255)
        cATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCActctct                                                                                                                                             <  2:12916764/150‑1 (MQ=255)
         aTCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCATGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCtt                                                                                                                                            <  2:14985219/150‑1 (MQ=255)
         aTCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCtt                                                                                                                                            <  2:3190281/150‑1 (MQ=255)
           ccGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTAt                                                                                                                                          <  2:10249119/150‑1 (MQ=255)
            cGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATc                                                                                                                                         <  2:13515287/150‑1 (MQ=255)
            cGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATc                                                                                                                                         <  1:9237914/150‑1 (MQ=255)
            cGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATGTTCACTCTCTTATc                                                                                                                                         >  2:4267929/1‑150 (MQ=255)
              cAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTAACTCTCTtatata                                                                                                                                       >  1:3608703/1‑150 (MQ=255)
              cAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTGTTATTTTTTAATTTTTTATGTTAACTATATTATCta                                                                                                                                       >  2:11751488/1‑150 (MQ=255)
                gCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCtata                                                                                                                                     <  2:2824328/150‑1 (MQ=255)
                      gAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTta                                                                                                                               >  1:1201177/1‑150 (MQ=255)
                       aGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTAATTTATTATGTTCACTCTCTGATCTATATATTtat                                                                                                                              >  2:7532105/1‑150 (MQ=255)
                        gCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTAATTTATTATGTTGACTCTCTTATATAAATATTtata                                                                                                                             >  1:11118663/1‑150 (MQ=255)
                         cAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATaa                                                                                                                            <  2:349033/150‑1 (MQ=255)
                          aattaattAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTGTTTTTTCATTTATTATGTTCACTCTGTTATATATATATTTATAAg                                                                                                                           >  1:3525508/1‑150 (MQ=255)
                           attaattaAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTAATTTTTTATGTTCACTCTCGTATCTATATATTTATAAgg                                                                                                                          >  1:12989119/1‑150 (MQ=255)
                              aattaaAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCaa                                                                                                                       <  2:5380894/150‑1 (MQ=255)
                               attaaAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAAt                                                                                                                      >  2:5559043/1‑150 (MQ=255)
                               attaaAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAAt                                                                                                                      >  2:2613290/1‑150 (MQ=255)
                                ttaaAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAAtt                                                                                                                     >  2:10789888/1‑150 (MQ=255)
                                  aaaaTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTAACTCTCTTATCTATATTTTTATAAGGGAATTaa                                                                                                                   >  1:7270288/1‑150 (MQ=255)
                                   aaaTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTaaa                                                                                                                  <  1:10471677/150‑1 (MQ=255)
                                   aaaTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTaaa                                                                                                                  >  2:6005888/1‑150 (MQ=255)
                                        aGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGaaa                                                                                                             <  2:8435077/150‑1 (MQ=255)
                                        aGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGaaa                                                                                                             <  2:8812219/150‑1 (MQ=255)
                                        aGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGaaa                                                                                                             <  2:9304670/150‑1 (MQ=255)
                                        aGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGaaa                                                                                                             <  1:15283899/150‑1 (MQ=255)
                                        aGTCCTGCGAAATATATTAAATTGCAGGACTTATCCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGaaa                                                                                                             <  2:1431393/150‑1 (MQ=255)
                                         gTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGaaaa                                                                                                            <  1:8201039/150‑1 (MQ=255)
                                         gTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATTTTCACTCTCATATATATATATTTTTAATGCCATTAAATGaaaa                                                                                                            >  2:4007375/1‑150 (MQ=255)
                                           ccTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAgg                                                                                                          <  2:4886811/150‑1 (MQ=255)
                                              gCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAAt                                                                                                       <  1:3711587/150‑1 (MQ=255)
                                                   atatatTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATttatt                                                                                                  <  1:2915710/150‑1 (MQ=255)
                                                      tatTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATc                                                                                               <  1:13247699/150‑1 (MQ=255)
                                                        ttAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAATGCAATTAAATGAAAAGGAATTTATTATCtt                                                                                             >  2:9051940/1‑150 (MQ=255)
                                                        ttAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATGTTAACTATCTTATATATATATTTATAAGGCAAATAAAATAAAAGGAAATTATTATAtt                                                                                             >  1:11427937/1‑150 (MQ=255)
                                                         tAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATTTTCACTATATTATCTATATATTTATAAGGCAATTAAATTAAAAGGAATTtattataata                                                                                            >  2:13792976/1‑146 (MQ=255)
                                                         tAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTc                                                                                            <  1:4330487/150‑1 (MQ=255)
                                                         tAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTGTTATTTTTTTATTTATTTTGTTGACTCTCTTTTATATATATTTATAAGGGAATTAAATGAAAATTAATttattttatta                                                                                            >  1:5724872/1‑144 (MQ=255)
                                                         tAAAGTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTc                                                                                            <  1:2258113/150‑1 (MQ=255)
                                                          aaaTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTcc                                                                                           <  1:1077148/150‑1 (MQ=255)
                                                            aTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGc                                                                                         <  2:2012160/150‑1 (MQ=255)
                                                            aTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGc                                                                                         <  1:1884617/150‑1 (MQ=255)
                                                                        aTTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGc                                                                             <  1:6948161/150‑1 (MQ=255)
                                                                        aTTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGc                                                                             <  1:13961971/150‑1 (MQ=255)
                                                                        aTTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTTTTTTAAATCTATTATCTATATATTTATAAGGCAAATAAATTAAAAGGAAATTATTATCTTAAGCGATtatattatc                                                                             >  1:8628989/1‑145 (MQ=255)
                                                                          tCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATTTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTACGCTATTATAGTCGCca                                                                           >  2:8656142/1‑150 (MQ=255)
                                                                              ttCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATGATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCcatcat                                                                       <  2:12019971/150‑1 (MQ=255)
                                                                                 gTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTc                                                                    <  1:10026286/150‑1 (MQ=255)
                                                                                 gTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATTTTCACACTCTTATCTATATATTTATAAGGCAATTAAATTAAAAGGAATTTATTATATTACGCTATTATAGTAGCCATCATGta                                                                    >  2:7266723/1‑149 (MQ=255)
                                                                                        ttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGt                                                             <  2:1391928/150‑1 (MQ=255)
                                                                                              ttttatttATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACg                                                       <  1:10181242/150‑1 (MQ=255)
                                                                                                 tatttatAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGt                                                    <  1:7406186/150‑1 (MQ=255)
                                                                                                   tttatAAACAACATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTtg                                                  <  2:14508338/150‑1 (MQ=255)
                                                                                                     tataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTtgtg                                                <  1:15543541/150‑1 (MQ=255)
                                                                                                     tataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTtgtg                                                <  1:2551169/150‑1 (MQ=255)
                                                                                                     tataAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCGATTATAGTAGCCATCATGTCCCTCGGTCTGACGGGGtgtg                                                >  2:13281535/1‑150 (MQ=255)
                                                                                                      ataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtga                                               <  2:12726732/150‑1 (MQ=255)
                                                                                                      ataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtga                                               <  2:4445382/150‑1 (MQ=255)
                                                                                                       taAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgat                                              <  2:1697499/150‑1 (MQ=255)
                                                                                                        aaaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgatg                                             <  2:9562511/150‑1 (MQ=255)
                                                                                                         aaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgatga                                            <  1:9950210/150‑1 (MQ=255)
                                                                                                         aaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtga                                               >  2:6246946/1‑147 (MQ=255)
                                                                                                            aaaaTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaa                                         <  2:12322683/150‑1 (MQ=255)
                                                                                                            aaaaTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaa                                         <  1:9521091/150‑1 (MQ=255)
                                                                                                            aaaaTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaa                                         <  2:14831743/150‑1 (MQ=255)
                                                                                                               aTGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaa                                      <  1:5057808/150‑1 (MQ=255)
                                                                                                                tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  2:3669907/150‑1 (MQ=255)
                                                                                                                tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  1:605617/150‑1 (MQ=255)
                                                                                                                tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  2:4604352/150‑1 (MQ=255)
                                                                                                                tGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  1:7715552/150‑1 (MQ=255)
                                                                                                                  cGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGc                                   <  1:11890634/150‑1 (MQ=255)
                                                                                                                      tAGTGGCGTTTTCTTTTTTAATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAAATAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGaa                               >  2:5164347/1‑150 (MQ=255)
                                                                                                                           gCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGa                          >  2:2255657/1‑150 (MQ=255)
                                                                                                                            cGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGaa                         <  1:12259374/150‑1 (MQ=255)
                                                                                                                             gTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGaaa                        <  2:11347982/150‑1 (MQ=255)
                                                                                                                              ttttCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAAc                       <  1:11628630/150‑1 (MQ=255)
                                                                                                                                   ttttttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTcc                  <  2:11489517/150‑1 (MQ=255)
                                                                                                                                    tttttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATATTACGCTATTATAGTCGCCATCATGTCCCACGGTCTGACGGGTTGTGATGATAAAAAATCCGAAACAGAAACGCTccc                 >  1:4010377/1‑150 (MQ=255)
                                                                                                                                      tttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGc               >  2:15517500/1‑150 (MQ=255)
                                                                                                                                       ttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTACGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGACGAAACAGAAACGCTCACGcc              >  2:11243272/1‑150 (MQ=255)
                                                                                                                                        tcttttATTATGTACACTCTCTTATATATATATTTATACGGAAATTAAATGAAAATTATTTTATTATCTTCCGCTATTATAGTCGCCATCTTGTCCCTCGGTCTGAGGGGTTGTGATTATAAAAAAGCCTAAACAGAAACGCTCCCGCCt             <  1:12573074/147‑1 (MQ=255)
                                                                                                                                          aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           >  2:5883148/1‑150 (MQ=255)
                                                                                                                                          aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           <  2:6871725/150‑1 (MQ=255)
                                                                                                                                          aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           >  1:6582266/1‑150 (MQ=255)
                                                                                                                                          aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           >  1:14190657/1‑150 (MQ=255)
                                                                                                                                          aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           >  1:12874595/1‑150 (MQ=255)
                                                                                                                                           tttATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGcc          >  2:6588468/1‑150 (MQ=255)
                                                                                                                                             tattatTTTCACTATCTTATCTATATATTTATAAGGCAATTAAATTAAAAGTAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGTTGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCaa        <  1:311222/150‑1 (MQ=255)
                                                                                                                                              attatGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAAt       >  1:14724707/1‑150 (MQ=255)
                                                                                                                                               ttatGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATa      >  1:14957147/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCA  >  NC_000913/4086818‑4087114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: