Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,462,345 C→G 5.9% V106L (GTC→CTC)  nrdD ← anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,462,3450CG 5.9% 298.4 / 6.5 82V106L (GTC→CTC) nrdDanaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (28/49);  new base G (3/2);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATGTGCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCATTGATATCC  >  NC_000913/4462197‑4462493
                                                                                                                                                    |                                                                                                                                                    
tCATGTGCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGTTACGATCCCGGCCAGCAGGTCGCACTGGGCTGGAATCACCTTGCTGACTTTGTTGGCGTTTTCGTGGAGTAACGAGGCGTTGGTCTGCTCgac                                                                                                                                                      >  2:220844/1‑149 (MQ=255)
tCATGTGCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCgac                                                                                                                                                      >  4:276365/1‑149 (MQ=255)
tCATGTGCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCgac                                                                                                                                                      >  3:613615/1‑149 (MQ=255)
    gtgCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCgaccaga                                                                                                                                                  <  1:216418/149‑1 (MQ=255)
     tgCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCgaccagac                                                                                                                                                 <  4:102501/149‑1 (MQ=255)
         tGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGGTCgagcagaccaac                                                                                                                                             <  2:104632/149‑3 (MQ=255)
              cacGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGctct                                                                                                                                        >  3:615351/1‑149 (MQ=255)
              cacGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGAtct                                                                                                                                        >  1:518636/1‑149 (MQ=255)
              cacGTCACGCGGCAGCAGGTGCCGACGTGCATAGTGTTTAGCCACGATCCCTGCCAGCAAGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTGTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGCAtct                                                                                                                                        >  4:386509/1‑149 (MQ=255)
               acGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGAtctc                                                                                                                                       <  1:383389/149‑1 (MQ=255)
                 gTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCt                                                                                                                                     >  3:171319/1‑149 (MQ=255)
                  tCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTg                                                                                                                                    <  3:494050/149‑1 (MQ=255)
                  tCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTg                                                                                                                                    <  1:443529/149‑1 (MQ=255)
                  tCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTg                                                                                                                                    <  4:258881/149‑1 (MQ=255)
                    aCGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGt                                                                                                                                  >  4:315906/1‑149 (MQ=255)
                      gcgGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTc                                                                                                                                <  4:123006/149‑1 (MQ=255)
                      gcgGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTc                                                                                                                                <  3:47794/149‑1 (MQ=255)
                        ggCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAg                                                                                                                              <  1:436971/149‑1 (MQ=255)
                          cagcagGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGc                                                                                                                            <  4:165428/149‑1 (MQ=255)
                                  gCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGt                                                                                                                    <  1:416215/149‑1 (MQ=255)
                                   cTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGtt                                                                                                                   <  2:135971/149‑1 (MQ=255)
                                    tGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGttt                                                                                                                  >  1:451787/1‑149 (MQ=255)
                                    tGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGttt                                                                                                                  <  2:199029/149‑1 (MQ=255)
                                     gACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGGTCgagcagaccaac                                                                                                                                             <  1:321077/121‑3 (MQ=255)
                                     gACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGGTCgagcagaccaac                                                                                                                                             >  2:321077/1‑119 (MQ=255)
                                      aCGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGttttt                                                                                                                <  1:364290/149‑1 (MQ=255)
                                        gTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCa                                                                                                              >  2:183022/1‑149 (MQ=255)
                                        gTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCa                                                                                                              <  2:166384/149‑1 (MQ=255)
                                                tgtTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAAt                                                                                                      >  3:12512/1‑149 (MQ=255)
                                                    tAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTcg                                                                                                  <  1:169157/149‑1 (MQ=255)
                                                    tAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTcg                                                                                                  >  4:598593/1‑149 (MQ=255)
                                                        cACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGAt                                                                                              >  3:541735/1‑149 (MQ=255)
                                                        cACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGAt                                                                                              <  4:613615/149‑1 (MQ=255)
                                                            aTCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTg                                                                                          >  3:288553/1‑149 (MQ=255)
                                                                     agcagGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTc                                                                                 <  4:457467/149‑1 (MQ=255)
                                                                        agGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGAt                                                                              <  1:123687/149‑1 (MQ=255)
                                                                        agGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGAt                                                                              <  1:622678/149‑1 (MQ=255)
                                                                        agGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGAt                                                                              <  1:5300/149‑1 (MQ=255)
                                                                        agGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGAt                                                                              <  1:115483/149‑1 (MQ=255)
                                                                         ggTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATg                                                                             >  1:126234/1‑149 (MQ=255)
                                                                            cgcgCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTaa                                                                          <  4:72672/149‑1 (MQ=255)
                                                                             gcgcTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAg                                                                         <  3:437175/149‑1 (MQ=255)
                                                                              cgcTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGc                                                                        >  3:455159/1‑149 (MQ=255)
                                                                              cgcTGGGTTGGAATCAACTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTTAATGTCGCGATCGTGACGGTACTCGATGTAAGc                                                                        <  3:392937/149‑1 (MQ=255)
                                                                                 tGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACg                                                                     <  3:598593/149‑1 (MQ=255)
                                                                                   ggTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAg                                                                   <  3:562425/149‑1 (MQ=255)
                                                                                     ttGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTaa                                                                          >  3:78754/1‑140 (MQ=255)
                                                                                     ttGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTaa                                                                          <  4:78754/140‑1 (MQ=255)
                                                                                       ggAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAg                                                               <  1:183022/149‑1 (MQ=255)
                                                                                        gAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTc                                                                                                               >  4:434160/1‑100 (MQ=255)
                                                                                        gAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTc                                                                                                               <  3:434160/100‑1 (MQ=255)
                                                                                           tCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGt                                                                            <  4:275223/132‑1 (MQ=255)
                                                                                           tCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGt                                                                            >  3:275223/1‑132 (MQ=255)
                                                                                           tCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAgttgt                                                           <  2:191696/149‑1 (MQ=255)
                                                                                           tCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAgttgt                                                           >  1:227149/1‑149 (MQ=255)
                                                                                           tCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAgttgt                                                           >  2:583945/1‑149 (MQ=255)
                                                                                            cACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAgttgtt                                                          <  3:276365/149‑1 (MQ=255)
                                                                                             aCCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGTGTTTTTCACGTTAAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgttt                                                         <  4:455159/149‑1 (MQ=255)
                                                                                             aCCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgttt                                                         <  3:548017/149‑1 (MQ=255)
                                                                                               cTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCgacc                                                                                                                                                     >  1:308723/1‑55 (MQ=255)
                                                                                               cTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCgacc                                                                                                                                                     <  2:308723/55‑1 (MQ=255)
                                                                                               cTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTaa                                                                          >  3:140239/1‑130 (MQ=255)
                                                                                               cTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTaa                                                                          <  4:140239/130‑1 (MQ=255)
                                                                                                ttGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGTGTTTTTCACGTTCAATGt                                                                                                    <  4:47106/103‑1 (MQ=255)
                                                                                                ttGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGt                                                                                                    >  3:47106/1‑103 (MQ=255)
                                                                                                  gCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATg                                                    <  4:16992/149‑1 (MQ=255)
                                                                                                      tCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgtttgt                                                       <  3:500829/142‑1 (MQ=255)
                                                                                                      tCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgtttgt                                                       >  4:500829/1‑142 (MQ=255)
                                                                                                      tCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgttt                                                         >  1:521072/1‑140 (MQ=255)
                                                                                                      tCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgttt                                                         <  2:521072/140‑1 (MQ=255)
                                                                                                      tCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgttt                                                         <  3:350877/140‑1 (MQ=255)
                                                                                                      tCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGttgttt                                                         >  4:350877/1‑140 (MQ=255)
                                                                                                            ttgGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACAt                                          <  2:277256/149‑1 (MQ=255)
                                                                                                             tgGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATc                                         <  2:18121/149‑1 (MQ=255)
                                                                                                                 gTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCt                                     >  3:569707/1‑149 (MQ=255)
                                                                                                                    ttCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGAt                                  <  4:569707/149‑1 (MQ=255)
                                                                                                                        ttGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTa                                                      >  4:212770/1‑125 (MQ=255)
                                                                                                                        ttGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTa                                                      <  3:212770/125‑1 (MQ=255)
                                                                                                                        ttGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTc                              <  4:541735/149‑1 (MQ=255)
                                                                                                                          gAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGtttt                                                                                                                 <  4:289030/64‑1 (MQ=255)
                                                                                                                          gAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGtttt                                                                                                                 >  3:289030/1‑64 (MQ=255)
                                                                                                                               aCGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCg                       >  3:37081/1‑149 (MQ=255)
                                                                                                                                     ccgtggtctggtcGAGCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTgg                 >  1:87386/4‑149 (MQ=255)
                                                                                                                                      cgtggtctggtcGAGCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGAc                                                 <  3:423384/114‑1 (MQ=255)
                                                                                                                                      cgtggtctggtcGAGCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGAc                                                 >  4:423384/3‑116 (MQ=255)
                                                                                                                                        tggtctggtcGAGCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGAtc              >  3:266182/1‑149 (MQ=255)
                                                                                                                                        tGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGAtc              <  2:27560/149‑1 (MQ=255)
                                                                                                                                             tggtCGAGCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCAtt         <  2:87386/146‑1 (MQ=255)
                                                                                                                                             tGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCAtt         <  2:518636/149‑1 (MQ=255)
                                                                                                                                               cTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCATTGa       <  4:21930/149‑1 (MQ=255)
                                                                                                                                                tCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCATTGat      >  1:570460/1‑149 (MQ=255)
                                                                                                                                                 cGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCATTGata     <  2:451787/149‑1 (MQ=255)
                                                                                                                                                    ccagaccaGGGATCTCCTGGTTCAGGCGACCGTGTTTTTCACGTTCAATGTCGCGATCGTGTCGGTACTCGATGTAAGCACGAGCGAGTTCTTTGTATGGACCTTACATCAGTTGATTTTCAACTGCGGTCTGGATCTCATTGATATcc  <  3:386509/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATGTGCCTGCACCACGTCACGCGGCAGCAGGTGCTGACGTGCATAGTGTTTAGCCACGATCCCGGCCAGCAGGTCGCGCTGGGTTGGAATCACCTTGCTGTCTTTGTTGGCGTTTTCGTTGAGTAACGAGGCGTTGGTCTGCTCGACCAGACCACGGATCTCCTGGTTCAGGCGACCGCGTTTTTCACGTTCAATGTCGCGATCGTGACGGTACTCGATGTAAGCACGAGCCAGTTGTTTGTATGGACCTGACATCAGCTGATTTTCAACTGCGGTCTGGATCTCATTGATATCC  >  NC_000913/4462197‑4462493

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