Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,992,052 T→A 5.3% D300E (GAT→GAA cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,0520TA 5.3% 423.8 / 7.3 129D300E (GAT→GAAcyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (60/62);  new base A (3/4);  total (63/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

TCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTCGCGCAG  >  NC_000913/3991904‑3992200
                                                                                                                                                    |                                                                                                                                                    
tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAt                                                                                                                                                      >  2:165884/1‑149 (MQ=255)
    cAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATc                                                                                                                                                     <  3:163262/146‑1 (MQ=255)
     aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACTCTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCAt                                                                                                                                                  >  6:343752/1‑148 (MQ=255)
         gTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTg                                                                                                                                              >  8:73910/1‑148 (MQ=255)
             cGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCTAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCatga                                                                                                                                         >  2:349113/1‑149 (MQ=255)
              gATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCatga                                                                                                                                         <  1:352155/148‑1 (MQ=255)
               aTTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCatgatg                                                                                                                                       >  6:249224/1‑149 (MQ=255)
                  ccccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGACCCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTg                                                                                                                                    >  5:319689/1‑149 (MQ=255)
                    ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTg                                                                                                                                    >  4:142272/1‑147 (MQ=255)
                    ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTg                                                                                                                                    >  7:154370/1‑147 (MQ=255)
                    ccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGa                                                                                                                                  <  7:41601/149‑1 (MQ=255)
                      aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGa                                                                                                                                  >  8:160393/1‑147 (MQ=255)
                      aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGa                                                                                                                                  <  7:160393/147‑1 (MQ=255)
                      aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGc                                                                                                                                >  8:209513/1‑149 (MQ=255)
                       tACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCg                                                                                                                               >  1:181332/1‑149 (MQ=255)
                        aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAt                                                                                                                                                      <  5:134259/125‑1 (MQ=255)
                         cAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtg                                                                                                                             <  7:279820/149‑1 (MQ=255)
                          aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtgt                                                                                                                            <  3:183664/149‑1 (MQ=255)
                          aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtgt                                                                                                                            >  4:340496/1‑149 (MQ=255)
                          aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtgt                                                                                                                            >  1:28131/1‑149 (MQ=255)
                          aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCgtgt                                                                                                                            >  5:37874/1‑149 (MQ=255)
                          aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGGTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtgt                                                                                                                            <  4:340744/149‑1 (MQ=255)
                             gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTAc                                                                                                                         <  1:165884/149‑1 (MQ=255)
                              cGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACt                                                                                                                        <  2:253420/149‑1 (MQ=255)
                               ggTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACt                                                                                                                        <  6:37874/148‑1 (MQ=255)
                                     gAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATAcc                                                                                                                 >  7:302714/1‑149 (MQ=255)
                                      aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTa                                                                                                                          <  2:286320/139‑1 (MQ=255)
                                      aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTa                                                                                                                          >  1:286320/1‑139 (MQ=255)
                                      aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCt                                                                                                                >  3:118405/1‑149 (MQ=255)
                                      aaaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGACGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATa                                                                                                                                                 >  8:65969/1‑116 (MQ=255)
                                        aaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGa                                                                                                              >  1:283087/1‑149 (MQ=255)
                                        aaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGa                                                                                                              <  5:218312/149‑1 (MQ=255)
                                         acacTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGAc                                                                                                             <  2:168645/149‑1 (MQ=255)
                                               ctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTTTTACTGAATACCTGACGGCGAt                                                                                                       <  4:133089/149‑1 (MQ=255)
                                                  ctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGa                                                                                                    <  4:266981/149‑1 (MQ=255)
                                                     gAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTg                                                                                                     >  7:206167/1‑145 (MQ=255)
                                                     gAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGGCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGa                                                                                                 <  4:341236/149‑1 (MQ=255)
                                                        gCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAttt                                                                                              <  7:209513/149‑1 (MQ=255)
                                                         ccTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGAtttt                                                                                             <  8:154370/149‑1 (MQ=255)
                                                          cTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                                            >  1:159593/1‑149 (MQ=255)
                                                           tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                           >  8:294042/1‑149 (MQ=255)
                                                           tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAc                                                                                           >  6:71006/1‑149 (MQ=255)
                                                                cTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGTAGAGTGTTACTGAATACATTACGTCGATTGAAGATTTTACCCGTc                                                                                      >  3:240126/1‑149 (MQ=255)
                                                                 tGGGAATACCAGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCt                                                                                     >  6:227668/1‑149 (MQ=255)
                                                                    gAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGa                                                                                  <  7:194995/149‑1 (MQ=255)
                                                                    gAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGa                                                                                  <  8:302714/149‑1 (MQ=255)
                                                                      aTACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGGGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGAtt                                                                                <  2:181332/149‑1 (MQ=255)
                                                                      aTACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGAtt                                                                                <  8:191969/149‑1 (MQ=255)
                                                                       tACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGAttt                                                                               <  5:317020/149‑1 (MQ=255)
                                                                         cccGAACCCTCGTCTGCTGTCTAAAGATTTCAAACAGCGTTTGCACGAGGGCGTGATTGTATCGTTTTTTCTCTATCCATTCTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTa                                                                              <  4:118405/148‑1 (MQ=255)
                                                                         cccGAACCCACTTCTGCTGTCTAAAGATATCAAACAGCGTTTGCACGACGGCTAGATTGTATCGTTTGGTCTCTATCCATACTTCATGATTCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTa                                                                              <  4:113299/148‑1 (MQ=255)
                                                                         cccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAg                                                                             <  5:71006/149‑1 (MQ=255)
                                                                          ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAg                                                                             >  4:143071/1‑148 (MQ=255)
                                                                               cccACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGAACCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTc                                                                       <  5:249224/149‑1 (MQ=255)
                                                                                  aCGTCTGCTGGCTAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTg                                                                                                                                    <  6:348811/85‑1 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtg                                                                                                                             >  5:188920/1‑92 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtg                                                                                                                             <  6:188920/92‑1 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                    >  2:116253/1‑149 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                    <  3:279948/149‑1 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                    >  8:245762/1‑149 (MQ=255)
                                                                                  aCGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCAATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCgct                                                                    <  2:283087/149‑1 (MQ=255)
                                                                                      ctgctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCtt                                                                <  5:74024/149‑1 (MQ=255)
                                                                                       tgctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtg                                                                                                                             >  8:388161/1‑87 (MQ=255)
                                                                                       tgctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCgtg                                                                                                                             <  7:388161/87‑1 (MQ=255)
                                                                                         ctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                                            >  6:309293/1‑118 (MQ=255)
                                                                                            gCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTAtt                                                           <  7:120842/148‑1 (MQ=255)
                                                                                            gCGAAAGATATCAAACAGCGTTGGCACGACGGCGAGATTGTATCGGTGGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTAttt                                                          <  1:300324/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                      <  6:227398/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                      <  5:206268/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                      <  8:115896/149‑1 (MQ=255)
                                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                      <  8:55737/149‑1 (MQ=255)
                                                                                                 aGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAg                                                     <  7:294042/149‑1 (MQ=255)
                                                                                                  gATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAgt                                                    <  3:142272/149‑1 (MQ=255)
                                                                                                    tatCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAgtgt                                                  >  8:40659/1‑149 (MQ=255)
                                                                                                    tatCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAgtgt                                                  >  5:143758/1‑149 (MQ=255)
                                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGCTTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGtgtg                                                 <  5:227668/148‑1 (MQ=255)
                                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGtgtg                                                 <  1:116253/148‑1 (MQ=255)
                                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGtgtg                                                 <  5:288667/148‑1 (MQ=255)
                                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGc                                                >  3:27503/1‑149 (MQ=255)
                                                                                                           cAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaaaa                                           <  2:367446/149‑1 (MQ=255)
                                                                                                                     acgacgGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGt                                 <  7:40659/149‑1 (MQ=255)
                                                                                                                      cgacgGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTg                                <  1:349113/149‑1 (MQ=255)
                                                                                                                            cGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcgcg                          >  1:45336/1‑149 (MQ=255)
                                                                                                                             gagaTTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGGGAAAAGCTCAGCCGTGAACgcgc                         >  1:371319/1‑149 (MQ=255)
                                                                                                                                 ttGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGc                     >  8:17496/1‑149 (MQ=255)
                                                                                                                                   gTATGGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGCCGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg                    <  4:240126/148‑1 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAATTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg                    >  7:109742/1‑147 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCg                    >  3:17109/1‑147 (MQ=255)
                                                                                                                                    tATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGt                   <  2:45336/148‑1 (MQ=255)
                                                                                                                                      tCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTaaaa                                                      <  4:86998/111‑1 (MQ=255)
                                                                                                                                      tCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGaaaa                                           >  7:90061/1‑122 (MQ=255)
                                                                                                                                        gTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGc               <  4:27503/148‑1 (MQ=255)
                                                                                                                                             gTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGt         >  8:159830/1‑149 (MQ=255)
                                                                                                                                                tcGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTcgc      >  8:63444/1‑149 (MQ=255)
                                                                                                                                                tcGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGCCTTGCGTagc      >  6:305/1‑149 (MQ=255)
                                                                                                                                                  gATCCATACTGGATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGGCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGGGTGCGAAAAGCTCAGCCGTGAACCCGCCTGCGTAGGCTGGCGTcg       <  3:364441/146‑1 (MQ=255)
                                                                                                                                                  gATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTcgc      <  7:104877/147‑1 (MQ=255)
                                                                                                                                                  gATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTcgc      >  8:104877/1‑147 (MQ=255)
                                                                                                                                                  gATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACATGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGAGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGTCTTcgct     >  5:315719/1‑147 (MQ=255)
                                                                                                                                                   aTCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTCGCGCa   <  3:123986/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTAGGCTGGCGTCGCGCAG  >  NC_000913/3991904‑3992200

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: