Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,683,472 G→T 9.1% intergenic (‑25/+158) dctA ← / ← yhjK C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,4720GT 9.1% 453.4 / 25.7 151intergenic (‑25/+158)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (70/67);  new base T (7/7);  total (77/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

Only 100 of 154 total aligned reads displayed.

CGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAG  >  NC_000913/3683324‑3683619
                                                                                                                                                    |                                                                                                                                                   
cGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGGCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAg                                                                                                                                                     <  7:208359/149‑1 (MQ=255)
     ccGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCtt                                                                                                                                                <  7:143300/149‑1 (MQ=255)
      cGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTg                                                                                                                                               <  3:232042/149‑1 (MQ=255)
      cGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTg                                                                                                                                               <  2:408965/149‑1 (MQ=255)
       gTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTg                                                                                                                                               >  6:232190/1‑148 (MQ=255)
        tCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGAACTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGt                                                                                                                                              >  4:20308/1‑148 (MQ=255)
          gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAg                                                                                                                                           >  7:376487/1‑149 (MQ=255)
               gCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGGGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTAt                                                                                                                                      <  5:232190/149‑1 (MQ=255)
                  gTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggg                                                                                                                                   >  8:233135/1‑149 (MQ=255)
                   tttCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATgggg                                                                                                                                  >  8:174835/1‑149 (MQ=255)
                     tCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggg                                                                                                                                   >  7:344017/1‑146 (MQ=255)
                     tCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATggg                                                                                                                                   <  8:344017/146‑1 (MQ=255)
                       aTTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAcc                                                                                                                              >  7:361087/1‑149 (MQ=255)
                        tttGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTg                                                                                                                            >  7:129417/1‑149 (MQ=255)
                             tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAATGGCTTTTAAACAGAGAGGTTTTCATTGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGc                                                                                                                        >  3:317177/1‑149 (MQ=255)
                             tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTAATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCgt                                                                                                                         >  3:275/1‑147 (MQ=255)
                               gCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGcc                                                                                                                       <  1:403571/148‑1 (MQ=255)
                               gCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCt                                                                                                                      >  8:297951/1‑149 (MQ=255)
                               gCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCt                                                                                                                      <  1:166739/149‑1 (MQ=255)
                                 cTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATg                                                                                                                                     >  1:119573/1‑132 (MQ=255)
                                    tttcttgATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGATAGTAAAGGCTTTTAACCAGAGAGGTTTTCATGGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaaaa                                                                                                                 <  2:158043/143‑1 (MQ=255)
                                    tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTaa                                                                                                                   <  3:156912/147‑1 (MQ=255)
                                    tttCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTGTTACACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCTGCCTTaaaa                                                                                                                 >  6:189349/1‑149 (MQ=255)
                                       cAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTAc                                                                                                                               >  1:319558/1‑131 (MQ=255)
                                       cAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTAc                                                                                                                               <  2:319558/131‑1 (MQ=255)
                                         ggATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAaca                                                                                                            <  7:123256/149‑1 (MQ=255)
                                           aTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                          <  8:102432/149‑1 (MQ=255)
                                           aTAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                          <  3:34384/149‑1 (MQ=255)
                                              gAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCa                                                                                                        <  3:350599/148‑1 (MQ=255)
                                               aaaTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGa                                                                                                      >  6:143064/1‑148 (MQ=255)
                                               aaaTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTCTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAACTGCGGCCTTAAAATAACACCCAGa                                                                                                      >  6:327305/1‑149 (MQ=255)
                                                    gCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGTCCTTAAAATAACACCCAGACAACAt                                                                                                <  8:129417/149‑1 (MQ=255)
                                                      cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGGGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATc                                                                                               <  8:268704/149‑1 (MQ=255)
                                                         ggAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca                                                                                            <  2:388291/149‑1 (MQ=255)
                                                           agaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGa                                                                                          <  3:25192/149‑1 (MQ=255)
                                                             aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaaa                                                                                        <  8:185669/149‑1 (MQ=255)
                                                             aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaaa                                                                                        <  8:361087/149‑1 (MQ=255)
                                                             aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaaa                                                                                        <  2:146874/149‑1 (MQ=255)
                                                               tACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCATAGGTCTTGTAAGGTTATggg                                                                                                                                   >  4:333329/1‑104 (MQ=255)
                                                                   aaTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTaaaa                                                                                                                 <  5:215647/117‑1 (MQ=255)
                                                                       gCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGa                                                                              >  6:347514/1‑149 (MQ=255)
                                                                       gCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTAAGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGa                                                                              <  1:38927/148‑1 (MQ=255)
                                                                          aTCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCAGAGGGGGTCCTTTAGTAAAACCACAGGTCTGGTAAGGTTATGGGGTACCTGCGGCCTATAATTAACACCCAGACAACATCACAGAAATGTACCTGGATCa                                                                           >  4:269221/1‑149 (MQ=255)
                                                                               tGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                      >  7:109019/1‑149 (MQ=255)
                                                                                   aGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGaa                                                                   <  7:213482/148‑1 (MQ=255)
                                                                                    ggACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGaaa                                                                  >  7:426599/1‑148 (MQ=255)
                                                                                         tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGTTCCTTCGGCCTTAAAATAAAACCCCGACAAAATAACAGAAATTTTCCTGGATCATAATTTAAACAAAAAg                                                           >  5:400717/1‑149 (MQ=255)
                                                                                           aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAg                                                           <  3:69793/148‑1 (MQ=255)
                                                                                             aGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGtt                                                        <  7:6415/149‑1 (MQ=255)
                                                                                              gTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTa                                                       >  1:121830/1‑148 (MQ=255)
                                                                                                aaaGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTAACTGCGGACTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTaa                                                      >  3:119358/1‑148 (MQ=255)
                                                                                                  aGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAg                                                  <  1:371065/149‑1 (MQ=255)
                                                                                                    gCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTCCCTGGATCATAATTGAAACAAAAAGGGTAAAAAGt                                                 >  2:297269/1‑149 (MQ=255)
                                                                                                     cTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGtt                                                <  8:109019/149‑1 (MQ=255)
                                                                                                      ttttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               >  7:142439/1‑148 (MQ=255)
                                                                                                      atttAAACAGATAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACATGGATCAGAATAGAAACAAAAAGGTTAAAAAGatt                                               >  4:290136/2‑148 (MQ=255)
                                                                                                       tttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               >  5:89464/1‑148 (MQ=255)
                                                                                                        ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               <  3:196904/147‑1 (MQ=255)
                                                                                                        ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               <  7:2552/147‑1 (MQ=255)
                                                                                                        ttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGCAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAt                                                                       <  4:51271/123‑1 (MQ=255)
                                                                                                         tAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAg                                            <  6:183566/148‑1 (MQ=255)
                                                                                                          aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                              <  6:130661/146‑1 (MQ=255)
                                                                                                            aCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTg                                         >  4:222353/1‑149 (MQ=255)
                                                                                                               gagagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCt                                      <  3:391082/149‑1 (MQ=255)
                                                                                                                agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCt                                      >  3:308749/1‑148 (MQ=255)
                                                                                                                agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTc                                     >  4:382080/1‑149 (MQ=255)
                                                                                                                agagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCg                                    >  7:128412/1‑149 (MQ=255)
                                                                                                                  agGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGTTCATAATTGAAACAAAAAGGGTAAAAAGTTTGAGCTGGCTcac                                   >  1:351691/1‑147 (MQ=255)
                                                                                                                    gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               <  6:333788/135‑1 (MQ=255)
                                                                                                                    gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                               >  5:333788/1‑135 (MQ=255)
                                                                                                                    gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGca                                  <  4:4800/148‑1 (MQ=255)
                                                                                                                    gTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGca                                  >  3:4800/1‑148 (MQ=255)
                                                                                                                         cATAGGGTGTCCTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCAGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAc                             <  8:376487/148‑1 (MQ=255)
                                                                                                                         cATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAc                             >  6:406773/1‑147 (MQ=255)
                                                                                                                          aTAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAAcc                            <  8:58797/148‑1 (MQ=255)
                                                                                                                           tAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAg                          <  8:219717/149‑1 (MQ=255)
                                                                                                                              ggTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGCTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCAc                       >  2:386163/1‑149 (MQ=255)
                                                                                                                                    cTTTAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                  >  3:92835/1‑148 (MQ=255)
                                                                                                                                    cTTTAGTAAAACCACAGTTCTTGTAAGGATATGGGGTACCTGCGGACTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttttt                  >  4:155485/1‑148 (MQ=255)
                                                                                                                                     tttAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                 >  8:401355/1‑148 (MQ=255)
                                                                                                                                     tttAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                 >  5:387957/1‑148 (MQ=255)
                                                                                                                                     tttAGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                 >  7:180655/1‑148 (MQ=255)
                                                                                                                                       tAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAGAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTaaa               <  1:60317/148‑1 (MQ=255)
                                                                                                                                        aGTAAAACCACAGTTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTaaat              <  7:297951/148‑2 (MQ=255)
                                                                                                                                         gTAAAACCACAGGTCTTGAAGGTTATGGGGTACCTGCGGCCTTAAAATAAAACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGttt           >  4:304944/1‑149 (MQ=255)
                                                                                                                                          tAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGCGCTGGCTCAGACAAACCAGCACTTTTTaatgggt           >  2:237454/1‑144 (MQ=255)
                                                                                                                                          tAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGttt           >  4:249637/1‑149 (MQ=255)
                                                                                                                                            aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg         <  6:89464/149‑1 (MQ=255)
                                                                                                                                            aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg         <  3:222353/149‑1 (MQ=255)
                                                                                                                                            aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTg         <  1:237454/149‑1 (MQ=255)
                                                                                                                                            aaaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGGTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGtttt          <  5:189349/148‑1 (MQ=255)
                                                                                                                                            aaaCCACAGGTCTTGTAAGGCTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTg         <  1:386163/149‑1 (MQ=255)
                                                                                                                                             aaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCTGCCTTAAAATAACACCCAGCCAACATAACATAAATGAACATGGATCATAATTGCAACAAAAAGGTTAAAAAGTTTGAGCTGTCTCAAACAAACCCGCAGTATTTAAAgagttc         >  6:372191/1‑143 (MQ=255)
                                                                                                                                               ccACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTaa      <  2:41668/149‑1 (MQ=255)
                                                                                                                                                 acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt            <  2:61250/141‑1 (MQ=255)
                                                                                                                                                 acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt            >  1:61250/1‑141 (MQ=255)
                                                                                                                                                 acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAAt     >  5:193260/1‑148 (MQ=255)
                                                                                                                                                 acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAACCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAAt     >  7:84300/1‑148 (MQ=255)
                                                                                                                                                  caGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTAAAATAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAATCa   <  5:112902/149‑1 (MQ=255)
                                                                                                                                                   aGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCACACAAACCAGCACTTTTTAAAGTTTTGTAATCAg  >  3:138248/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAG  >  NC_000913/3683324‑3683619

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