Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,958,683 A→T 5.1% V195D (GTC→GAC)  cutC ← copper homeostasis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,958,6830AT 5.1% 303.5 / 3.6 96V195D (GTC→GAC) cutCcopper homeostasis protein
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base T (2/3);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCTGATGTCAG  >  NC_000913/1958540‑1958831
                                                                                                                                               |                                                                                                                                                    
cAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTc                                                                                                                                                  <  1:163880/148‑1 (MQ=255)
cAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTc                                                                                                                                                  <  8:474777/148‑1 (MQ=255)
cAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTc                                                                                                                                                  <  6:113079/148‑1 (MQ=255)
 aaTGATTCCTTTGATTTCAGCAAGCGCCGCCCCGTCTACGATATTACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGTTTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTcc                                                                                                                                                 <  3:217547/148‑1 (MQ=255)
 aaTGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTc                                                                                                                                                  >  1:277868/1‑147 (MQ=255)
   tGATTCCTTTCATTTCAGCAACCGCCGCCCCGCCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAg                                                                                                                                               >  1:470104/1‑148 (MQ=255)
         cTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATCCTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCg                                                                                                                                        <  3:39358/149‑1 (MQ=255)
          tttcatttcaGCAACCGCCGCCCCGTCTACGATATAACGCGCATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCgg                                                                                                                                       >  1:331268/1‑149 (MQ=255)
          tttcatttcaGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCgg                                                                                                                                       >  8:128761/1‑149 (MQ=255)
             catttaaGCACCCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCAt                                                                                                                                    >  6:203010/1‑149 (MQ=255)
              atttcaGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATc                                                                                                                                   >  5:55631/1‑149 (MQ=255)
               tttcaGCAACCGCCGACCCCTCTACGATATAACGAGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATTCCCCATCGGTGAGGCTTGCCACGCTCCCGCGCAGCTATGCACTTCCAGCACTCCGCCATc                                                                                                                                   >  6:302509/1‑148 (MQ=255)
                 tcaGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCgg                                                                                                                                 <  3:216900/148‑2 (MQ=255)
                  caGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAg                                                                                                                                >  4:197110/1‑148 (MQ=255)
                   aGCAACCGCCGCCCCGTCTACGATATACCGCGAATACTCGTCCGCGTGTTCATCTGAAGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAgg                                                                                                                               <  5:302509/148‑1 (MQ=255)
                    gCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAgg                                                                                                                               >  3:337453/1‑147 (MQ=255)
                    gCAACCGCCGCCCCGCCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGaa                                                                                                                             <  2:470104/149‑1 (MQ=255)
                        ccgccgCCCCGTCTTCGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAgtgg                                                                                                                         <  3:430777/149‑1 (MQ=255)
                             gccCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAg                                                                                                                    >  1:331523/1‑149 (MQ=255)
                                ccGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGtt                                                                                                                 <  7:115293/149‑1 (MQ=255)
                                 cGTCTACGATATATCTCGAATACTGGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTTCGTTAACGGAGGGGTGAGGCTTGAGAGGGTCCCGCGGAGCTATGGTCTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAg                                                                                                                    <  1:162754/145‑1 (MQ=255)
                                 cGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGttt                                                                                                                >  3:43063/1‑149 (MQ=255)
                                 cGTCAACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAAACCTTGATAACGATAACGCATCCGTGAGACATCCCACGCTCCCGCGGAGCTATGGACATCCAGCACACCGGAATAGAGGAAGACGAGCAggag                                                                                                                 >  2:162754/1‑146 (MQ=255)
                                   tCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTc                                                                                                              <  7:128761/149‑1 (MQ=255)
                                     tACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAg                                                                                                                                >  6:22516/1‑129 (MQ=255)
                                     tACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAg                                                                                                                                <  5:22516/129‑1 (MQ=255)
                                     tACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTg                                                                                                            >  4:194267/1‑149 (MQ=255)
                                        gATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCAc                                                                                                         >  3:217218/1‑149 (MQ=255)
                                          tataACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAg                                                                                                                                >  4:234257/1‑124 (MQ=255)
                                          tataACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAg                                                                                                                                <  3:234257/124‑1 (MQ=255)
                                          tataACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACg                                                                                                        >  6:458762/1‑148 (MQ=255)
                                            taACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGAc                                                                                                     <  6:440844/149‑1 (MQ=255)
                                                          gTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGcc                                                                                        <  5:458762/148‑1 (MQ=255)
                                                                     cATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTgg                                                                             <  8:191242/148‑1 (MQ=255)
                                                                          gatgaCATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCAt                                                                        >  2:257122/1‑148 (MQ=255)
                                                                           atgaCATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGc                                                                          <  7:260808/145‑1 (MQ=255)
                                                                           atgaCATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCa                                                                         <  6:55631/146‑1 (MQ=255)
                                                                           atgaCATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCAAGCACTCCGGCATCGAGGAAGTGCTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGCGCATCa                                                                      >  1:375119/1‑149 (MQ=255)
                                                                             gaCATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCat                                                                     >  1:247460/1‑147 (MQ=255)
                                                                                aTGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATg                                                                 >  6:121153/1‑149 (MQ=255)
                                                                                 tGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAACACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                <  7:74274/149‑1 (MQ=255)
                                                                                  ggACAATCCTAGATTACGATAACGCATCGGTGAGGCTTGACACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATggg                                                               >  8:70154/1‑149 (MQ=255)
                                                                                   gACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATggg                                                               >  7:360858/1‑148 (MQ=255)
                                                                                   gACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGc                                                              <  8:360858/149‑1 (MQ=255)
                                                                                   gACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGc                                                              <  3:197110/149‑1 (MQ=255)
                                                                                    aCAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCATCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGc                                                              >  1:238544/1‑148 (MQ=255)
                                                                                      aaTCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                >  8:128412/1‑144 (MQ=255)
                                                                                      aaTCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                <  7:128412/144‑1 (MQ=255)
                                                                                      aaTCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCaa                                                            >  8:67778/1‑148 (MQ=255)
                                                                                        tCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCa                                                                                                                                                <  6:188086/62‑1 (MQ=255)
                                                                                        tCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCa                                                                                                                                                >  5:188086/1‑62 (MQ=255)
                                                                                        tCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAAt                                                           >  7:381826/1‑147 (MQ=255)
                                                                                           ttGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGtt                                                      >  4:261071/1‑149 (MQ=255)
                                                                                             gATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGtt                                                      <  8:505021/147‑1 (MQ=255)
                                                                                             gATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTc                                                     <  3:194267/148‑1 (MQ=255)
                                                                                              aTTACGATAACGCAACGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCa                                                   >  7:485144/1‑149 (MQ=255)
                                                                                                   gATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAAtt                                              >  2:395511/1‑149 (MQ=255)
                                                                                                      aaCGCATCGGGGAGGCTTGCCACGCTCCCGCGGAGCTATGGACATCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAAttttt                                           <  1:409834/149‑1 (MQ=255)
                                                                                                       aCGCTTCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCATCACTCCGGCATCGAGGAAGTGGTGCAGGTGTTCTGCAAGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAAttttt                                           <  2:375119/148‑1 (MQ=255)
                                                                                                       aCGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAAttttt                                           <  4:43063/148‑1 (MQ=255)
                                                                                                         gCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAAttttt                                           <  5:415855/146‑1 (MQ=255)
                                                                                                         gCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAAttttt                                           >  6:415855/1‑146 (MQ=255)
                                                                                                         gCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGa                                         >  2:279112/1‑148 (MQ=255)
                                                                                                         ccATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTACCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCTTAATGATTGGAGCATCAAGATGCCCAATAAGTTACATAATATATGa                                         >  6:8956/2‑148 (MQ=255)
                                                                                                            tCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATaaa                                     <  1:257122/149‑1 (MQ=255)
                                                                                                              ggTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAAc                                    <  8:381826/148‑1 (MQ=255)
                                                                                                                  aGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTg                                <  4:217218/148‑1 (MQ=255)
                                                                                                                    gCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAg                             >  3:141487/1‑149 (MQ=255)
                                                                                                                     cTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGt                            <  2:277868/149‑1 (MQ=255)
                                                                                                                      ttGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                <  1:408678/112‑1 (MQ=255)
                                                                                                                      ttGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                >  2:408678/1‑112 (MQ=255)
                                                                                                                      ttGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGt                            >  8:26496/1‑148 (MQ=255)
                                                                                                                        gCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGt                            <  7:26496/146‑1 (MQ=255)
                                                                                                                        gCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGc                          <  5:242790/148‑1 (MQ=255)
                                                                                                                        gCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCg                         >  3:114275/1‑149 (MQ=255)
                                                                                                                          cACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGc                          >  6:242790/1‑146 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGTCTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATg                                                                 >  8:424099/1‑102 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGTCTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATg                                                                 >  8:416974/1‑102 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGTCTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATg                                                                 <  7:416974/102‑1 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGTCTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATg                                                                 <  7:424099/102‑1 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                >  2:461955/1‑103 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATgg                                                                <  1:461955/103‑1 (MQ=255)
                                                                                                                               tCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGAtt                  >  6:33538/1‑149 (MQ=255)
                                                                                                                                 ccGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGAtttt                >  5:243590/1‑149 (MQ=255)
                                                                                                                                  cgcgGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGAttttt               <  1:279112/149‑1 (MQ=255)
                                                                                                                                  cgcgGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTACTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGAttttt               <  5:121153/149‑1 (MQ=255)
                                                                                                                                   gcgGAGGTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTGCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTg              <  6:112214/149‑1 (MQ=255)
                                                                                                                                       aGCTATGGACTTCCAGCACTCCGGCATCGCGGAAGGGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCt          <  2:331268/149‑1 (MQ=255)
                                                                                                                                        gCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCt          >  3:200815/1‑148 (MQ=255)
                                                                                                                                           aTGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCTGATg      >  2:94656/1‑149 (MQ=255)
                                                                                                                                            tGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCTGATGt     <  4:141487/149‑1 (MQ=255)
                                                                                                                                             ggACTTCCAGCACTCCGGCATCGAGGAAGTGGTGc                                                                                                                      <  4:41608/35‑1 (MQ=255)
                                                                                                                                             ggACTTCCAGCACTCCGGCATCGAGGAAGTGGTGc                                                                                                                      >  3:41608/1‑35 (MQ=255)
                                                                                                                                             ggACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTgg                                                                             >  1:354556/1‑76 (MQ=255)
                                                                                                                                             ggACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTgg                                                                             <  2:354556/76‑1 (MQ=255)
                                                                                                                                              gACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCTGATGTCa   <  7:440900/149‑1 (MQ=255)
                                                                                                                                               aCTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCTGATGTCAg  <  7:67778/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGACAATCCTTGATTACGATAACGCATCGGTGAGGCTTGCCACGCTCCCGCGGAGCTATGGACTTCCAGCACTCCGGCATCGAGGAAGTGGTGCAGGTTTTCTGCACGGACTCCTGCTCCGGCCATAATGATTGGAGCATCACGATGGGCAATAAGTTCCATAATTTTTGATAAACCTTGCAGTGCGTCTGATTTTTGCCCTGATGTCAG  >  NC_000913/1958540‑1958831

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