Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,252 T→A 8.9% intergenic (+16/+216) rhlE → / ← ybiA ATP‑dependent RNA helicase/DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,2520TA 8.9% 278.5 / 14.4 90intergenic (+16/+216)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base A (4/4);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAG  >  NC_000913/832104‑832399
                                                                                                                                                    |                                                                                                                                                   
ccGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCt                                                                                                                                                     <  3:379418/149‑1 (MQ=255)
  gTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                   <  5:81218/149‑1 (MQ=255)
     gggAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGc                                                                                                                                                 <  2:259456/148‑1 (MQ=255)
      ggAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGccc                                                                                                                                               >  8:178286/1‑149 (MQ=255)
       gAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              >  2:78163/1‑149 (MQ=255)
       gAAGGTGGCGCAAAATCCGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              >  6:121078/1‑149 (MQ=255)
         aGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGc                                                                                                                                            <  3:231094/149‑1 (MQ=255)
          ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                   >  1:144528/1‑141 (MQ=255)
          ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                   <  2:144528/141‑1 (MQ=255)
          ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              >  5:238222/1‑146 (MQ=255)
          ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              <  6:238222/146‑1 (MQ=255)
            tGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATc                                                                                                                                         <  5:346292/149‑1 (MQ=255)
              gcgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                        <  3:429398/148‑1 (MQ=255)
              gcgcAAAATCAGCAAGCGCGAAACCTCCAGAAAAACCGTCTCGCCGCCTCGCCGATGCCAAACCCGCTGGCGCACAACAACCTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAg                                                                                                                                       >  5:100319/1‑149 (MQ=255)
               cgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              <  2:179214/141‑1 (MQ=255)
               cgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              >  1:179214/1‑141 (MQ=255)
                gcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGc                                                                                                                                     <  5:204982/149‑1 (MQ=255)
                     aTCTGCAATGGCGAAACCTGCAGATAAACGGTCTCGCCGCCTCGGCGATGCCTAACCGGCAGGCGAAGAACAACGTCGGCGCCGTCCGCGTATACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtgatg                                                                                                                                <  6:2773/149‑1 (MQ=255)
                        tGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           >  2:50546/1‑135 (MQ=255)
                        tGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           <  1:50546/135‑1 (MQ=255)
                          cAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           >  8:213521/1‑133 (MQ=255)
                          cAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           <  7:213521/133‑1 (MQ=255)
                           aaGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaaa                                                                                                                                      <  2:148614/137‑5 (MQ=255)
                           aaGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaaa                                                                                                                                      >  1:148614/1‑133 (MQ=255)
                               gcgAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCa                                                                                                                      <  7:178286/149‑1 (MQ=255)
                                 gAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATCAACAAACGCaaa                                                                                                                    <  1:172052/149‑1 (MQ=255)
                                    aCCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACt                                                                                                                 <  7:177007/149‑1 (MQ=255)
                                      cTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                               <  3:92891/149‑1 (MQ=255)
                                       tGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                  <  8:390637/129‑1 (MQ=255)
                                       tGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                  >  7:390637/1‑129 (MQ=255)
                                       tGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                              >  3:249205/1‑149 (MQ=255)
                                       tGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                               >  4:222057/1‑148 (MQ=255)
                                         cAGATAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCTGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                              <  6:100319/147‑1 (MQ=255)
                                              aaaCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATgcg                                                                                                       >  4:130901/1‑149 (MQ=255)
                                                    tctcGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACga                                                                                                 <  7:330559/149‑2 (MQ=255)
                                                        gccgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtgat                                                                                                                                 <  6:11583/113‑1 (MQ=255)
                                                        gccgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtgat                                                                                                                                 >  5:11583/1‑113 (MQ=255)
                                                         ccgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTAt                                                                                             >  7:46161/1‑148 (MQ=255)
                                                          cgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTAt                                                                                             >  6:437174/1‑147 (MQ=255)
                                                              tCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                       >  7:202714/1‑149 (MQ=255)
                                                              tCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                       >  8:143610/1‑149 (MQ=255)
                                                                ggCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                       <  5:121078/147‑1 (MQ=255)
                                                                 gCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                    <  6:444835/149‑1 (MQ=255)
                                                                   gATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                  <  5:306380/149‑1 (MQ=255)
                                                                        cAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCaaaa                                                                                                                   <  7:348659/111‑1 (MQ=255)
                                                                        cAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCaaaa                                                                                                                   >  8:348659/1‑111 (MQ=255)
                                                                         aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAActc                                                                            >  1:141392/1‑149 (MQ=255)
                                                                         aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAATAATCTTTTATGCCGGGCTATGCCCGGCa                                                                                                                                           <  4:419056/86‑1 (MQ=255)
                                                                         aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAATAATCTTTTATGCCGGGCTATGCCCGGCa                                                                                                                                           >  3:419056/1‑86 (MQ=255)
                                                                                 aGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAata                                                                     >  2:273445/1‑148 (MQ=255)
                                                                                  ggCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAatata                                                                   <  1:390235/149‑1 (MQ=255)
                                                                                   gCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAatatat                                                                  >  8:416716/1‑149 (MQ=255)
                                                                                    cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                  >  6:150508/1‑132 (MQ=255)
                                                                                    cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                                  <  5:150508/132‑1 (MQ=255)
                                                                                    cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGAAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                 <  3:222057/149‑1 (MQ=255)
                                                                                     gAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                 >  4:17794/1‑148 (MQ=255)
                                                                                      aacaacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAg                                                                                                                                       >  5:452/1‑77 (MQ=255)
                                                                                      aacaacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAg                                                                                                                                       <  6:452/77‑1 (MQ=255)
                                                                                        caacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAt                                                             <  5:437174/149‑1 (MQ=255)
                                                                                         aacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAtt                                                            >  7:86883/1‑149 (MQ=255)
                                                                                         aacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaa                                                                                                                                       <  1:293217/74‑4 (MQ=38)
                                                                                         aacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaa                                                                                                                                       >  2:293217/1‑71 (MQ=38)
                                                                                            aacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGc                                                         <  3:130901/149‑1 (MQ=255)
                                                                                             acGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCa                                                        <  1:78163/149‑1 (MQ=255)
                                                                                                 cgccgcCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCt                                                    >  4:69181/1‑149 (MQ=255)
                                                                                                 cgccgcCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCt                                                    >  4:63187/1‑149 (MQ=255)
                                                                                                   ccgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCttt                                                  >  6:299033/1‑149 (MQ=255)
                                                                                                    cgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCtttt                                                 <  8:46161/149‑1 (MQ=255)
                                                                                                    cgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCtttt                                                 <  7:42296/149‑1 (MQ=255)
                                                                                                       cgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTa                                              <  2:141392/149‑1 (MQ=255)
                                                                                                        gTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAg                                             <  3:69181/149‑1 (MQ=255)
                                                                                                        gTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAg                                             <  3:63187/149‑1 (MQ=255)
                                                                                                        gTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGAATGCTTTTGTAg                                             >  4:301534/1‑149 (MQ=255)
                                                                                                           cgcgTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGc                                           <  4:165080/148‑1 (MQ=255)
                                                                                                            gcgTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCg                                         >  7:325988/1‑149 (MQ=255)
                                                                                                              gTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGa                                       >  1:394877/1‑149 (MQ=255)
                                                                                                              gTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGa                                       <  7:416716/149‑1 (MQ=255)
                                                                                                                 aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATaa                                    <  3:17794/149‑1 (MQ=255)
                                                                                                                 aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATaa                                    <  8:325988/149‑1 (MQ=255)
                                                                                                                  aCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAg                                   >  1:69174/1‑149 (MQ=255)
                                                                                                                    cTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGc                                 <  4:249205/149‑1 (MQ=255)
                                                                                                                             cgcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAAcgcc                        >  6:299591/1‑149 (MQ=255)
                                                                                                                             cgcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                         <  7:143610/148‑1 (MQ=255)
                                                                                                                                 gTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                    >  7:442287/1‑149 (MQ=255)
                                                                                                                                 gTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                    >  7:8446/1‑149 (MQ=255)
                                                                                                                                    aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgcc                        >  5:237703/1‑142 (MQ=255)
                                                                                                                                    aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGAGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgcc                        <  6:237703/142‑1 (MQ=255)
                                                                                                                                     tCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                >  3:243439/1‑149 (MQ=255)
                                                                                                                                       ttttATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa              >  4:316292/1‑149 (MQ=255)
                                                                                                                                         ttATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc               <  5:299591/146‑1 (MQ=21)
                                                                                                                                             gCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTca         >  3:68143/1‑148 (MQ=255)
                                                                                                                                                    tATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAg  >  5:259363/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAG  >  NC_000913/832104‑832399

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