Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,933 G→A L343L (CTA→TTA)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,9330GA88.6% 116.8 / NA 44L343L (CTA→TTA) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/7);  minor base C (0/5);  total (32/12)
Fisher's exact test for biased strand distribution p-value = 7.29e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AGCACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACC  >  NC_000913/6700‑7173
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
aGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGTGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTc                                                                                                                                                                                                                                    <  1:550022/250‑1 (MQ=255)
aGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAACCCAGGAACCCATT‑AACACCCGCATCGCCTTCt                                                                                                                                                                                                                                   >  1:19329/1‑251 (MQ=255)
aGCCATGCAGGCCATTATGATAGCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCAACGCCTTc                                                                                                                                                                                                                                    <  2:2197/250‑1 (MQ=255)
       cAGGCCATTACGATCTCTGCCGGTTGCCACATCAGCGGAATACTTAGCAAGGTGACGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGCTCGCTTTAGGGATGTTCCGGCGTACAAAGCAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGAAAACAGAATAACCACGAGGGGAACAAACTCAGCACCCCAGGAACCCATT‑AACACCGGCATCGCCTTCTGGA‑Ta                                                                                                                                                                                                                             >  2:304487/1‑249 (MQ=255)
                     atCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCAGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGAAAAAAACTAAGCACCCCAGGA‑CCCATTTAACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTcc                                                                                                                                                                                                              >  1:572259/1‑250 (MQ=255)
                         gCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGAGCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCgg                                                                                                                                                                                                         <  1:399113/251‑1 (MQ=255)
                           cAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCa                                                                                                                                                                                                       <  1:155584/251‑1 (MQ=255)
                            aGTTGCCACATCAGCGGACGCCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCAt                                                                                                                                                                                                      >  1:472022/1‑251 (MQ=255)
                            aGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCAt                                                                                                                                                                                                      >  2:394866/1‑251 (MQ=255)
                                     aTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTtgt                                                                                                                                                                                             >  1:136287/1‑251 (MQ=255)
                                                acTTAGAAAGGCGCCGCCGCTGACCGTAGCGCAGGGGCAGAACCGCAAAGACCCGCTCGCTTCAGGGTTGTTCAGGCATAAAAAGAAGAGATTGCTTTA‑GGCGTAAATGAAGTTGGCAACGATGGAGATGAAGGCAAACAGAATAACCACGAGGGAAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGATTGTGCCGTT‑AACCg                                                                                                                                                                                 >  2:410081/2‑251 (MQ=255)
                                                      cAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGt                                                                                                                                                                            >  2:219594/1‑250 (MQ=255)
                                                               gccgATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGc                                                                                                                                                                   >  1:306763/1‑250 (MQ=255)
                                                                      aCCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGAATGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCt                                                                                                                                                           >  1:169411/1‑251 (MQ=255)
                                                                         gTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGc                                                                                                                                                        >  1:224047/1‑251 (MQ=255)
                                                                            gCGAAGGGGCAGACCCGCAACCAACAGCGCGTTTTGGGGTTTTTCAGGGGTAAAAAGAAGAGATGGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCAAAAGGCAAACAAAATAACCACGAGGGTAACAAACTCACCCCCCCATGACCCCATT‑ACCCCCCGCATCGCCTTCGGAA‑AAAGCTGAATCCCTTCCAGCGGCAGGTAGGTTGTGCCGTT‑CCCCGCCAGTAATATCAGCATGGCGCTTGCCGt                                                                                                                                                     <  1:582244/251‑1 (MQ=255)
                                                                               aaGGTGCAGATCCGAAAACACAGGTTCTTTTTGGGGTTGTCCAGGCGAAAAAAGAAGAGTTTTTTTGC‑GGCTTAAATGTAGTTGGCAACGATGGAGCGGAGGGAAAACAGAATAACCCCGGGGGTACCAAACCCAGCCCCCCAGGACCCCTTT‑ACCCCCCGCAGCGCCTTCTGGA‑TAAGCTGAATCCCTCCCAGCGGCATGTAGGTTGGGCCGTT‑CCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACa                                                                                                                                                  <  1:394866/251‑1 (MQ=255)
                                                                                            gCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAACAGAAGAGCTTGTTTTC‑GGCGTCAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGACTAACCACGAGGGTAACACACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGgtg                                                                                                                                     >  2:394948/1‑251 (MQ=255)
                                                                                               aaCACCAGACCGTTTTAGGGTGGTTCGGGCAAAAAAAAAAGAAATTTTTTCC‑GGCGTAAATGAAGTGGGCAACGATGGAGCGAAAGGCAAAAAAAATACCCACGAGGGTAACAAACCCACCCCCCCAGGACCCCATT‑ACCCCCCGCATCCCTTTCTGGA‑AAAGCAGAATCCCTCCCAGCGGCATGTGGGTTGGGCCGTT‑CCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCg                                                                                                                                  <  1:468409/251‑1 (MQ=255)
                                                                                                  cccAGACCGTTTTAGGGTTGTTCGGGCAAAAAAAAAAGAGAT‑GTTTTCCGGCGAAAATGTAGTGGGAAACGAGGGAGCGGAAGGAAAACAAAATAACCACGGGGATAACAACCTCAGCACCCCAGAACCCCATT‑ACCCCCCGACCCGCCTCCGGGA‑AAACCTAAATCCCTCCCGGCGGCATGTGGGTTGTGCCGTT‑CCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATa                                                                                                                               <  2:230600/250‑1 (MQ=255)
                                                                                                      gATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATaaaa                                                                                                                            >  1:519143/1‑250 (MQ=255)
                                                                                                               aGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatc                                                                                                                   >  1:85980/1‑250 (MQ=255)
                                                                                                                 ggTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatcatc                                                                                                                >  1:114360/1‑251 (MQ=255)
                                                                                                                       tCAGGCGAAAAAAGAAGAGATTGTTTGC‑GGCGAAATTGTAGTTGGCAACGAGGGAGCGGAAGGCAAACAGAATAACCACGAGGGTAACAACCTCAGCACCCCAGGAACCCTTT‑AACCCCCGCATCTCCTTCGTGA‑TAAGCTGAATCCCTCCCAGCGGCATGGGGGTTGTGCCTTC‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACTAGGGTGTCGATAAAAATGCCAATCATCTGGAc                                                                                                           <  1:446226/250‑1 (MQ=255)
                                                                                                                                aaaaGAAGCGATAGTTTTC‑GACGTAAATGTAGTTGCCAACGATGGAGCTGAAGGCACACAGAATAACCACGAGGGTAACAAACTCCGCACCCAAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGCATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACACGCCAGTACTATCAGCCTGGCGCTTGCCGTACAGATGACCAGGGGGCCGATAAAACTGCCAATCATCTGGACAATCCCTTgc                                                                                                 >  2:319605/1‑251 (MQ=255)
                                                                                                                                           ttGTTTTC‑GGCGTAAATGTAGTTGGCAACGCTGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATg                                                                                      >  2:259555/1‑251 (MQ=255)
                                                                                                                                           ttGTTTTC‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCCGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGAGTCGATAACAATGCTAATCATCTGGACAATCCCTTTCGCTGCCGGATg                                                                                      >  2:221938/1‑251 (MQ=255)
                                                                                                                                            ttttttC‑GGCGTAAATGTAGTGGGCAAGGAGGGAGCTGAAGGCAAAAAGAATAACCACGAGGGAAACAAACTCCGCCCCCCAGGAACCCATT‑AACCCCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCGGCGGAATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATg                                                                                      <  1:484517/248‑1 (MQ=255)
                                                                                                                                                  c‑GGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGCCCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGc                                                                               >  1:56730/1‑251 (MQ=255)
                                                                                                                                                     gcgtAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGAAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGcc                                                                              >  2:541696/4‑250 (MQ=255)
                                                                                                                                                         aaaTGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCAGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGa                                                                          <  1:345721/250‑1 (MQ=255)
                                                                                                                                                                   ggaaaagagggAGCTGAAGTCAAACAAAATACCCACGGGGGTAAAAAACCCCGCACCCCGGAAACCCATT‑ACCCCCCCCCCCGCCTT‑TGGAATAAGCAGAATCCCTCCCAGGGGATTGTAGGTTG‑GCCCTTGCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACgccgctgcc                                                                <  2:492718/241‑1 (MQ=255)
                                                                                                                                                                     cAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACgccgctgccgct                                                             >  1:157195/1‑251 (MQ=255)
                                                                                                                                                                      aaCGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGGGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGccgctgccgctg                                                            >  2:168602/1‑251 (MQ=255)
                                                                                                                                                                              aGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGttt                                                    >  1:597039/1‑251 (MQ=255)
                                                                                                                                                                              aGCTGAAGGCAAACAGAACAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACAACCGCATAGACTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGATGTACCGTT‑ACCCGCCCGTAATATCAGCAAGCCGCTTGCCGTACAGATGACCAGGGCGTCGATAAAAATGCCAATACTCTGTACAATCCCTCGCGCCGCCGGAGTAGGAGGACATGACGCCGATGCACCAGCCGCGttt                                                    >  2:73588/1‑251 (MQ=255)
                                                                                                                                                                                 tGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTCGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGAGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAAAGCCAATCCTCTGGAAAAACACTTGCGCTGCCGGATGCGGACGCCAGGACGCCGCTGCCGCCGCCGCGTTTGGc                                                 >  2:484529/1‑251 (MQ=255)
                                                                                                                                                                                  gAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCCGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATTCGGAGGCCAGGACGCCGCTGCCGCTGCCCCGTTTTGCg                                                >  2:343921/1‑251 (MQ=255)
                                                                                                                                                                                   aaGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGGGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGATGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGt                                               >  2:1393/1‑251 (MQ=255)
                                                                                                                                                                                           cAGAATAACCCCGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACgccgctg                                                                  >  2:584739/1‑224 (MQ=255)
                                                                                                                                                                                           cAGAATAACCACGAGGGTAAAAAAATCAGCACCCCAGGCACCCATT‑AACAACCGCAGCGCCTTCTGGC‑TACGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAAAATCCGCATGGAGCTTGGCGTACAGATTACCAGGGGGTCGATAAAAATGCAAATACTCTGGACAATCAATTGCGCTGCCGGAGGAGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCgcccccc                                       >  2:444803/1‑245 (MQ=255)
                                                                                                                                                                                                     acgaGGGAAACAAACTCCGCACCCCAGGAACCCATT‑AACACCCGCACAGCATTATGGA‑TAAGCTGAATAACTTCCAGCGGCATGTAGGTTGTGCCGTT‑AACCGCCAGTAATATCAGAATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCACTTGCGCTGCCGGATGCGGCGGCCAGGACGCAGTTGCCGCTGCCCAGTTTGGCGTccgaacggaagagcgtcg                             >  2:420982/4‑234 (MQ=255)
                                                                                                                                                                                                                aaaCTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAAc                   >  2:172526/1‑250 (MQ=255)
                                                                                                                                                                                                                                gAACCCATT‑AACACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑AACCGAAAGTAATATCAGAATGGCGCTTGCCGTACAGCTGAACAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGACCCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAAcc  >  1:591513/1‑251 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
AGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTC‑GGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AGCACCCGCATCGCCTTCTGGA‑TAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTT‑ACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGCGTCGAACCCATTCCCGCCTCATTGGAAAACATACTGCGCTGAAAACC  >  NC_000913/6700‑7173

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