Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,253 A→T 8.7% intergenic (+17/+215) 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,2530AT 8.7% 281.1 / 12.5 91intergenic (+17/+215)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base T (4/4);  total (44/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 = 9.73e-01

GTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAG  >  NC_000913/832106‑832399
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                                              >  1:179214/1‑141 (MQ=255)
             cgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                              <  2:179214/141‑1 (MQ=255)
              gcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGc                                                                                                                                     <  5:204982/149‑1 (MQ=255)
                   aTCTGCAATGGCGAAACCTGCAGATAAACGGTCTCGCCGCCTCGGCGATGCCTAACCGGCAGGCGAAGAACAACGTCGGCGCCGTCCGCGTATACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtgatg                                                                                                                                <  6:2773/149‑1 (MQ=255)
                      tGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           <  1:50546/135‑1 (MQ=255)
                      tGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           >  2:50546/1‑135 (MQ=255)
                        cAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           >  8:213521/1‑133 (MQ=255)
                        cAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                           <  7:213521/133‑1 (MQ=255)
                         aaGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaaa                                                                                                                                      >  1:148614/1‑133 (MQ=255)
                         aaGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCAtaaaa                                                                                                                                      <  2:148614/137‑5 (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                                                                                                                   >  8:348659/1‑111 (MQ=255)
                                                                      cAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCaaaa                                                                                                                   <  7:348659/111‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
GTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAG  >  NC_000913/832106‑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: