Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,274,117 G→A 6.0% A112T (GCC→ACC)  ychO → putative invasin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,274,1170GA 6.0% 209.5 / 2.3 66A112T (GCC→ACC) ychOputative invasin
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TGAAAGATTTCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTGGTGAGCAAT  >  NC_000913/1273971‑1274260
                                                                                                                                                  |                                                                                                                                               
tgAAAGATTTCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGc                                                                                                                                                >  7:367058/1‑148 (MQ=255)
   aaGATTTCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAgt                                                                                                                                            <  4:157523/149‑1 (MQ=255)
   aaGATTTCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAgt                                                                                                                                            >  8:278265/1‑149 (MQ=255)
     gATTTCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAgtgt                                                                                                                                          <  4:129518/149‑1 (MQ=255)
       tttCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGtt                                                                                                                                         <  7:372162/148‑1 (MQ=255)
         tCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGa                                                                                                                                       <  8:37842/148‑1 (MQ=255)
            gCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGTAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTc                                                                                                                                   <  4:83682/149‑1 (MQ=255)
                  ccAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTg                                                                                                                             <  8:367058/149‑1 (MQ=255)
                     gTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGAt                                                                                                                          <  6:62108/149‑1 (MQ=255)
                         gaatgaTAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAAc                                                                                                                       >  3:149674/1‑148 (MQ=255)
                         gaatgaTAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGAACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAAc                                                                                                                       >  3:327363/1‑148 (MQ=255)
                                  cGGGCTGGGTACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGTAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACAtt                                                                                                             >  2:30179/1‑149 (MQ=255)
                                  cGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGGTACTCAGCACGTAGAGTCCTGGCTAGCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACAtt                                                                                                             <  5:175468/149‑1 (MQ=255)
                                        ggATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCg                                                                                                       >  7:274404/1‑149 (MQ=255)
                                         gATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCgg                                                                                                      >  2:58108/1‑149 (MQ=255)
                                             cTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGt                                                                                                  >  7:155089/1‑149 (MQ=255)
                                              tGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAgtg                                                                                                                                           <  2:269682/107‑1 (MQ=255)
                                              tGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAgtg                                                                                                                                           >  1:269682/1‑107 (MQ=255)
                                                gCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCCGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGt                                                                                               >  5:10175/1‑149 (MQ=255)
                                                gCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGt                                                                                               <  5:350993/149‑1 (MQ=255)
                                                  gAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTgg                                                                                             <  4:225112/149‑1 (MQ=255)
                                                           aaaGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGtt                                                                                    >  6:283835/1‑149 (MQ=255)
                                                                tttCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGc                                                                               >  7:223340/1‑149 (MQ=255)
                                                                      aTTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTa                                                                          >  1:97280/1‑148 (MQ=255)
                                                                      aTTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTAc                                                                         <  8:223340/149‑1 (MQ=255)
                                                                           gAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGa                                                                     >  5:260948/1‑148 (MQ=255)
                                                                             aaaGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATaa                                                                  <  7:278265/149‑1 (MQ=255)
                                                                                gTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAg                                                                      <  2:65686/142‑1 (MQ=255)
                                                                                gTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAg                                                                      >  1:65686/1‑142 (MQ=255)
                                                                                gTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGa                                                               >  8:119025/1‑149 (MQ=255)
                                                                                 tCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGCAGGACATTTCACCGGCAGTCGTGGAAGCTGGTGTGTGCCGTTACAAGATAATGAt                                                              >  2:339108/1‑149 (MQ=255)
                                                                                   cgcgACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCg                                                            <  5:283835/149‑1 (MQ=255)
                                                                                      gACGCGCTTAGTCAACAGGTTAATCCGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATACCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCCCCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGAGAATGAGCGTTa                                                         >  1:190073/1‑149 (MQ=255)
                                                                                      gACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCCCCGTGGGGAAATACCAGTGTTGACGTCAAAGTGGATACCGAAGGCCATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTa                                                         >  1:184512/1‑149 (MQ=255)
                                                                                               aGTCACCATGTTAATCATCACGTATAGTCCTGTCTATCACCGTGTCTAAAAGCCAGTGGGGACGGCAAATTGTATATCGAAGTACATTTCAACGGCATTCGTGGAAGCTGGTTTCTGCCGTTACAAGATAATGATCGTTATCTCACCt                                                 <  5:262464/148‑1 (MQ=255)
                                                                                               aGTCAACAGGGTACTCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTGGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCAcc                                                  <  8:343663/147‑1 (MQ=255)
                                                                                                 tCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGa                                              <  8:155089/149‑1 (MQ=255)
                                                                                                 tCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGa                                              <  8:144400/149‑1 (MQ=255)
                                                                                                   aaCAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTCTCTCACCTGGAGc                                            >  3:131674/1‑149 (MQ=255)
                                                                                                   aaCAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGc                                            >  1:47490/1‑149 (MQ=255)
                                                                                                   aaCAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGGGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGc                                            <  8:292389/149‑1 (MQ=255)
                                                                                                     cAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCa                                          <  6:260948/149‑1 (MQ=255)
                                                                                                         ttAATCAGCACGTAGAGTCCAGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCtt                                      >  6:56935/1‑149 (MQ=255)
                                                                                                           aaTCAGCAAGTCGAGTCATTGCTATCACCGTGTGGAAATACCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTTTAGCCAGCTTg                                     <  2:184512/148‑1 (MQ=255)
                                                                                                           aaTCAGCAAGTAGTGTACTGGCTATCACCGTTGGGAAATACCAGTGTTGACGTCAAAGTGGATCACGAAGGACATTTCCCCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTg                                     <  2:190073/148‑1 (MQ=255)
                                                                                                             tCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTc                                  >  7:257389/1‑149 (MQ=255)
                                                                                                             tCAGCACGTAGAGTCCGTGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTc                                  <  4:327363/149‑1 (MQ=255)
                                                                                                                gCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTa                               <  4:278057/149‑1 (MQ=255)
                                                                                                                      agagTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTcag                          >  8:258093/1‑148 (MQ=255)
                                                                                                                        agTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGa                                              <  3:68847/126‑1 (MQ=255)
                                                                                                                        agTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGa                                              >  4:68847/1‑126 (MQ=255)
                                                                                                                        agTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTcagcag                       <  8:257389/149‑1 (MQ=255)
                                                                                                                          tCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGcc                                           <  7:237747/127‑1 (MQ=255)
                                                                                                                          tCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGcc                                           >  8:237747/1‑127 (MQ=255)
                                                                                                                          tCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGa                     <  6:10175/149‑1 (MQ=255)
                                                                                                                            cTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTg                                                >  8:123624/1‑120 (MQ=255)
                                                                                                                            cTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTg                                                <  7:123624/120‑1 (MQ=255)
                                                                                                                            cTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGAt                    >  4:313140/1‑148 (MQ=255)
                                                                                                                             tGGCTATCACCGTGGGGAAATGCCACTGTTGACGTCAAAGTGGAGAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTGGTGCCGTTACAAGATCAAGATCGGTATCTACCCAGGAGCCAGCTTGGTCTTACTCAGCAGGATaa                  >  2:264412/1‑149 (MQ=255)
                                                                                                                              ggCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAAt                 >  7:10194/1‑149 (MQ=255)
                                                                                                                                   tCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATggg              <  4:149674/147‑1 (MQ=255)
                                                                                                                                    cACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTg           >  1:194128/1‑149 (MQ=255)
                                                                                                                                      ccGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTGGt         <  2:47490/149‑1 (MQ=255)
                                                                                                                                      ccGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTGGt         <  5:56935/149‑1 (MQ=255)
                                                                                                                                        gTGGGGAAATGCCAGTGTTGCCGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTGGTGa       <  6:387774/149‑1 (MQ=255)
                                                                                                                                             gAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTGGTGAGCAAt  <  7:258093/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGAAAGATTTCGGCGAAACCAGTATGAATGATAACGGGCTGGATACTGGCGAGCAGGCAAAAGCTTTCGCATTGGGAAAAGTCCGCGACGCGCTTAGTCAACAGGTTAATCAGCACGTAGAGTCCTGGCTATCACCGTGGGGAAATGCCAGTGTTGACGTCAAAGTGGATAACGAAGGACATTTCACCGGCAGTCGTGGAAGCTGGTTTGTGCCGTTACAAGATAATGATCGTTATCTCACCTGGAGCCAGCTTGGTCTTACTCAGCAGGATAATGGGTTGGTGAGCAAT  >  NC_000913/1273971‑1274260

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: ‑