Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,186 C→T 5.3% G225D (GGC→GAC)  cmtA ← putative mannitol‑specific PTS IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,1860CT 5.3% 382.4 / 8.1 113G225D (GGC→GAC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base C (50/57);  new base T (3/3);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 120 total aligned reads displayed.

ATGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAAC  >  NC_000913/3078039‑3078329
                                                                                                                                                   |                                                                                                                                               
aTGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGcc                                                                                                                                                >  3:435565/1‑149 (MQ=255)
aTGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGcc                                                                                                                                                >  3:392862/1‑149 (MQ=255)
     aataatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGccc                                                                                                                                               >  2:105371/1‑145 (MQ=255)
     aataatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGccc                                                                                                                                               <  1:105371/145‑1 (MQ=255)
     aataatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGAc                                                                                                                                           >  4:374294/1‑149 (MQ=255)
     aataatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGAc                                                                                                                                           <  3:793233/149‑1 (MQ=255)
     aataatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGAc                                                                                                                                           <  3:28690/149‑1 (MQ=255)
         atGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                       <  2:58144/149‑1 (MQ=255)
          tGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTggg                                                                                                                                      >  2:700540/1‑149 (MQ=255)
           ggTCAGCGGCTTCATCAGCTCATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGt                                                                                                                                     <  4:544990/149‑1 (MQ=255)
           ggTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTgggg                                                                                                                                     <  1:271996/149‑2 (MQ=255)
           ggTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGt                                                                                                                                     <  1:96680/149‑1 (MQ=255)
             tCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTa                                                                                                                                   <  3:38013/149‑1 (MQ=255)
               aGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTaga                                                                                                                                 <  1:81988/149‑1 (MQ=255)
                gCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTagag                                                                                                                                >  2:782175/1‑149 (MQ=255)
                    cTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                                                             <  3:87302/148‑1 (MQ=255)
                     ttCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                                                                           >  1:536752/1‑149 (MQ=255)
                     ttCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                                                                           >  3:419102/1‑149 (MQ=255)
                        atcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTg                                                                                                                                        >  2:326331/1‑133 (MQ=255)
                        atcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTg                                                                                                                                        <  1:326331/133‑1 (MQ=255)
                        atcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCAt                                                                                                                        >  3:429327/1‑149 (MQ=255)
                         tcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATa                                                                                                                       >  2:402498/1‑149 (MQ=255)
                          caGCACATACGGGAAATACAGTTCGAGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaa                                                                                                                      <  1:61681/149‑1 (MQ=255)
                              acaTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaagaa                                                                                                                  <  2:605604/149‑1 (MQ=255)
                              acaTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaagaa                                                                                                                  <  1:745845/149‑1 (MQ=255)
                                aTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                                                                    <  1:83993/129‑3 (MQ=255)
                                aTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                                                                    >  2:83993/1‑127 (MQ=255)
                                aTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAagaaga                                                                                                                <  1:601127/149‑1 (MQ=255)
                                 tACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATGAAGAAGAt                                                                                                               <  2:136370/149‑1 (MQ=255)
                                  aCGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGccc                                                                                                                                               <  3:567328/116‑1 (MQ=255)
                                  aCGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGccc                                                                                                                                               >  4:567328/1‑116 (MQ=255)
                                    gggAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGa                                                                                                            <  1:75665/149‑1 (MQ=255)
                                     ggAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGAt                                                                                                           <  1:350349/149‑1 (MQ=255)
                                         aTACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTAc                                                                                                       <  4:337818/149‑1 (MQ=255)
                                             aGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGtt                                                                                                   >  1:645518/1‑149 (MQ=255)
                                             aGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGtt                                                                                                   >  4:365423/1‑149 (MQ=255)
                                             aGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGtt                                                                                                   >  3:452072/1‑149 (MQ=255)
                                               ttCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                 >  1:108633/1‑149 (MQ=255)
                                               ttCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                 >  2:176354/1‑149 (MQ=255)
                                               ttCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                 >  1:37835/1‑149 (MQ=255)
                                                tCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                                                                <  4:419102/149‑1 (MQ=255)
                                                   tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGa                                                                                             <  4:435565/149‑1 (MQ=255)
                                                    ggATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGaa                                                                                            <  4:392862/149‑1 (MQ=255)
                                                     gATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAg                                                                                           >  2:676167/1‑149 (MQ=255)
                                                         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTg                                                                                       <  4:452072/149‑1 (MQ=255)
                                                          cACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGt                                                                                                                                     >  3:394835/1‑102 (MQ=255)
                                                           aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGtt                                                                                     <  1:532068/149‑1 (MQ=255)
                                                           aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGtt                                                                                     <  2:108633/149‑1 (MQ=255)
                                                                ggAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCAtt                                                                                >  2:341614/1‑149 (MQ=255)
                                                                ggAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCAtt                                                                                <  1:402498/149‑1 (MQ=255)
                                                                 gAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTc                                                                               <  1:289610/149‑1 (MQ=255)
                                                                  aaGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGTGTTAGAGGCCACCATaaa                                                                                                                     >  1:13368/1‑110 (MQ=255)
                                                                  aaGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCACCAGCCCCAGGCCCGGACCTGTGTTAGAGGCCACTATaaa                                                                                                                     <  2:13368/110‑1 (MQ=255)
                                                                    gTGAATAATCATCGCCTCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGAt                                                                                                               >  4:207331/1‑114 (MQ=255)
                                                                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGAt                                                                                                               <  3:207331/114‑1 (MQ=255)
                                                                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGAt                                                                                                               >  1:283294/1‑114 (MQ=255)
                                                                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGAt                                                                                                               <  2:283294/114‑1 (MQ=255)
                                                                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                 >  2:729384/1‑128 (MQ=255)
                                                                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                 <  1:729384/128‑1 (MQ=255)
                                                                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCa                                                                            >  3:793997/1‑149 (MQ=255)
                                                                      gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGc                                                                          >  3:271510/1‑149 (MQ=255)
                                                                       aataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCg                                                                         >  1:210768/1‑149 (MQ=255)
                                                                       aataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCg                                                                         <  4:271510/149‑1 (MQ=255)
                                                                        ataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGaa                                                                                            >  4:680745/1‑129 (MQ=255)
                                                                        ataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGaa                                                                                            <  3:680745/129‑1 (MQ=255)
                                                                              atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                                                                  >  4:162762/1‑149 (MQ=255)
                                                                              atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                                                                  >  1:124453/1‑149 (MQ=255)
                                                                                cGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGa                                                                <  2:210768/149‑1 (MQ=255)
                                                                                 gCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGAc                                                               <  3:603758/149‑1 (MQ=255)
                                                                                 gCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGAc                                                               <  2:43964/149‑1 (MQ=255)
                                                                                              aaCGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                                                                <  1:48522/103‑1 (MQ=255)
                                                                                              aaCGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGat                                                                      <  3:538486/129‑1 (MQ=255)
                                                                                                 gTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCa                                               <  3:374294/149‑1 (MQ=255)
                                                                                                     aCTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGa                                                                                             >  2:297686/1‑99 (MQ=255)
                                                                                                       tCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTAc                                                                                                       >  4:467592/1‑87 (MQ=255)
                                                                                                         aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                                                                >  4:787192/1‑92 (MQ=255)
                                                                                                         aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAg                                       <  3:388883/149‑1 (MQ=255)
                                                                                                         aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAg                                       <  2:319824/149‑1 (MQ=255)
                                                                                                         aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAg                                       >  2:411211/1‑149 (MQ=255)
                                                                                                            ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaaa                                                                                                                     <  1:735339/68‑1 (MQ=255)
                                                                                                            ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaaa                                                                                                                     >  2:735339/1‑68 (MQ=255)
                                                                                                                   cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACt                             >  1:491340/1‑149 (MQ=255)
                                                                                                                   cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACt                             <  4:429327/149‑1 (MQ=255)
                                                                                                                   cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACt                             >  4:60368/1‑149 (MQ=255)
                                                                                                                    aaaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACtt                            >  2:155345/1‑149 (MQ=255)
                                                                                                                     aaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataa                                                                    <  1:385311/108‑1 (MQ=255)
                                                                                                                     aaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataa                                                                    >  2:385311/1‑108 (MQ=255)
                                                                                                                             aaaCGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                        >  2:59039/1‑128 (MQ=255)
                                                                                                                             aaaCGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                        <  1:59039/128‑1 (MQ=255)
                                                                                                                                 gCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCAtt               <  3:217443/149‑1 (MQ=255)
                                                                                                                                  ccAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATtg              >  1:646157/1‑149 (MQ=255)
                                                                                                                                    agcagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATtgat            <  1:700540/149‑1 (MQ=255)
                                                                                                                                     gcagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATtgatg           >  2:460481/1‑149 (MQ=255)
                                                                                                                                       ctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgac         >  2:532063/7‑149 (MQ=255)
                                                                                                                                        tctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                        >  1:127944/6‑117 (MQ=255)
                                                                                                                                        gcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgact        >  3:470024/1‑149 (MQ=255)
                                                                                                                                         cagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGcc                       >  1:491745/1‑133 (MQ=255)
                                                                                                                                         cagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGcc                       <  2:491745/133‑1 (MQ=255)
                                                                                                                                             cccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacaa   >  3:442271/1‑149 (MQ=255)
                                                                                                                                              ccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacaac  >  1:276316/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ATGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAAC  >  NC_000913/3078039‑3078329

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: