Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,903,497 T→C I139T (ATC→ACC)  manY → mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,903,4970TC98.3% 225.4 / ‑4.0 60I139T (ATC→ACC) manYmannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/32);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

AGAGCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGC  >  NC_000913/1903350‑1903642
                                                                                                                                                   |                                                                                                                                                 
agagCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATcc                                                                                                                                                 >  2:698398/1‑150 (MQ=255)
 gagCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCAc                                                                                                                                               >  2:81619/1‑151 (MQ=255)
  agCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACg                                                                                                                                              <  2:1019243/151‑1 (MQ=255)
     aTTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGttt                                                                                                                                           >  2:328331/1‑151 (MQ=255)
            cAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTctct                                                                                                                                    <  2:801564/151‑1 (MQ=255)
            cAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTctct                                                                                                                                    <  2:258202/151‑1 (MQ=255)
               gTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGtt                                                                                                                                 <  2:823725/151‑1 (MQ=255)
                tATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTc                                                                                                                                >  2:158832/1‑151 (MQ=255)
                tATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTc                                                                                                                                >  1:717833/1‑151 (MQ=255)
                   cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCt                                                                                                                              <  1:811631/150‑1 (MQ=255)
                   cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCAACATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCt                                                                                                                              <  1:383179/150‑1 (MQ=255)
                    gCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgc                                                                                                                            >  1:255306/1‑151 (MQ=255)
                        tGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagc                                                                                                                        >  1:840829/1‑151 (MQ=255)
                         ggCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaag                                                                                                                         <  1:188723/149‑1 (MQ=255)
                         ggCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaag                                                                                                                         <  2:174433/149‑1 (MQ=255)
                           cAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAAt                                                                                                                     <  2:76649/151‑1 (MQ=255)
                                    tGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGc                                                                                                            >  2:499465/1‑151 (MQ=255)
                                     ggCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCt                                                                                                           <  1:811918/151‑1 (MQ=255)
                                     ggCCGCTGCGGGCCAGGTACTGACCATCAGCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCt                                                                                                           <  1:413047/151‑1 (MQ=255)
                                            gCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTcc                                                                                                      <  2:439331/149‑1 (MQ=255)
                                            gCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCg                                                                                                     <  2:229510/150‑1 (MQ=255)
                                             cGGGCCAGGTACTGACCATCATCGTTCGTTCTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGc                                                                                                   <  1:480060/151‑1 (MQ=255)
                                             cGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGc                                                                                                   >  1:310457/1‑151 (MQ=255)
                                             cGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGc                                                                                                   >  2:779361/1‑151 (MQ=255)
                                                gCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGt                                                                                                >  1:654767/1‑151 (MQ=255)
                                                gCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGt                                                                                                <  1:994708/151‑1 (MQ=255)
                                                gCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGt                                                                                                >  1:190225/1‑151 (MQ=255)
                                                 ccAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTc                                                                                               >  2:177319/1‑151 (MQ=255)
                                                       aCTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGtt                                                                                         >  2:277823/1‑151 (MQ=255)
                                                        cTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTg                                                                                        >  1:999716/1‑151 (MQ=255)
                                                           aCCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcgc                                                                                     >  1:151597/1‑151 (MQ=255)
                                                            ccATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGct                                                                                    <  2:565588/151‑1 (MQ=255)
                                                             catcatCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctg                                                                                   <  1:992181/151‑1 (MQ=255)
                                                              atcatcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgt                                                                                  <  2:478690/151‑1 (MQ=255)
                                                                catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtct                                                                                <  1:433383/151‑1 (MQ=255)
                                                                catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtct                                                                                <  1:657260/151‑1 (MQ=255)
                                                                   cgttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGtt                                                                             <  1:641969/151‑1 (MQ=255)
                                                                             aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAg                                                                     <  1:562921/149‑1 (MQ=255)
                                                                                aCCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGaa                                                                 <  2:689892/150‑1 (MQ=255)
                                                                                     tGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACa                                                            <  1:707871/150‑1 (MQ=255)
                                                                                         ttCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGaa                                                         <  1:678094/149‑1 (MQ=255)
                                                                                              gCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                  <  2:221633/151‑1 (MQ=255)
                                                                                                 cGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGaa                                                <  1:830345/150‑1 (MQ=255)
                                                                                                 cGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGa                                                 >  2:559800/1‑149 (MQ=255)
                                                                                                     gCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCg                                            <  2:840918/150‑1 (MQ=255)
                                                                                                        gATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTc                                        <  1:417800/151‑1 (MQ=255)
                                                                                                        gATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTc                                        <  1:909290/151‑1 (MQ=255)
                                                                                                           aaGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                     >  2:563931/1‑151 (MQ=255)
                                                                                                             ggCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                    >  2:636042/1‑150 (MQ=255)
                                                                                                                 gctgATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgt                                 <  1:225370/149‑1 (MQ=255)
                                                                                                                  ctgATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                              >  1:88958/1‑151 (MQ=255)
                                                                                                                     aTAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                           >  2:1021018/1‑151 (MQ=255)
                                                                                                                             aaCCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                   >  1:49934/1‑151 (MQ=255)
                                                                                                                               ccTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTg                 >  2:221660/1‑151 (MQ=255)
                                                                                                                                   aCAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat              >  2:470215/1‑150 (MQ=255)
                                                                                                                                    cAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAatat            >  2:439667/1‑151 (MQ=255)
                                                                                                                                        gATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc        >  2:220942/1‑151 (MQ=255)
                                                                                                                                          tttCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCgg      >  2:825293/1‑151 (MQ=255)
                                                                                                                                            tCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGt     <  1:580674/150‑1 (MQ=255)
                                                                                                                                               tGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc  >  1:568141/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGAGCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGC  >  NC_000913/1903350‑1903642

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