Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,483,245 G→C 7.7% S15T (AGT→ACT) ‡ yheT → UPF0017 family putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,483,2450GC 7.7% 221.6 / 11.1 78S15T (AGT→ACT) ‡yheTUPF0017 family putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base C (3/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

GCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGAGTGAAA  >  NC_000913/3483097‑3483394
                                                                                                                                                    |                                                                                                                                                     
gcAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGt                                                                                                                                                      >  2:216037/1‑150 (MQ=255)
       cAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAAt                                                                                                                                               <  2:757476/150‑1 (MQ=255)
       cAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAt                                                                                                                                               <  2:1137484/150‑1 (MQ=255)
         aGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTc                                                                                                                                             >  2:1141095/1‑150 (MQ=255)
         aGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACAGATGGCGCAGATACCGCCGCCAGCTGCCAATGAATTCAGCCGCAGTGCTGAAATc                                                                                                                                             >  2:98708/1‑150 (MQ=255)
            cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAtt                                                                                                                                          <  1:537177/150‑2 (MQ=255)
             aGCGCCAAATCCGGCCTGTAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGTAGCAGCTTGAGCAGATGCTGCTGGAAGGACAAAGCAACTGATGGCGCAGATAACGACGACAGATGCCAATGAATTCAGCACTGCTGCTGAATTCAttg                                                                                                                                         <  2:860155/150‑3 (MQ=255)
              gcgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATccc                                                                                                                                        <  1:29114/150‑1 (MQ=255)
               cgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATcccc                                                                                                                                       >  2:685520/1‑150 (MQ=255)
               cgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATcccc                                                                                                                                       >  2:345928/1‑150 (MQ=255)
                gcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCa                                                                                                                                      <  1:732690/150‑1 (MQ=255)
                 ccAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCAt                                                                                                                                     >  2:417379/1‑150 (MQ=255)
                  cAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTTGACGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTCATGGCGCTGAAAACGACGCCCGATGCCAATGAATTCCGCAGCAGTGCTGAATTCATCCCCATg                                                                                                                                    >  2:1057180/1‑150 (MQ=255)
                       ccGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGCTGCTGCTGGAAGGCCAAAGCAACTGATGGCGAAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCgg                                                                                                                               >  1:690807/1‑150 (MQ=255)
                             tGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAg                                                                                                                         >  1:742971/1‑150 (MQ=255)
                             tGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAg                                                                                                                         <  2:378657/150‑1 (MQ=255)
                                 agagTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAATGCAACTGATGGCGCAGATAACGACGCCCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                                     >  2:227655/1‑150 (MQ=255)
                                  gagTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAtt                                                                                                                    >  1:406504/1‑150 (MQ=255)
                                      gCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                                >  1:1015834/1‑150 (MQ=255)
                                      gCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                                <  1:539984/150‑1 (MQ=255)
                                      gCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGAAGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                                <  1:845922/150‑1 (MQ=255)
                                       cGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGAAGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCAt                                                                                                               >  2:673293/1‑150 (MQ=255)
                                          actggcatggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGACAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTg                                                                                                            <  2:790141/148‑1 (MQ=255)
                                          aaTGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGATGGAAGGCCAAAGCAACGGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGATTTCATCCCCCTGCGCGGCTTTAGCAATTGTCATCTg                                                                                                            >  1:1177537/1‑150 (MQ=255)
                                              gcatggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCaaa                                                                                                        >  2:1176194/1‑150 (MQ=255)
                                               catggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAc                                                                                                       <  2:848621/150‑1 (MQ=255)
                                                    cTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCAtgc                                                                                                  >  1:682266/1‑150 (MQ=255)
                                                     tGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCAtgct                                                                                                 <  1:1154281/150‑1 (MQ=255)
                                                         atcccAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAAGTGATGGCGCAGATAACGACGACCGATGCCAATGATTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcg                                                                                             <  1:675987/148‑1 (MQ=255)
                                                           cccAGGAGCAGCTTGAGAAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCGGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcgcg                                                                                           <  1:605130/150‑1 (MQ=255)
                                                            ccAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGCTGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGGCGCGt                                                                                          >  1:188112/1‑150 (MQ=255)
                                                            ccAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGt                                                                                          <  1:474686/150‑1 (MQ=255)
                                                                  gcagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                                    <  2:336539/150‑1 (MQ=255)
                                                                  gcagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                                    <  1:798138/150‑1 (MQ=255)
                                                                   cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                                   >  1:752251/1‑150 (MQ=255)
                                                                    agcTTGAGCAGATGCTGCTGGATGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCCGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcg                                                                                  >  1:184483/1‑150 (MQ=255)
                                                                     gcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgt                                                                                 >  1:268621/1‑150 (MQ=255)
                                                                     gcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgt                                                                                 >  2:654857/1‑150 (MQ=255)
                                                                      cTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtc                                                                                >  2:893193/1‑150 (MQ=255)
                                                                        tGAGCAGATGCTGCTGGAAGGCAAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGACATCTGCAAACCATGCAGCCGCGTCTGTTTcgtcgt                                                                              >  1:731626/1‑150 (MQ=255)
                                                                           gCAGTTGCTGCGGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAg                                                                           <  1:95181/150‑1 (MQ=255)
                                                                           acAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAAATCAGCAGCAGTGCGGATATCATACCCATGCGCGGCTATAGAAAATGTCATCTGCAAACCATGCTGCCGGGTCTGTTAGGTCGTCAg                                                                           >  2:702462/2‑150 (MQ=255)
                                                                              gATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCAGGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTg                                                                        <  1:763440/150‑1 (MQ=255)
                                                                               aTGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGa                                                                       <  2:712399/150‑1 (MQ=255)
                                                                                    tctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAGCATGGTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTc                                                                  <  2:534233/149‑1 (MQ=255)
                                                                                    gctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGGCTGTTTCGTCGTCAGGTGAAATTc                                                                  <  2:316876/150‑1 (MQ=255)
                                                                                       ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGCCCGATGCCAATGAATTCAGCAGCAGTGCTGAAGTCATCCGCATGCGCGGCTTTAGCTATTGTGATGGGCACACCATGCTGCCGCGTCTGTTTCGTCGTGCGGTGAAATTCAcc                                                               >  2:822280/1‑150 (MQ=255)
                                                                                       ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAcc                                                               <  1:756616/150‑1 (MQ=255)
                                                                                        gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAccc                                                              <  2:1051884/150‑1 (MQ=255)
                                                                                         aaGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAcccc                                                             <  2:468859/150‑1 (MQ=255)
                                                                                           ggCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGt                                                           >  1:1129542/1‑150 (MQ=255)
                                                                                            gCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTa                                                          >  2:874935/1‑150 (MQ=255)
                                                                                            gCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTa                                                          >  2:301679/1‑150 (MQ=255)
                                                                                              cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTGTCGTCGTCAGGTGAAATTCACCCCGTACt                                                        >  2:801258/1‑150 (MQ=255)
                                                                                                 aGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGc                                                     <  2:394106/150‑1 (MQ=255)
                                                                                                 aGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATACCCATGCGCGGCTTTAGCAATTGTCATCTGCACACAATGCTGCCGCGACTGTTTCGTCGTCAGGTGAAATTCACCCATTACTGGc                                                     >  2:738390/1‑150 (MQ=255)
                                                                                                     aCTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCg                                                 <  1:448762/150‑1 (MQ=255)
                                                                                                          tGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTgg                                            <  2:28304/150‑1 (MQ=255)
                                                                                                            gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAg                                          >  2:159797/1‑150 (MQ=255)
                                                                                                             cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCCCCCCGTACTGGCAGCGGCTGGAGt                                         >  2:716923/1‑150 (MQ=255)
                                                                                                             cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGAGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATCCTGCCCCGGCGGTTTCGGGGTCAGGTGAAATTACCCCCGGAGAGGCAGCGGCTGgtgt                                         >  1:1263474/1‑150 (MQ=255)
                                                                                                                aGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGc                                      <  2:520864/150‑1 (MQ=255)
                                                                                                                 caatacgacgACCGAAGCCAATGAATTCCGCCGCAGTGCTGAATTCATCCCCATGCGCGGCTGTAGCAATTGTCATCTGCAACCCATGCTGCCGCGTCTGCTTCGTCGTCAGGTGAAATTCACCCCGGACTGGCAGCGGCTGGAGTTGcc                                     <  1:862027/146‑1 (MQ=255)
                                                                                                                   tAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCg                                   <  2:174268/150‑1 (MQ=255)
                                                                                                                         cgacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCg                             <  1:267678/150‑1 (MQ=255)
                                                                                                                                tGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGt                      <  1:646399/150‑1 (MQ=255)
                                                                                                                                 gCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTc                     <  1:866475/150‑1 (MQ=255)
                                                                                                                                  ccAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGGTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGAATTTGTCg                    >  2:301277/1‑150 (MQ=255)
                                                                                                                                  ccAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg                    >  1:10327/1‑150 (MQ=255)
                                                                                                                                   ggaTGAATTCAGCACTGCTGCTGAATTCATTCCCATGCGCGGCTTTAGCTATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGa                   >  1:688181/3‑150 (MQ=255)
                                                                                                                                   cAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACACCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGa                   >  2:886552/1‑150 (MQ=255)
                                                                                                                                    gaTGAATTCAGCTCTGCTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt                  >  1:626551/2‑150 (MQ=255)
                                                                                                                                    aaTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt                  <  1:801258/150‑1 (MQ=255)
                                                                                                                                      tGAATTCAGCACTGCTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCt                >  1:1165737/1‑150 (MQ=255)
                                                                                                                                        aaTTCAGCAGCAGTGCTGAATTCATCCCCAAGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTg              >  2:646780/1‑150 (MQ=255)
                                                                                                                                         aTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGc             <  1:1240315/150‑1 (MQ=255)
                                                                                                                                          ttCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGTCTGGAGTTGCCCGACGGCGATTTTGTGGATCTTGCg            >  2:346190/1‑150 (MQ=255)
                                                                                                                                          ttCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGCTGAAATTCACCCCGTACTGGCTGCGGCGGGAGTTGCCCGACGGCCATATTGTCGATCATGCg            >  2:606602/1‑150 (MQ=255)
                                                                                                                                            cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTg          <  1:654857/150‑1 (MQ=255)
                                                                                                                                            cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTg          <  1:1230791/150‑1 (MQ=255)
                                                                                                                                            aagcatcagTGCTGAATTCATCCCCATGCGCGTCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGGGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCTATCTTGCGTg          <  1:622946/149‑1 (MQ=255)
                                                                                                                                             agcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTgg         <  1:685520/150‑1 (MQ=255)
                                                                                                                                               cagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGAg       >  2:1289558/1‑150 (MQ=255)
                                                                                                                                                 gcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGAGTg     <  1:93092/150‑1 (MQ=255)
                                                                                                                                                    gTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCGGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGAGTGaaa  >  1:228253/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTGGAGTGAAA  >  NC_000913/3483097‑3483394

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: