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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,483,2450GC 7.3% 305.7 / 10.4 93S15T (AGT→ACT) ‡yheTputative hydrolase YheT
Reads supporting (aligned to +/- strand):  ref base G (41/45);  new base C (3/4);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

TGCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTG  >  NC_000913/3483096‑3483386
                                                                                                                                                     |                                                                                                                                             
tgcAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTcagcagcag                                                                                                                                               >  2:1973/1‑150 (MQ=255)
   aaCAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGc                                                                                                                                            <  2:321046/150‑1 (MQ=255)
   aaCAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGc                                                                                                                                            <  2:170525/150‑1 (MQ=255)
    aCAGCAAGCCAGTGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCt                                                                                                                                           >  2:16540/1‑150 (MQ=255)
      aGCAAGCCACCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGa                                                                                                                                         >  2:184807/1‑150 (MQ=255)
       gCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGaa                                                                                                                                        >  2:65956/1‑150 (MQ=255)
        cAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAAt                                                                                                                                       <  1:966818/150‑1 (MQ=255)
        cAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAt                                                                                                                                       <  1:959444/150‑1 (MQ=255)
         aaGCCAGCGCCAAATCCGGCGTGGAAGAGTGCGAAATGGCTTGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGTCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAtt                                                                                                                                      <  2:1074808/150‑1 (MQ=255)
         aaGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAAtt                                                                                                                                      <  2:1355278/150‑1 (MQ=255)
         aaGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAtt                                                                                                                                      <  2:1289918/150‑1 (MQ=255)
          aGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCGGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATGCAGCAGCAGTGCTGAATTc                                                                                                                                     <  2:1456904/150‑1 (MQ=255)
           gCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCa                                                                                                                                    >  1:949619/1‑150 (MQ=255)
             cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGCCGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                  >  1:289726/1‑150 (MQ=255)
             cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                  >  2:1417101/1‑150 (MQ=255)
             cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                  <  1:636971/150‑1 (MQ=255)
             cAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                  >  1:1282261/1‑150 (MQ=255)
               gcgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTCAttgg                                                                                                                                <  1:988537/150‑4 (MQ=255)
                cgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATcccc                                                                                                                               <  2:206818/150‑1 (MQ=255)
                 gcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCa                                                                                                                              >  2:1517983/1‑150 (MQ=255)
                 gcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCa                                                                                                                              >  1:540854/1‑150 (MQ=255)
                   cAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATg                                                                                                                            <  1:165379/150‑1 (MQ=255)
                   cAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATg                                                                                                                            <  1:804187/150‑1 (MQ=255)
                     aaTCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATgcg                                                                                                                          >  1:915315/1‑150 (MQ=255)
                         cGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGc                                                                                                                      >  2:1163112/1‑150 (MQ=255)
                            ccTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCttt                                                                                                                   <  1:1381553/150‑1 (MQ=255)
                             cTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTa                                                                                                                  >  2:129339/1‑150 (MQ=255)
                               ggAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGc                                                                                                                <  1:439643/150‑1 (MQ=255)
                               ggAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGc                                                                                                                <  2:1052043/150‑1 (MQ=255)
                                gAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTTCTGAATTCATCCCCATGCGCGGCTTTAGCa                                                                                                               <  2:52000/150‑1 (MQ=255)
                                    agTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTg                                                                                                           <  1:1232321/150‑1 (MQ=255)
                                     gTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGt                                                                                                          >  2:1326153/1‑150 (MQ=255)
                                            atggcatggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGc                                                                                                   <  2:1424079/150‑1 (MQ=255)
                                                 atggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAcc                                                                                              <  2:215509/150‑1 (MQ=255)
                                                     cTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCAtgc                                                                                          >  2:761275/1‑150 (MQ=255)
                                                       ggAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCAtgctg                                                                                        >  2:12481/1‑150 (MQ=255)
                                                        gAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCAtgctgc                                                                                       <  2:385277/150‑1 (MQ=255)
                                                          aGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcg                                                                                     <  1:138441/150‑1 (MQ=255)
                                                            cccAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcgcg                                                                                   <  1:263097/150‑1 (MQ=255)
                                                                   tcagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                            <  1:1163112/149‑1 (MQ=255)
                                                                      gcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgt                                                                         >  2:1252762/1‑150 (MQ=255)
                                                                      gcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgt                                                                         >  1:466637/1‑150 (MQ=255)
                                                                        ttGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtcg                                                                       <  2:465537/150‑1 (MQ=255)
                                                                        ttGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtcg                                                                       <  1:127456/150‑1 (MQ=255)
                                                                         tGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTcgtcgt                                                                      >  2:1423355/1‑150 (MQ=255)
                                                                           aGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCa                                                                    >  2:1131940/1‑150 (MQ=255)
                                                                                 tgctgcTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaa                                                              >  1:426619/1‑150 (MQ=255)
                                                                                  gctgctGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaaa                                                             >  2:1039267/1‑150 (MQ=255)
                                                                                   ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAAt                                                            >  1:834188/1‑150 (MQ=255)
                                                                                   ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAAt                                                            <  1:496516/150‑1 (MQ=255)
                                                                                   ctgctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAAt                                                            <  1:1417101/150‑1 (MQ=255)
                                                                                     gctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTc                                                          >  2:1216413/1‑150 (MQ=255)
                                                                                         gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAccc                                                      <  2:949619/150‑1 (MQ=255)
                                                                                           aGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                           >  2:839849/1‑127 (MQ=255)
                                                                                           aGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtc                                                                           <  1:839849/127‑1 (MQ=255)
                                                                                             gCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTa                                                  >  1:914697/1‑150 (MQ=255)
                                                                                             gCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTa                                                  >  1:131467/1‑150 (MQ=255)
                                                                                               cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACt                                                <  1:184807/150‑1 (MQ=255)
                                                                                                aaaGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTg                                               <  1:65956/150‑1 (MQ=255)
                                                                                                     aaCTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGc                                          >  1:608475/1‑150 (MQ=255)
                                                                                                        tGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGc                                       <  1:1039267/150‑1 (MQ=255)
                                                                                                           tGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTgg                                    <  2:1412009/150‑1 (MQ=255)
                                                                                                            ggCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGTAGCGGCTGGa                                   >  2:243180/1‑150 (MQ=255)
                                                                                                            ggCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGa                                   >  1:662807/1‑150 (MQ=255)
                                                                                                            ggCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGa                                   <  1:1973/150‑1 (MQ=255)
                                                                                                             gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAg                                  >  2:1000458/1‑150 (MQ=255)
                                                                                                             gcgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAg                                  <  2:927940/150‑1 (MQ=255)
                                                                                                              cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGGCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGt                                 >  2:665393/1‑150 (MQ=255)
                                                                                                                  gATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGcc                             <  2:510611/150‑1 (MQ=255)
                                                                                                                    tAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCg                           <  1:1517983/150‑1 (MQ=255)
                                                                                                                       cgacgaCCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACg                        >  1:344436/1‑150 (MQ=255)
                                                                                                                       cgacgaCCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACg                        >  2:478105/1‑150 (MQ=255)
                                                                                                                       cgacgaCCCATGCCAATGAATTCAGCATCAGTGCTCAATTCATCCCCATGCGCGGCTTTAGCAATTTTTTTCTGCAAACCTTCGTGCCTCGTGTGTTTCTTCGTGAGTTTATATTTACCGCTTACTGGCTGAGTCTGTAGTTgtcatacg                        <  2:517326/150‑8 (MQ=255)
                                                                                                                        gacgacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACgg                       <  1:1131940/150‑1 (MQ=255)
                                                                                                                         acgacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGc                      <  1:1505013/150‑1 (MQ=255)
                                                                                                                         acgacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGc                      >  1:260687/1‑150 (MQ=255)
                                                                                                                          cgacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCg                     >  1:427059/1‑150 (MQ=255)
                                                                                                                           gacCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGa                    <  2:834188/150‑1 (MQ=255)
                                                                                                                            acCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGAt                   >  2:362541/1‑150 (MQ=255)
                                                                                                                             ccGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGAtt                  <  1:1216413/150‑1 (MQ=255)
                                                                                                                               gATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGAtttt                >  1:484650/1‑150 (MQ=255)
                                                                                                                                aTGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTg               <  2:484650/150‑1 (MQ=255)
                                                                                                                                  gCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTc             <  1:129339/150‑1 (MQ=255)
                                                                                                                                   ccAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg            <  2:426619/150‑1 (MQ=255)
                                                                                                                                   ccAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCg            <  2:1282261/150‑1 (MQ=255)
                                                                                                                                    ggaTGAATTCAGCACTGCTGCTGCATTCATCCCCATGCGCGGCGTTCGCAATTGTCATCTGCACACCATGCTGCCGCGGCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACGGGCAGCGGCTGGAGTTGCCCGCCGGGGATTTTGTCgg           >  2:511599/3‑149 (MQ=255)
                                                                                                                                     aaTGAATTCAGCAGAAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt          >  2:67775/1‑150 (MQ=255)
                                                                                                                                       tGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCt        <  2:408139/150‑1 (MQ=255)
                                                                                                                                       tGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGACTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCt        >  1:1049344/1‑150 (MQ=255)
                                                                                                                                        gAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCtt       >  2:1002282/1‑150 (MQ=255)
                                                                                                                                        gAATTCAGCACTGCTGCTGAATTCATTCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCAGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGACCGACGGCGATTTTGTCGATCtt       >  1:877063/1‑150 (MQ=255)
                                                                                                                                         aaTTCAGCACTGCTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTg      >  1:1329255/1‑150 (MQ=255)
                                                                                                                                          aTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGc     <  2:427059/150‑1 (MQ=255)
                                                                                                                                          aTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGc     >  1:975481/1‑150 (MQ=255)
                                                                                                                                             cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTg  >  2:639736/1‑150 (MQ=255)
                                                                                                                                             cagcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTg  <  2:260687/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TGCAACAGCAAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTG  >  NC_000913/3483096‑3483386

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