Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,543,119 C→A 5.7% W137L (TGG→TTG)  cysP ← thiosulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,543,1190CA 5.7% 315.3 / 5.2 89W137L (TGG→TTG) cysPthiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base A (2/4);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTCGGTCACCTGG  >  NC_000913/2542970‑2543265
                                                                                                                                                     |                                                                                                                                                  
aCTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTcc                                                                                                                                                    <  1:145495/150‑1 (MQ=255)
aCTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTcc                                                                                                                                                    <  1:145528/150‑1 (MQ=255)
   gggTCATAAACTGTTCGGTTTTGCCTTTGTCCCCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAAt                                                                                                                                                 >  2:610674/1‑150 (MQ=255)
    ggTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATc                                                                                                                                                <  2:344345/150‑1 (MQ=255)
      tCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGt                                                                                                                                              <  2:609906/150‑1 (MQ=255)
      tCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGt                                                                                                                                              <  1:692679/150‑1 (MQ=255)
        aTAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTgg                                                                                                                                            >  2:818321/1‑150 (MQ=255)
              tGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAATTCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttc                                                                                                                                      >  1:441633/1‑150 (MQ=255)
              tGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttc                                                                                                                                      >  2:518636/1‑150 (MQ=255)
              tGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttc                                                                                                                                      >  1:441630/1‑150 (MQ=255)
              tGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttc                                                                                                                                      <  2:690728/150‑1 (MQ=255)
                ttCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttctt                                                                                                                                    >  1:755938/1‑150 (MQ=255)
                ttCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttctt                                                                                                                                    >  1:540929/1‑150 (MQ=255)
                  cGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCg                                                                                                                                  <  1:936533/150‑1 (MQ=255)
                   ggTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCgg                                                                                                                                 >  2:423698/1‑150 (MQ=255)
                    gTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCggg                                                                                                                                >  2:330712/1‑150 (MQ=255)
                    gTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCggg                                                                                                                                <  1:1309196/150‑1 (MQ=255)
                     ttttGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGt                                                                                                                               >  2:64922/1‑150 (MQ=255)
                       ttGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTa                                                                                                                             >  2:1166510/1‑150 (MQ=255)
                          ccTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTAccc                                                                                                                          >  2:1001769/1‑150 (MQ=255)
                          ccTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTAccc                                                                                                                          >  2:1001732/1‑150 (MQ=255)
                           cTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCt                                                                                                                         >  1:757310/1‑150 (MQ=255)
                            tttGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCtt                                                                                                                        >  1:631286/1‑150 (MQ=255)
                             ttGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTa                                                                                                                       <  2:38847/150‑1 (MQ=255)
                             ttGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTa                                                                                                                       <  1:571760/150‑1 (MQ=255)
                               gTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACg                                                                                                                     <  2:393967/150‑1 (MQ=255)
                               gTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACg                                                                                                                     <  2:393914/150‑1 (MQ=255)
                                    caccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCAccc                                                                                                                >  1:130822/1‑149 (MQ=255)
                                    caccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCAAATCGTGGATATTCTTCGGGTTACCCTTACGCACCa                                                                                                                >  1:519684/1‑150 (MQ=255)
                                    caccGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGATTCACGTCGGAGCGCACCAGGTCGTTCAAATCGTGGATATTCTTCGGGTTACCCTTACGCACCa                                                                                                                >  1:194496/1‑150 (MQ=255)
                                         tCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAg                                                                                                           >  2:571341/1‑150 (MQ=255)
                                            gCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGccc                                                                                                        <  2:1179455/150‑1 (MQ=255)
                                            gCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGccc                                                                                                        <  2:148385/150‑1 (MQ=255)
                                             cTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCa                                                                                                       >  2:772580/1‑150 (MQ=255)
                                                    cGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGa                                                                                                <  2:1038804/150‑1 (MQ=255)
                                                    cGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGa                                                                                                <  1:1075602/150‑1 (MQ=255)
                                                    cGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGa                                                                                                <  1:1075605/150‑1 (MQ=255)
                                                         cgcCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGa                                                                                           >  2:1001482/1‑150 (MQ=255)
                                                          gcCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGaa                                                                                          <  1:32258/150‑1 (MQ=255)
                                                          gcCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTGACCCTTACGCACCCGGAAGCCCATGGGGGAGTAGaa                                                                                          >  2:287879/1‑150 (MQ=255)
                                                           ccccATGGCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAAc                                                                                         <  2:1300787/150‑1 (MQ=255)
                                                           ccccATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAAc                                                                                         >  2:1273988/1‑150 (MQ=255)
                                                            cccATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACg                                                                                        >  2:513272/1‑150 (MQ=255)
                                                            ccaaTGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGACGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACg                                                                                        >  2:1108178/4‑150 (MQ=255)
                                                                   cgccAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCt                                                                                 >  2:1045416/1‑150 (MQ=255)
                                                                   cgccAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCt                                                                                 <  1:454422/150‑1 (MQ=255)
                                                                         aTAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTattatt                                                                           <  2:31985/150‑1 (MQ=255)
                                                                          tAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTc                                                                          <  1:512692/150‑1 (MQ=255)
                                                                          tAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTc                                                                          <  1:512699/150‑1 (MQ=255)
                                                                           aGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCg                                                                         >  2:683930/1‑150 (MQ=255)
                                                                              tataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCa                                                                      <  2:648355/150‑1 (MQ=255)
                                                                              tataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCa                                                                      <  2:1203878/150‑1 (MQ=255)
                                                                              tataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGACTAGAACGGCGAGCTATTATTCGGCa                                                                      >  2:336429/1‑150 (MQ=255)
                                                                              tataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAGGGTGGAGTAGAACGGCGAGCTATTATTCGGCa                                                                      >  2:336435/1‑150 (MQ=255)
                                                                               ataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAg                                                                     >  2:649541/1‑150 (MQ=255)
                                                                               ataACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAg                                                                     >  2:649497/1‑150 (MQ=255)
                                                                                 aaCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGGTATTCTTCGGGTTACCCTTACGCACCAGGGAGCCCATGGTGGAGTAGAACGGCGCGCTATTATTCGGCAGgc                                                                   >  1:311933/1‑150 (MQ=255)
                                                                                 aaCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCGTCACGTCGGAGCGCACCAGGTCGTTCCAAACGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgc                                                                   >  1:1016931/1‑150 (MQ=255)
                                                                                  aCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgcg                                                                  <  1:262852/150‑1 (MQ=255)
                                                                                  aCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgcg                                                                  <  1:979956/150‑1 (MQ=255)
                                                                                   cgcgcgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgcgc                                                                 <  2:958254/150‑1 (MQ=255)
                                                                                   cgcgcgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGgcgc                                                                 <  2:309866/150‑1 (MQ=255)
                                                                                     cgcgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGa                                                               <  2:1296709/150‑1 (MQ=255)
                                                                                      gcgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGAc                                                              >  1:725581/1‑150 (MQ=255)
                                                                                      gcgTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGAc                                                              >  1:742244/1‑150 (MQ=255)
                                                                                         ttACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGc                                                           <  1:932350/150‑1 (MQ=255)
                                                                                           aCCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCa                                                         <  2:1079415/150‑1 (MQ=255)
                                                                                           aCCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCa                                                         <  2:1079413/150‑1 (MQ=255)
                                                                                             ccGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACACTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGt                                                       >  2:754940/1‑150 (MQ=255)
                                                                                              cGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTc                                                      >  1:1097196/1‑150 (MQ=255)
                                                                                                 cGTTTTCGGGTTCGGTAAAATCAGCTTCAGGTGGGAGCGCACCAGGTCGTTCCAATCTTGGATATTCTTGGGGTTACCCTTACGCACCAGGAAGCCCTTGGTGTTGTAGAAGGGCGAGCTATTATTCGGCTGGCGCGACTGCCAGTCGGc                                                   <  1:818321/150‑1 (MQ=255)
                                                                                                     ttcgggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCggg                                               >  2:800422/1‑150 (MQ=255)
                                                                                                     ttcgggttcgggAAAATCAGCTTCACGTCGCAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGAAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCggg                                               >  2:800522/1‑150 (MQ=255)
                                                                                                      tcgggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGa                                              >  2:200860/1‑150 (MQ=255)
                                                                                                      tcgggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGa                                              >  1:584354/1‑150 (MQ=255)
                                                                                                        gggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCAAATCGTGTATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGTGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATc                                            <  1:915124/150‑1 (MQ=255)
                                                                                                         ggttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCa                                           >  2:520708/1‑150 (MQ=255)
                                                                                                           ttcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGc                                         <  1:906430/150‑1 (MQ=255)
                                                                                                            tcgggAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCt                                        >  2:614257/1‑150 (MQ=255)
                                                                                                             cggcAAAATCAGCTTCACGTCGGCGCGCACCCGGTCGTTCAAATCGTGTATATTCTTGGGTTTACCCTTCCGCACCCGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCtt                                       <  1:359513/150‑1 (MQ=255)
                                                                                                                    aTCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTa                                <  2:1097196/150‑1 (MQ=255)
                                                                                                                     tCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTAt                               <  2:235948/150‑1 (MQ=255)
                                                                                                                     tCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTAt                               <  2:1394548/150‑1 (MQ=255)
                                                                                                                        gCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGt                            <  1:29774/150‑1 (MQ=255)
                                                                                                                              cGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAGGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGa                      <  1:1095156/150‑1 (MQ=255)
                                                                                                                               gTCGGAGCGCACCAGGTCGTTCCAATCGTGGTTATTCTTAGGGTTAACCTTACGCACCAGGAATCCCATGGTGGAGTAGAACGGCGAGCTATTATTAGGCAGGCGCGACTGCCAGTCGGCCGGGATAAGCTTGCCTTTATAGTTCAGGAt                     >  1:251704/1‑150 (MQ=255)
                                                                                                                                    aGCGCACCAGGTCTTTCCACTCGTGGTTATTCGTCGGGTTACCCGTCCGCACCAGTATGCCCATGGTGGAGTTGAACGGCGAGCTATTATTCGGCTGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTa                <  1:1345900/150‑1 (MQ=255)
                                                                                                                                       gcACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGt             <  2:550846/150‑1 (MQ=255)
                                                                                                                                                  ttCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTCGGTCACCTgg  <  2:134970/150‑1 (MQ=255)
                                                                                                                                                  ttCAAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTCGGTCACCTgg  <  2:519684/150‑1 (MQ=255)
                                                                                                                                                  ttCAAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTCGGTCACCTgg  <  2:194496/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGCTATTATTCGGCAGGCGCGACTGCCAGTCGGCCGGGATCAGCTTGCCTTTATCGTGCAGGATTTGTACGTCGGTCACCTGG  >  NC_000913/2542970‑2543265

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