Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,978,113 G→A 5.1% E197K (GAG→AAG)  wecF → TDP‑Fuc4NAc:lipidIIFuc4NAc transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,978,1130GA 5.1% 322.2 / 4.3 101E197K (GAG→AAG) wecFTDP‑Fuc4NAc:lipidIIFuc4NAc transferase
Reads supporting (aligned to +/- strand):  ref base G (42/52);  new base A (3/2);  total (45/56)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

Only 100 of 101 total aligned reads displayed.

TTTTTGCCAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCAGCGAAGAAAATCTAC  >  NC_000913/3977967‑3978260
                                                                                                                                                  |                                                                                                                                                   
tttttGCCAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAg                                                                                                                                                   >  5:184701/1‑149 (MQ=255)
 ttttGCCAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGa                                                                                                                                                    <  2:10519/147‑1 (MQ=255)
 ttttGCCAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAg                                                                                                                                                   <  5:140752/148‑1 (MQ=255)
  tttGCCAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGc                                                                                                                                                  <  5:342764/148‑1 (MQ=255)
       cAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATtg                                                                                                                                            >  6:116428/1‑149 (MQ=255)
        aaaaCGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATtgc                                                                                                                                           >  5:22890/1‑149 (MQ=255)
         aaaCGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATtgct                                                                                                                                          <  8:71989/149‑1 (MQ=255)
            cGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGcc                                                                                                                                       <  2:297521/149‑1 (MQ=255)
            cGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGcc                                                                                                                                       <  4:320942/149‑1 (MQ=255)
              cACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCtt                                                                                                                                     <  1:110148/149‑1 (MQ=255)
                cccAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTgc                                                                                                                                   >  5:122385/1‑149 (MQ=255)
                cccAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATATGCTCATtgctgcggtcg                                                                                                                                   <  5:23715/149‑6 (MQ=255)
                 ccAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTgcg                                                                                                                                  >  3:323865/1‑149 (MQ=255)
                  cAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATtg                                                                                                                                            <  7:253809/138‑1 (MQ=255)
                  cAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATtg                                                                                                                                            >  8:253809/1‑138 (MQ=255)
                   aaaGGTGCGGGGCGAACTGCTGTACTTCCCGAAGCGGATGGACCCTTCGCTCAATACGATGGAGAACGAACGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTgcgcg                                                                                                                                >  5:314999/1‑149 (MQ=255)
                     aGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGcc                                                                                                                              <  7:370896/149‑1 (MQ=255)
                       gTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGt                                                                                                                            <  5:170066/149‑1 (MQ=255)
                           ggggCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcat                                                                                                                        <  4:134009/149‑1 (MQ=255)
                           ggggCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcat                                                                                                                        <  7:124243/149‑1 (MQ=255)
                            gggCGAACTGCTGTTCTTCCCGACGCGTATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATTAGCATATTGCTGCCTTGCGCGCCGTtcatc                                                                                                                       <  2:278574/149‑1 (MQ=255)
                            gggCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcatc                                                                                                                       <  8:147775/149‑1 (MQ=255)
                             ggCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcatca                                                                                                                      <  1:346288/149‑1 (MQ=255)
                               cGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGc                                                                                                                    <  5:150910/149‑1 (MQ=255)
                                gAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAGCGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtc                                                                                                                          <  8:381758/142‑1 (MQ=255)
                                gAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAGCGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtc                                                                                                                          >  7:381758/1‑142 (MQ=255)
                                    tgctgTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAAttt                                                                                                               >  5:72328/1‑149 (MQ=255)
                                      ctgTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTgg                                                                                                             <  4:323865/149‑1 (MQ=255)
                                      ctgTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTgg                                                                                                             >  5:298568/1‑149 (MQ=255)
                                      ctgTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGCTCATCAGCAATTTgg                                                                                                             <  5:116428/149‑1 (MQ=255)
                                       tgTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGc                                                                                                            >  5:187975/1‑149 (MQ=255)
                                        gTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGTGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAg                                                                                                                                                   <  6:296183/109‑1 (MQ=255)
                                        gTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGAGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAg                                                                                                                                                   >  5:296183/1‑109 (MQ=255)
                                         ttcttcCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGATCTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGc                                                                                                            <  1:61571/147‑1 (MQ=255)
                                         ttcttcCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGa                                                                                                          >  4:109753/1‑149 (MQ=255)
                                          tcttcCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGa                                                                                                          >  7:28209/1‑148 (MQ=255)
                                            ttcCCGACGCGGATGGACCCTTCGCTCAATACGCTGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATAc                                                                                                       >  4:36103/1‑149 (MQ=255)
                                                cGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGt                                                                                                    >  6:96500/1‑148 (MQ=255)
                                                  aCGCGGATGTACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAATGGAAAATGACCATTCTGGTGGTGAACTCCGGCGACCGCAGCAATGAGCATATtgctgc                                                                                                                                        <  1:99045/110‑1 (MQ=255)
                                                    gcgGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAg                                                                                               >  3:218134/1‑149 (MQ=255)
                                                      ggATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcat                                                                                                                        <  6:208987/122‑1 (MQ=255)
                                                      ggATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcat                                                                                                                        >  5:208987/1‑122 (MQ=255)
                                                      ggATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtg                                                                                             <  6:72328/149‑1 (MQ=255)
                                                      ggATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtg                                                                                             <  5:234589/149‑1 (MQ=255)
                                                       gATGGACCCTTTGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtgg                                                                                            >  7:180838/1‑149 (MQ=255)
                                                       gATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtgg                                                                                            >  1:158671/1‑149 (MQ=255)
                                                         tGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtggtg                                                                                          <  2:158729/149‑1 (MQ=255)
                                                            aCCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcatca                                                                                                                      >  1:307174/1‑118 (MQ=255)
                                                            aCCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTtcatca                                                                                                                      <  2:307174/118‑1 (MQ=255)
                                                            aCCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtggtggtg                                                                                       >  8:209158/1‑149 (MQ=255)
                                                            aCCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtggtggtg                                                                                       >  5:279984/1‑149 (MQ=255)
                                                             cccTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtggtgg                                                                                         <  2:209445/146‑1 (MQ=255)
                                                             cccTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAgtggtgg                                                                                         >  1:209445/1‑146 (MQ=255)
                                                              ccTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGcc                                                                                     >  7:142971/1‑149 (MQ=255)
                                                                ttCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCAACAGCAATTTGGCCATACGGTAAAAGTGGTGGTGCCGa                                                                                   >  8:295031/1‑149 (MQ=255)
                                                                  cGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGAt                                                                                  <  5:287970/148‑1 (MQ=255)
                                                                   gCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATgg                                                                                <  3:109753/149‑1 (MQ=255)
                                                                   gCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTGCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATgg                                                                                <  5:42804/149‑1 (MQ=255)
                                                                     tCAATACGTTGGCGAACGATCTGCATCGTGTAGGGAATATGACCATTCTGGTGGGTAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCATCAATTTGGCGATACGGTAAAAGTGGTTGTGCCGATgg                                                                                <  6:314999/147‑1 (MQ=255)
                                                                       aaTACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATAc                                                                                                       >  7:34883/1‑122 (MQ=255)
                                                                       aaTACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATAc                                                                                                       <  8:34883/122‑1 (MQ=255)
                                                                       aaTACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGata                                                                            >  7:95339/1‑149 (MQ=255)
                                                                       aaTACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGa                                                                              <  7:209158/147‑1 (MQ=255)
                                                                         tACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATc                                                                          <  6:184701/149‑1 (MQ=255)
                                                                            gATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGc                                                                       >  5:334783/1‑149 (MQ=255)
                                                                                gCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATAAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCtaa                                                                   >  3:173988/1‑149 (MQ=255)
                                                                                 cGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCtaat                                                                  <  2:158671/149‑1 (MQ=255)
                                                                                  gAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCtaat                                                                  >  6:165928/1‑148 (MQ=255)
                                                                                      gATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACg                                                              >  2:40528/1‑148 (MQ=255)
                                                                                       aTCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGaa                                                            <  1:21679/149‑1 (MQ=255)
                                                                                        tCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAg                                                           <  1:7346/149‑1 (MQ=255)
                                                                                            cAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGTATATCCGCCTAATAACGAAGCGTa                                                       <  6:319824/149‑1 (MQ=255)
                                                                                            cAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTa                                                       <  1:241091/149‑1 (MQ=255)
                                                                                             aaCGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTa                                                       >  6:37201/1‑148 (MQ=255)
                                                                                              aCGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACa                                                     <  8:28209/149‑1 (MQ=255)
                                                                                                gTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACAtt                                                   >  6:169258/1‑149 (MQ=255)
                                                                                                   aaGGGAAAATGACCATTCTGGTGGGTAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAg                                                <  6:22890/149‑1 (MQ=255)
                                                                                                   aaGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGa                                                 <  5:165928/148‑1 (MQ=255)
                                                                                                    aGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAgg                                               <  5:37201/149‑1 (MQ=255)
                                                                                                      ggAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTTGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGaa                                             >  2:208895/1‑149 (MQ=255)
                                                                                                      ggAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATCAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGaa                                             <  2:299059/149‑1 (MQ=255)
                                                                                                        aaaaTGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGt                                           >  5:159643/1‑149 (MQ=255)
                                                                                                         aaaTGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGtt                                          <  4:218134/149‑1 (MQ=255)
                                                                                                          aaTGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTc                                         <  6:122385/149‑1 (MQ=255)
                                                                                                          aaTCACCATTCTGGTGGGGAACTCCGGCGAGCGCAGCTATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAATGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGttg                                         <  7:295031/149‑2 (MQ=255)
                                                                                                                cATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAgg                                   >  2:113739/1‑149 (MQ=255)
                                                                                                                    cTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCgg                                <  6:187975/148‑1 (MQ=255)
                                                                                                                          gggAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTg                           >  2:120633/1‑147 (MQ=255)
                                                                                                                          gggAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGa                         <  3:36103/149‑1 (MQ=255)
                                                                                                                           ggAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTTCGCGCAGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAg                        >  6:271010/1‑149 (MQ=255)
                                                                                                                            gAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGt                       <  6:279984/149‑1 (MQ=255)
                                                                                                                             aaCTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGt                       >  3:63177/1‑148 (MQ=255)
                                                                                                                                  cGGCGACCGCAGCAATGAGCATATTGCTGCATTGCGCGCCGTTCATCAGCAATGTGGAGATCCGGTAAAAGTGGTGGTGCGGATGGGATATCCGCCTAATAACGAAGCGTAAATTGAGGAAGTTCGTCAGTCGGTGCTGGAGTTAtt                   >  2:136342/1‑147 (MQ=255)
                                                                                                                                  cGGCGACCGCAGCAATAAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCa                 <  4:173988/149‑1 (MQ=255)
                                                                                                                                         ggcagcaATATGCTCATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCAGCgaag           >  4:134231/2‑148 (MQ=255)
                                                                                                                                         ggcagcaATATGCTCATAGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCAGCgaag           >  4:125176/2‑148 (MQ=255)
                                                                                                                                            agcaATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCGGGAGTTATTCAGCGAAGaaaa       <  5:266724/149‑1 (MQ=255)
                                                                                                                                              caATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGCGAAAAGTGGGGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCAGCGAAGAAAATc     >  3:208782/1‑149 (MQ=255)
                                                                                                                                                 tGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCTATACGGTAAAAGTGGTGGGGCCTATGGGATCTCCGCCTAATAACGATGTGTACATTGATGAAGTTCGTCAGGCGGGTCTGGAGTTATTCAGCGATGAAAATCTAc  <  1:136342/149‑1 (MQ=255)
                                                                                                                                                 tGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCAGCGAAGAAAATCTa   >  5:134082/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTTTTGCCAAAACGCACCCAAAGGTGCGGGGCGAACTGCTGTTCTTCCCGACGCGGATGGACCCTTCGCTCAATACGATGGCGAACGATCGGCAACGTGAAGGGAAAATGACCATTCTGGTGGGGAACTCCGGCGACCGCAGCAATGAGCATATTGCTGCCTTGCGCGCCGTTCATCAGCAATTTGGCGATACGGTAAAAGTGGTGGTGCCGATGGGATATCCGCCTAATAACGAAGCGTACATTGAGGAAGTTCGTCAGGCGGGGCTGGAGTTATTCAGCGAAGAAAATCTAC  >  NC_000913/3977967‑3978260

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