Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,245,862 (C)6→7 coding (1515/1698 nt) treA ← periplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,245,8561.C94.2% 357.3 / 15.8 103coding (1521/1698 nt)treAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (45/52);  total (48/55)
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.70e-01

Only 100 of 115 total aligned reads displayed.

TGAGGGTTGCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCG  >  NC_000913/1245706‑1246000
                                                                                                                                                       |                                                                                                                                                 
tGAGGGTTGCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCAccgccg                                                                                                                                                    >  1:388794/1‑151 (MQ=255)
tGAGGGTTGCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCAccgccg                                                                                                                                                    >  2:570716/1‑151 (MQ=255)
 gAGGGTTGCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACcgccg‑c                                                                                                                                                  >  1:126337/1‑151 (MQ=255)
 gAGGGTTGCGTCGTTGCTGACTTAACAGTAGAAAGAGAAGCCGAAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCATGTAATGAATAATCGCCACcgccg‑c                                                                                                                                                  >  1:1136356/1‑151 (MQ=255)
       tGCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑ccccccc                                                                                                                                            >  1:1169620/1‑150 (MQ=255)
        gCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑ccccccc                                                                                                                                            <  4:174438/150‑2 (MQ=255)
         cgtcgtTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑cccccccgt                                                                                                                                          <  1:1076659/151‑4 (MQ=255)
            cgtTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTcc                                                                                                                                       <  1:578446/151‑1 (MQ=255)
            cgtTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTcc                                                                                                                                       <  1:315539/151‑1 (MQ=255)
            cgtTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTcc                                                                                                                                       >  1:464427/1‑151 (MQ=255)
            cgtTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑cccccccgttc                                                                                                                                        <  3:261891/150‑6 (MQ=255)
             gtTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCAGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCC‑GTTCCgg                                                                                                                                     <  2:834850/151‑1 (MQ=255)
                gCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCggtg                                                                                                                                   <  1:993131/151‑1 (MQ=255)
                gCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCggt                                                                                                                                    <  2:611121/150‑1 (MQ=255)
                gCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCggt                                                                                                                                    <  4:449483/150‑1 (MQ=255)
                  tGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCggtggt                                                                                                                                 >  1:190915/1‑151 (MQ=255)
                   gACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCggtggt                                                                                                                                 <  4:149623/150‑1 (MQ=255)
                   gACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGgtggtg                                                                                                                                <  1:346062/151‑1 (MQ=255)
                       tAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTg                                                                                                                             <  2:1236873/150‑1 (MQ=255)
                          cGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGTGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGAcat                                                                                                                         <  3:11566/151‑1 (MQ=255)
                            gTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGAcatca                                                                                                                       <  3:352009/151‑1 (MQ=255)
                             tCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGAcatcat                                                                                                                      >  4:133871/1‑151 (MQ=255)
                              cagACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGAcatcat                                                                                                                      >  2:1137346/3‑150 (MQ=255)
                               ggACGCGTCGCCGGAACATTGTCCCACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCatat                                                                                                                    <  4:299919/151‑1 (MQ=255)
                                gACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATAtt                                                                                                                   >  1:297564/1‑151 (MQ=255)
                                       cGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATCGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTcca                                                                                                             <  2:974989/150‑1 (MQ=255)
                                          cGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTccacc                                                                                                           >  4:340917/1‑149 (MQ=255)
                                           ggAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGc                                                                                                        >  1:1191432/1‑151 (MQ=255)
                                              aCATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCttt                                                                                                     >  1:246074/1‑151 (MQ=255)
                                               cATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCtttt                                                                                                    >  2:916406/1‑151 (MQ=255)
                                               cATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCtttt                                                                                                    <  2:224641/151‑1 (MQ=255)
                                               cATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCtttt                                                                                                    >  4:47728/1‑151 (MQ=255)
                                               cATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTGCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCtttt                                                                                                    >  1:537774/1‑151 (MQ=255)
                                                aTTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCttttt                                                                                                   <  1:973957/151‑1 (MQ=255)
                                                aTTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCttttt                                                                                                   >  1:865628/1‑151 (MQ=255)
                                                    tCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTcc                                                                                               >  1:1138025/1‑151 (MQ=255)
                                                        acGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGt                                                                                           <  2:35190/151‑1 (MQ=255)
                                                         cGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGt                                                                                           <  2:99720/150‑1 (MQ=255)
                                                         cGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGt                                                                                           <  1:1271123/150‑1 (MQ=255)
                                                         cGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTc                                                                                          >  3:41811/1‑151 (MQ=255)
                                                         cGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTc                                                                                          <  1:494406/151‑1 (MQ=255)
                                                            ttGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCAt                                                                                        >  2:69390/1‑150 (MQ=255)
                                                             tGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATa                                                                                       <  1:882637/150‑1 (MQ=255)
                                                                tctTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGgtg                                                                                   >  2:30943/1‑151 (MQ=255)
                                                                   ttCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGgtgt                                                                                  <  1:98754/149‑1 (MQ=255)
                                                                   ttCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGc                                                                                >  4:297491/1‑151 (MQ=255)
                                                                     tgggCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTg                                                                              >  3:376451/2‑151 (MQ=255)
                                                                       ggCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGaa                                                                            <  2:765877/151‑1 (MQ=255)
                                                                         cAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAAc                                                                           >  3:95642/1‑150 (MQ=255)
                                                                          aGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACAt                                                                         >  4:409263/1‑151 (MQ=255)
                                                                           gATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACAtt                                                                        >  3:213533/1‑151 (MQ=255)
                                                                             tCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTg                                                                       >  2:838836/1‑150 (MQ=255)
                                                                              cAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGt                                                                     >  2:17426/1‑151 (MQ=255)
                                                                              cAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGt                                                                     >  1:348157/1‑151 (MQ=255)
                                                                               aaaTCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTc                                                                    <  1:70476/151‑1 (MQ=255)
                                                                                aaTCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTc                                                                    >  2:847615/1‑150 (MQ=255)
                                                                                   ccAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAgg                                                                 >  1:216342/1‑150 (MQ=255)
                                                                                   ccAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGa                                                                <  2:728884/150‑1 (MQ=255)
                                                                                    cAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGaa                                                               <  2:278445/151‑1 (MQ=255)
                                                                                     aGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAg                                                              <  2:1196288/151‑1 (MQ=255)
                                                                                     aGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACAGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGa                                                                >  1:128920/1‑149 (MQ=255)
                                                                                        aTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGc                                                           >  3:51147/1‑151 (MQ=255)
                                                                                         ttttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGcc                                                          >  4:18893/1‑151 (MQ=255)
                                                                                          tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCCAGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGcc                                                          >  2:1099863/1‑150 (MQ=255)
                                                                                           ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATGCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCa                                                         <  1:1276937/150‑1 (MQ=255)
                                                                                           ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGcc                                                          >  2:304525/1‑149 (MQ=255)
                                                                                           ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGcc                                                          >  4:289223/1‑149 (MQ=255)
                                                                                           ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCa                                                         <  3:138348/150‑1 (MQ=255)
                                                                                           ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCa                                                         <  3:409257/150‑1 (MQ=255)
                                                                                            tCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAg                                                        <  2:263467/150‑1 (MQ=255)
                                                                                            tCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAg                                                        <  1:886255/150‑1 (MQ=255)
                                                                                            tCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGc                                                       >  1:868947/1‑151 (MQ=255)
                                                                                             cAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGc                                                       <  4:242537/150‑1 (MQ=255)
                                                                                               gCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTaa                                                    <  3:404329/151‑1 (MQ=255)
                                                                                                   cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTc                                                >  1:1034272/1‑151 (MQ=255)
                                                                                                   cACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTc                                                <  1:1185051/151‑1 (MQ=255)
                                                                                                   cACGACATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCCCCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGt                                                 <  2:209527/150‑1 (MQ=255)
                                                                                                     cGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCa                                              <  1:50433/151‑1 (MQ=255)
                                                                                                      gCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCAt                                             >  3:450546/1‑151 (MQ=255)
                                                                                                         aTTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGc                                          >  2:820594/1‑151 (MQ=255)
                                                                                                               ccAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACctc                                    >  1:205791/1‑151 (MQ=255)
                                                                                                                 aGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCtt                                  <  2:308718/151‑1 (MQ=255)
                                                                                                                  gCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCttt                                 <  4:393759/151‑1 (MQ=255)
                                                                                                                    cAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTg                              >  3:128783/1‑151 (MQ=255)
                                                                                                                     aaaGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTg                              >  2:601266/1‑151 (MQ=255)
                                                                                                                       aGCCATCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGcc                            >  2:903748/1‑151 (MQ=255)
                                                                                                                           aTCCTGTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGt                         >  4:194113/1‑150 (MQ=255)
                                                                                                                                gTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGttt                    <  2:297520/150‑1 (MQ=255)
                                                                                                                                gTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGttt                    <  2:190889/150‑1 (MQ=255)
                                                                                                                                gTAATGGATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGtt                     <  1:570827/149‑1 (MQ=255)
                                                                                                                                   aTGGATATTCGCCACCGCCG‑CCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTaa               <  2:828609/151‑1 (MQ=255)
                                                                                                                                     ggATATTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAAt              <  2:1212685/151‑1 (MQ=255)
                                                                                                                                       atatTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATcc            >  2:1015989/1‑151 (MQ=255)
                                                                                                                                         atTCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCtt          <  2:127161/151‑1 (MQ=255)
                                                                                                                                         atTCGCCACCGCCG‑CCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTc         >  2:1000191/1‑151 (MQ=255)
                                                                                                                                          ttCGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTc         >  4:335500/1‑151 (MQ=255)
                                                                                                                                            cGCCACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTg       <  3:194110/151‑1 (MQ=255)
                                                                                                                                              ccACCGCCGCCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTc     >  4:5663/1‑151 (MQ=255)
                                                                                                                                                 ccgccgcCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTcgcg  <  2:398219/151‑1 (MQ=255)
                                                                                                                                                  cgccgcCCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTcgcg  >  3:72704/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TGAGGGTTGCGTCGTTGCTGACTTAACGGTCGGACGCGTCGCCGGAACATTGTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑CCCCCC‑GTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCG  >  NC_000913/1245706‑1246000

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