Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,952,942 A→G 57.1% T154A (ACC→GCC)  ilvE → branched‑chain amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,952,9420AG57.1% 60.1 / 98.6 105T154A (ACC→GCC) ilvEbranched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base G (30/29);  total (54/51)
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 106 total aligned reads displayed.

GATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGTTTGA  >  NC_000913/3952795‑3953087
                                                                                                                                                   |                                                                                                                                                 
gATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAacac                                                                                                                                                  <  1:504132/149‑1 (MQ=255)
  tGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCCCCAAACACCa                                                                                                                                                >  2:405208/1‑149 (MQ=255)
  tGTTGGCATGGGAGTAAACCCGCCAGCGGAATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGTAGCACCAAACACCa                                                                                                                                                >  1:472873/1‑149 (MQ=255)
   gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCGTCCTGGAACCGCGCAGCACCAAACACCAt                                                                                                                                               <  1:459457/149‑1 (MQ=255)
   gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                                                                                               <  1:146354/149‑1 (MQ=255)
   gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                                                                                               <  1:217111/149‑1 (MQ=255)
    ttGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATc                                                                                                                                              >  2:680020/1‑149 (MQ=255)
      ggCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGTAGGGGATCGATGGGAGGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATccc                                                                                                                                            >  1:325960/1‑149 (MQ=255)
      ggCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATccc                                                                                                                                            >  1:197355/1‑149 (MQ=255)
       gCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCg                                                                                                                                           >  1:359177/1‑149 (MQ=255)
                  aaCCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCaaa                                                                                                                                <  2:6410/149‑1 (MQ=255)
                    cccGCCAGCGGGATACTCAACCGACGTGATGATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGCTGCGATGGTTTCCTCCGGGAACCGCGCAGCACCAACCCCCAGCCCCACGGCGGCAAAAg                                                                                                                              >  2:519274/1‑149 (MQ=255)
                     ccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCGGTGGGGAGCGTATATGGGCGCAGAAGCGCTGGAGCAGGGGATCGATACGATGGTTTCCTCGTGGAACCGCGCAGCACCAAAAGCCATCCCCACGACGGCAAGAGc                                                                                                                             >  2:262336/1‑149 (MQ=255)
                     ccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGc                                                                                                                             <  2:297202/149‑1 (MQ=255)
                     ccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGc                                                                                                                             <  2:297241/149‑1 (MQ=255)
                       gCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCCCCAAACGCCATCCCGACGGCGGCAAAAGCCg                                                                                                                           >  2:402037/1‑149 (MQ=255)
                       gCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAACAGCCg                                                                                                                           >  2:504626/1‑149 (MQ=255)
                          atccggATAATCAACCGTCGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCggtg                                                                                                                        <  1:137408/145‑1 (MQ=255)
                          agtgggATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAACAGGGGATCTATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                                                                        <  2:149092/146‑1 (MQ=255)
                          aGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCATGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                                                                        <  2:491218/149‑1 (MQ=255)
                            cGGGATACTCAACAGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCCGCCCCAACCGCCATCCCGACGGCGGCCAAAGCCgggggt                                                                                                                      >  2:457323/1‑149 (MQ=255)
                              ggATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTac                                                                                                                    >  2:493671/1‑148 (MQ=255)
                               gATACTCCACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGTTGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAAc                                                                                                                   <  2:36703/149‑1 (MQ=255)
                               gATACTCAATCTACGTGATTAGCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAATGCCGGTGGTAAc                                                                                                                   <  2:639933/149‑1 (MQ=255)
                                   cTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                                                                                                               <  2:130463/149‑1 (MQ=255)
                                    tCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCCGCCCCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAccc                                                                                                              >  2:16713/1‑148 (MQ=255)
                                    tCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTacctacct                                                                                                              >  2:171773/1‑149 (MQ=255)
                                    tCAACCGACGTGATTATCGCTGCTTTCCCGGGGGGAGCGTATCTGGGCGCAGAAGCGCCGGAGCCGGGGATCGATGCGAGGGCTTCCTCCTGGAACCGCGCAGCCCCAAACACCGTCCGCCCGGCCGCACAAGCCGGCGGTAAccccct                                                                                                              >  2:16865/1‑144 (MQ=255)
                                      aaCCGACGTGATTATCGCTGCTTTCCCGTGGGGCGCGTATCTGGGCGCAGAAGCGATGGAGCCGGGGAACGATGCGATGGTTTCCTCCTGGACACGCGCAGCACCAAACGCCATACCGAAGGCGGCAGAAGCCGGGGGTAACTCCctct                                                                                                            >  2:104291/1‑149 (MQ=255)
                                        ccGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctct                                                                                                          <  2:359177/149‑1 (MQ=255)
                                           aCGGGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCATCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTcc                                                                                                       <  2:649141/149‑1 (MQ=255)
                                             gTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCCCTCTTCCct                                                                                                     >  2:548696/1‑149 (MQ=255)
                                              tGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                                                                                    <  1:312840/149‑1 (MQ=255)
                                               gATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGAACGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgc                                                                                                   >  2:554552/1‑149 (MQ=255)
                                                 ttctCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgctg                                                                                                 <  2:78137/146‑1 (MQ=255)
                                                   ataaatcctTGCCCGTGGGGAGCGTAACGGGACTCATAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggt                                                                                               <  1:544422/142‑1 (MQ=255)
                                                    tCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                                                                                              <  1:674873/149‑1 (MQ=255)
                                                     cGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtgg                                                                                             <  1:496061/149‑1 (MQ=255)
                                                       ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                                                                                           >  1:346748/1‑149 (MQ=255)
                                                       ctgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGCCGGCAAAAGCCGGTGGCAACTACCCCCCTCCCCTGCCGGTGGGt                                                                                           >  1:9236/1‑149 (MQ=255)
                                                        tgctTTCCCGGGGGGAGGGGATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAGGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTa                                                                                          <  2:176836/149‑1 (MQ=255)
                                                         gctTTCCCGTGGGGAGCGGATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                                                                                         <  1:138088/149‑1 (MQ=255)
                                                           tttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGGAGCg                                                                                       >  1:485690/1‑149 (MQ=255)
                                                             tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                                                                                     >  1:494345/1‑149 (MQ=255)
                                                             tCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGaa                                                                                     <  1:680020/149‑1 (MQ=255)
                                                              ccccTGGGGAGCGTACCTCGGCGTAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCACCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                                                                                    <  2:472454/149‑1 (MQ=255)
                                                              cccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGTTGGTGGGTAGCGAAg                                                                                    <  1:487980/149‑1 (MQ=255)
                                                               ccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgc                                                                                   <  2:331073/149‑1 (MQ=255)
                                                                 gTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgc                                                                                 <  1:327095/149‑1 (MQ=255)
                                                                    gggAGCGTATCGGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgcc                                                                              <  1:27440/149‑1 (MQ=255)
                                                                     ggggCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccg                                                                             <  2:378894/146‑1 (MQ=255)
                                                                      gAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGcgccgc                                                                            >  1:626381/1‑149 (MQ=255)
                                                                               tGGGCGCAGAAGCGCTGGAGCAGGGGATCGAGGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAt                                                                   <  1:688028/149‑1 (MQ=255)
                                                                               tGGGCCCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAt                                                                   <  1:542187/149‑1 (MQ=255)
                                                                                     cAGAAGCGCTGGGGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGaa                                                             <  1:364232/149‑1 (MQ=255)
                                                                                     cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCCCCGAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTAGCAGGaa                                                             >  1:344752/1‑149 (MQ=255)
                                                                                      aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTGCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAg                                                            >  2:468487/1‑149 (MQ=255)
                                                                                      aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAg                                                            >  1:352305/1‑149 (MQ=255)
                                                                                      aGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAg                                                            >  1:163597/1‑149 (MQ=255)
                                                                                         aGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACTAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTa                                                         >  2:403189/1‑149 (MQ=255)
                                                                                           cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATc                                                       >  1:611082/1‑149 (MQ=255)
                                                                                           cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATc                                                       >  1:148831/1‑149 (MQ=255)
                                                                                                  ggaGGGGACCGATGCGATGGTTTCCTCCTGGAACCGCGCTGCACCAAACACCATCCCGACGGCGGCAAAAGCCAGTATTATCTATCTCCGGTCCCTGCCGGTAGATATCTAAGCGCGCCGCCTCGGTTATCAGTAAGGTATCGCGCTgg                                                <  2:344752/147‑1 (MQ=255)
                                                                                                      gggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGAtgt                                            <  1:635088/149‑1 (MQ=255)
                                                                                                      gggATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAACGAAGCGCGCCGCCACGGTTGTCAGGAAGGTATCGCGCTGGAtgt                                            >  1:18247/1‑149 (MQ=255)
                                                                                                        ggtCGGTGCGGTGGTTTCCTCCTGGAACCGGGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGGTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGa                                          <  1:329793/147‑1 (MQ=255)
                                                                                                          tCGAGGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAAc                                        <  1:725723/149‑1 (MQ=255)
                                                                                                            gatgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACgg                                      <  1:442798/149‑1 (MQ=255)
                                                                                                             atgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGAAGCAATGCAGATCAGGAAGGTATCGCGCTGGATGTGAACGGt                                     <  1:554552/149‑1 (MQ=255)
                                                                                                              tgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGtt                                    <  2:738951/149‑1 (MQ=255)
                                                                                                              tgcgatgGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGATGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGGGAACGGtt                                    >  2:546865/1‑149 (MQ=255)
                                                                                                                    ggTTTCCTCCTGGAACCGCGCAGCACCAATCACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATAtct                              >  1:457719/1‑149 (MQ=255)
                                                                                                                     gtatccaccTGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATAtctc                             <  2:428265/146‑1 (MQ=255)
                                                                                                                     gTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATAtctc                             <  1:341436/149‑1 (MQ=255)
                                                                                                                      tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATctct                            <  2:550409/149‑1 (MQ=255)
                                                                                                                      tttCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGGGAACGGGTATATctct                            >  1:43281/1‑149 (MQ=255)
                                                                                                                        tcctccTGGAACCGCGCAGCACCAAACACAACCGCAACGCCCGCAAGAGCCCGTGGTAACGACCTTCCGTCCCTTGTTGTGCGTGTGTAAGCGCTGCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGa                          <  2:438552/149‑1 (MQ=255)
                                                                                                                        ccctcctAGCACCCCGCAGCACCAACCCCCATCCCGACGGCGGCAAAAGCCCGTAGTAACTCCCCTTCTTCCCTCCGTGTTGCTAGCTAACCGCGCCGCCACGGTTATCACGAAGGTATCGCTCTCGATGTAAACAGTTAAAGctctca                          <  2:304457/148‑3 (MQ=255)
                                                                                                                         cctcctGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGaa                         >  1:562904/1‑149 (MQ=255)
                                                                                                                         cctcctGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGaa                         <  2:742000/149‑1 (MQ=255)
                                                                                                                         cctccgGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGaa                         <  1:306874/149‑1 (MQ=255)
                                                                                                                          ctcctGCAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAg                        <  1:454174/149‑1 (MQ=255)
                                                                                                                           tcctGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGCACCTACCTCTCTTCCCTGCTGGTGGGTAGCGACGCGCTGTGCGTCTGGTAGCAGGAAGGTATCGCGCTGGATGTGAACGGGTACATCTCTGAAgg                       >  1:738876/1‑149 (MQ=255)
                                                                                                                            cctGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGCAGgc                      >  1:43528/1‑149 (MQ=255)
                                                                                                                            cctGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGgc                      <  1:632975/149‑1 (MQ=255)
                                                                                                                             ctGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGGGAACGGTTATATCTCTGAAGgcg                     >  1:274484/1‑149 (MQ=255)
                                                                                                                              tGGAACCGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGgcgc                    >  1:482718/1‑149 (MQ=255)
                                                                                                                               ggAACCGCGGACCACCAAACGCCAGCCCGACGGCGGCAACAGCCGGTGGTAACTTCCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCa                   <  1:104291/149‑1 (MQ=255)
                                                                                                                                gAAACGCGCAGGACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGGACGGTTATATCTCTGAAGGCGCAg                  >  2:381864/1‑149 (MQ=255)
                                                                                                                                 aaCCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCCCTCTTCCCTGCTGGTGGGTAGCGAAGAGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAgg                 <  1:623270/149‑1 (MQ=255)
                                                                                                                                   ccGCGCAGCACCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGGGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCg               >  2:380375/1‑149 (MQ=255)
                                                                                                                                     gcgcAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGaa             >  1:130573/1‑149 (MQ=255)
                                                                                                                                       gcagcaCCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGaaaa           <  1:634001/149‑1 (MQ=255)
                                                                                                                                        cagcaCCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAAc          >  2:47008/1‑149 (MQ=255)
                                                                                                                                           caCCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTg       >  2:68245/1‑149 (MQ=255)
                                                                                                                                           caCCAAACGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGAGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATAGCTCTGAAGGCGCAGGCGAAAACCTg       >  1:463785/1‑149 (MQ=255)
                                                                                                                                           caCCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCGGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTAGCGCGCTGGTTGTGAACGGTTATATCTCGGAAGGCACAGGCGACAACCTg       >  1:38372/1‑149 (MQ=255)
                                                                                                                                            aCCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGCTATACCTCTGAAGGCGCAGGCGAAAACCTGt      <  1:405208/149‑1 (MQ=255)
                                                                                                                                             ccAAACACCATCCCGACGGCGGCAAAAGCCGGAGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGtt     >  1:442669/1‑149 (MQ=255)
                                                                                                                                                aaCGCCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGTTTGa  >  2:422482/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGATGTGAACGGTTATATCTCTGAAGGCGCAGGCGAAAACCTGTTTGA  >  NC_000913/3952795‑3953087

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