Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,368,956 C→T 20.0% L98L (CTA→TTA)  fxsA → suppressor of F exclusion of phage T7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,368,9560CT20.0% 144.3 / 5.6 70L98L (CTA→TTA) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base C (18/36);  new base T (12/3);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 2.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTG  >  NC_000913/4368721‑4369194
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
tttttttGAGGTTTCCTTTGTTTTGGTGGTATTGCTGCCCATCGGGCTTGTTAAATTCACGGCAGTTATCGGTATGTCACCGGTACGTAACCAGGGCTTTAAGAATTTCGTGCAGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑c                                                                                                                                                                                                                                                   <  2:101647/238‑1 (MQ=255)
          gTTGCCCATGTATTGGGGGAATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAACAATGGAGGCGGGTGAAAACCCAGCGTCCGGGGAGAGTTAAAAGGGTTTC‑CTGTGCATTGTTGGTTTGTGTGTTTTTTTACCGTGCGTTTTTACTCTCACCCTC‑GGCCTTTTACCTTTTTTTgc                                                                                                                                                                                                                                          >  2:526732/1‑239 (MQ=255)
               ccATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGGGATGATTAAAAGGGTTTCGCTGATCATTGCTTGGTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTTCT‑CTTTTATTtccgcccg                                                                                                                                                                                                                                    >  1:38883/1‑237 (MQ=255)
                    ttttGGGGGTTTTGGTTGCCCTCTTGTTGGTTTTTTTTACGGTATTTATCGGTATGGCACTGGTACGGAACCAGGGGTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCaa                                                                                                                                                                                                                               <  2:274777/237‑1 (MQ=255)
                    ttttGGGGGTTTTGCTGCCCCGTGTGTTGGTTTTTTTCCAGTCATTTATCGGGATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGTTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCaa                                                                                                                                                                                                                               <  2:678271/237‑1 (MQ=255)
                     aTTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCaa                                                                                                                                                                                                                               <  1:647178/238‑1 (MQ=255)
                              aTTGCTGACCCTGGTGCTGGTTATAGTCACGTCAGTTACCGGTATGCCACTGATACGTAACCAGGGATTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGTTGAAATCCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGa                                                                                                                                                                                                                    <  1:526732/239‑1 (MQ=255)
                              aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGGGATGATTAAAAAGATTTTGCTGATCTTTTCTGCGTTGCTGCCTTTTTTTACGGGCGTTTTTTCCCCCTTCCCC‑GGGCTTTTT‑CTTTTTTTTccgccgtggtcaaaacacccg                                                                                                                                                                                                                      >  2:647178/1‑223 (MQ=255)
                              aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGGTCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGGCTTTTA‑CTTTTTTTTCCGCCCGTGCAAAA‑ACATCTg                                                                                                                                                                                                                     >  1:533747/1‑238 (MQ=255)
                                    taCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGa                                                                                                                                                                                                              <  1:4972/238‑1 (MQ=255)
                                     aCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGGTGATTTAAAAAGTTTCGCCGATCATTGCTGGGTTGCTGCTTTTTTTTCCGGGCTTTTTTATCGACTTCCTC‑GGGCTTTTT‑CTTTTTTTTCCCCCCGGGCAAAA‑GAACCTGaaatggaa                                                                                                                                                                                                             >  2:90846/1‑232 (MQ=255)
                                          cGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCG‑GCGGAGGTGATTTAAAATGGTTTGCTGGTCATTGCTGGTTTGGTGCTTTTATTTACGGGCTTTTTTTACGCCTTTCTC‑GGGCTTTTC‑CTTTTTTTTGCCCCGGTCGTAAA‑ACACTTGCCAGTGA‑AGGTgg                                                                                                                                                                                                       >  2:571291/1‑238 (MQ=255)
                                             ggtGGTTTTTTTCACGTAAGTTTTGGGTATGGTATGTGGACGTAACCAGGGCTTTAAGAATTGCGTGCGGATGCAGCAAAAAATGGCGGCGGGTGAAAGCCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTGTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGc                                                                                                                                                                                                    <  2:20919/237‑1 (MQ=255)
                                                 gTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAAGGTTTTGCTGATCATTGCTGGGTTGCTGCTTTTTTTACCGGGGTTTTTTACCGACTTCCCC‑CGTCTTTTA‑CTTTTTTTGCCCCCGGTGCCAAA‑GCATCTGACAGTGG‑AGTTTATGCCGCa                                                                                                                                                                                                >  1:660842/1‑239 (MQ=255)
                                                       tttccGTTAGTTTTCGGTTTGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGAAAAAAAGGGCGGGGGGTGAAAACCCAGGG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTCCCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTgcg                                                                                                                                                                                          <  2:281888/235‑1 (MQ=255)
                                                           cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTTCCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAACACAT‑TGTCAGTGA‑AGTTTATGCCGCATTTTGGCttt                                                                                                                                                                                      >  1:911200/1‑239 (MQ=255)
                                                           cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGAGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCCGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTTCCCCCGGTGCAAAA‑ACATCTGCCCGTGA‑AGTTGATTCCGCCTTTtcggctt                                                                                                                                                                                      >  1:281888/1‑235 (MQ=255)
                                                           cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGAGTTTCGCTGATCATTGCTGGTTTGCTGCCTTTATTTACGGGCCTTTTTACCGACTTCCTC‑GGTCTTTTA‑CTTTTATTGCCCCCGGTGCAAAA‑ACATCTTACAGTGA‑AGTTGATTCCGCATTTTCGCttt                                                                                                                                                                                      >  1:825937/1‑239 (MQ=255)
                                                               aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTgccgcc                                                                                                                                                                                                                                      >  1:642852/1‑189 (MQ=255)
                                                               aGTTATCGGTATGGCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTgccgcc                                                                                                                                                                                                                                      <  2:642852/189‑1 (MQ=255)
                                                                 ttATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAATGGTTCGCTGATCATTGATGGTTTTGTGCTTTTATTACCCGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTTCCGCCCGTGCCAAA‑AAATCTTCCAATGA‑AGTTTGTGTCGCAATTTGGCtttttttc                                                                                                                                                                                 >  2:323600/1‑234 (MQ=255)
                                                                      ggTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGGTTAAAAAAGTTTCGCTGATCATTGCTGGTTTGCGGCTTTTATTACCGGGCTTTTTTTCCTCCTTCCTC‑GGGCTTCTA‑CTTTTTTTTCCGCCGGTGCAAAA‑ACACCTGACACTGA‑AGGTCATGCCCCATTTCTGCTTTTTTCGCAAGcc                                                                                                                                                                           >  2:160767/1‑239 (MQ=255)
                                                                           gTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCCTTTTTACCGACTTCCTC‑GGTCTTTTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGGGA‑AGGTGATGCCGCATTTTGGCTTTtcttcccttcctgggg                                                                                                                                                                      >  1:464176/1‑226 (MQ=255)
                                                                              aCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGTAAAACCCCGCG‑GCGGAGGTGGTAAAAAGGTGTTCGCTGATCATTGCTGGTTTGTTGCTTTTATTACTGGGCTTTTTTACCGACTTCCCC‑GGGCTTTTA‑CTTTTTTTTCCGCCGGGGCAAAA‑ACATTTGCCCGTGA‑AGTTGATACCCCATTTCTGCTTTTCTTCCATGcgctgggggg                                                                                                                                                                   >  2:802067/1‑230 (MQ=255)
                                                                                    aCGTAACCAGGGCTTTAAGAATTTCGGGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTa                                                                                                                                                            <  1:571291/239‑1 (MQ=255)
                                                                                             gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGGCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACCGTGA‑AGTTTATGCCGCACTTTCGCTTTTTTCGCATGCCTGGCGGC‑GGGTTTTGAG‑CCggg                                                                                                                                                   >  1:162169/1‑238 (MQ=255)
                                                                                             gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCG‑GCGGAGATGATTAAAAGAGGTTCGCTGGTCATTGCTGGTTTTCTGCTTTTATTACCGGGCGTTTTTACCGACTTCCTC‑GGTCTTATA‑CTTTTATTTCCGCCGTGGTAAAA‑ACATCTGCAAGGGAGATTTAATCCCCACTTTGGGCTTTT‑TCGCGTGCCCGGCGGC‑GGGGTTATCG‑CCcggg                                                                                                                                                  >  2:522794/1‑238 (MQ=255)
                                                                                              ggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGGTTGCTGCTTTTTTTATCGGGCTTTTTTACCGGCTTCCTC‑GGGCTTTTA‑CTTTTTTTTCCGCCGGTGCCAAA‑GCACCTTGCAGGGA‑ATGTGAGTCCCCATTTTCGCTTTTCTCGCATGCCCGGCGGCGGGTTTAAGCG‑CCGGGa                                                                                                                                                  >  2:31356/1‑239 (MQ=255)
                                                                                               gCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGa                                                                                                                                                  <  1:249991/237‑1 (MQ=255)
                                                                                                  tttataatttggggtggATGCAGCAAAAATGGGTGGGGGCGGAAACACCAGGT‑GCGGAGATGATTAAAAGTGTGTAGCTGATAGTGGCCGGTTTGCTGCTTTTATGACCGGGCTTTTATGACGACTTCCTC‑GGCCTTCTA‑CTTGTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGGGCTTTTCTCGCATGACTGGCGGC‑GGTTTTAGCG‑CCGTGACCGGt                                                                                                                                             <  2:635937/223‑1 (MQ=255)
                                                                                                  ttAACAATTGCGTGTTGAGGCACAAAAAACTGGAGGCGGCTGAAAAGCCAGGG‑GATGTGATGTTGAAAAGAGTGTCGTTGTTCCTTGCTGGTTTGCGTTTTTTATTTCCGGGCTTTTTGACCGCCTTCATC‑GGTCTTCTT‑CTTTTATTGCCGCCGGTGCAAAA‑ACATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGt                                                                                                                                             <  2:531153/239‑1 (MQ=255)
                                                                                                   tAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCG‑GCGGAGATGATTAAAAGAGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTTTTACCGGGCTTTTTTTCCGACTTCCTC‑GGGCTTTTA‑CTTTTATTTCCGCCCGTGCAAAA‑ACACCTGCCAGGGA‑AGTTGATGCCCCATTTTCGCCTTTTTCCCATGCCTGGGGGG‑GGGTTTATAG‑CCCGGGCCgggg                                                                                                                                            >  2:183625/1‑237 (MQ=255)
                                                                                                        aTTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTg                                                                                                                                                                                                               <  1:127530/170‑1 (MQ=255)
                                                                                                        aTTCCGTGCTGAGGCAGAAAAAAATGGCGGGGGGGGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTTCCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGt                                                                                                                                       <  1:737931/239‑1 (MQ=255)
                                                                                                         tttCGTGGTGATGAAGCAAAAAATGGCGGCGGGTGAAACCCCAGCG‑GCGGAGATGATTAAAAGTGTTTTGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTa                                                                                                                                      <  1:837620/239‑1 (MQ=255)
                                                                                                           tcgtgcgtgtGCAGAAAAAAATGGGGGCGGGTGAAACCCCGGGG‑GGGGAGATGATTAAAAGTGTTTCCCTGACCATTGCTGGTTTGCTGTTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTT‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAAt                                                                                                                                    <  2:464176/230‑1 (MQ=255)
                                                                                                            cGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGGTTTTGCGCCCGGGACCGTGGGGGGTata                                                                                                                                    >  1:173737/1‑237 (MQ=255)
                                                                                                                  gaggcagcaAAAAATGGCGGCGGGTGAAACCCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTg                                                                                                                             <  1:160767/236‑1 (MQ=255)
                                                                                                                       tcaaaaaatTGGCGGCGGGTTAAAACCCGGGG‑GCGGAGTTTATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        <  2:660842/236‑1 (MQ=255)
                                                                                                                        gaaaaaaaTGGGGGGGGGTGAAACCCCGGCG‑GGGGATATGATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        <  2:247574/236‑1 (MQ=255)
                                                                                                                            aaaaTGGGGGCGGGTGAAAACCCGGGG‑GTGGAGATGATTAAAAGTTTTTCCGTGATTATTGCGTGTTTTCTGTTTTTATTACCGGGCTTTTTGACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTa                                                                                                                   <  2:162169/239‑1 (MQ=255)
                                                                                                                             aaaTGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTAc                                                                                                                  <  1:653680/239‑1 (MQ=255)
                                                                                                                                  gggggggggtaaacCCCCGCG‑GTGGAGTTGATTAAAGGTGTTCTGCTGAGCATTGCTGGTTTTCTGCTTTTATTACCGGGTTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                                                                                             <  2:38883/229‑1 (MQ=255)
                                                                                                                                  gcggcgGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                                                                                             <  1:183625/239‑1 (MQ=255)
                                                                                                                                   cgggggGTGAAACCCCATGG‑GGGGGTATGTTTAAAAGTTTTTCCCTTATCATTGCTGTGTTTCTGCTTGTATTACCGGGTTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTTCAAAA‑GCATCTTACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                                                                                             <  2:204441/238‑1 (MQ=255)
                                                                                                                                                   aGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTTACG‑CCGGGAGCCGGGGCGGGGATACTTTTTATGGTGAGGACCAGCGAAAGGGAGATGAggggg                                                                                            >  2:517579/1‑235 (MQ=255)
                                                                                                                                                    ggggggggttttgtttaaagttttttCGCTGATCTTTGCTGGTTTTCTGTTTTTATTCCCGGGCTTTTTGCCGGACTTCCTC‑GGTCTTTTA‑TTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGGTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACTAGCGAAAGGATGATGAGcgcg                                                                                            <  2:173737/216‑1 (MQ=255)
                                                                                                                                                    gcg‑gcgGGTATTATTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                           <  1:473478/239‑1 (MQ=255)
                                                                                                                                                    gcg‑gcgGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                                                                                           <  1:132704/239‑1 (MQ=255)
                                                                                                                                                             gattatTAAAAGTTTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCtt                                                                                    <  1:323600/238‑1 (MQ=255)
                                                                                                                                                                           tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        >  1:352647/1‑188 (MQ=255)
                                                                                                                                                                           tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGt                                                                                                                        <  2:352647/188‑1 (MQ=255)
                                                                                                                                                                           tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCCT‑AAAGACGa                                                                     >  1:632078/1‑238 (MQ=255)
                                                                                                                                                                              cGCTGATCGTTGCTTGTTTGCTGCTTTTATTACCGGGCTTTTTTCCCGATTTCCTC‑GGTCTGCTG‑CTTTTATGGCCTCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CGGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGc                                                                 <  1:66480/239‑1 (MQ=255)
                                                                                                                                                                               gCTGATCTTGGCTGGTTTGCGGCTTTTTTTCCCGGGTTTTTTTACCGACTTCCGC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTCAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGcc                                                                <  1:31356/239‑1 (MQ=255)
                                                                                                                                                                                   aTCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGGGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAgag                                                            >  1:867231/1‑237 (MQ=255)
                                                                                                                                                                                     cATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGTGTACCAGCGAAAGGATGATGACGGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGGtt                                                          >  1:840498/1‑239 (MQ=255)
                                                                                                                                                                                        tGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGTTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGCCCGGTGGCGGGAATACTTTTTTTGGGGAGTACCATCGAAAGGATGAGTAGCGCGCCCGCCCTTATCAC‑AAAAAAGAGCGCCCggggtaag                                                       >  2:751622/1‑233 (MQ=255)
                                                                                                                                                                                             ttttgttttttttttTGCGGGGTTTTTTTCCCGTCTTCCTC‑GGTCTTCTA‑CTTTTATTCCCGCCGGTGCAAAA‑GGATCTGACAGGGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCg                                                  <  2:533747/229‑1 (MQ=255)
                                                                                                                                                                                             gTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGAAGATGAGCGGGACCGCCCTGATCAT‑AAAGACGATCGCCAGGGTTAATG‑TCGa                                                 >  1:43801/1‑239 (MQ=255)
                                                                                                                                                                                                tgctgcTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCCGGATTAATG‑TCGaaa                                               >  1:715531/1‑238 (MQ=255)
                                                                                                                                                                                                 gttttttttttttCCGGGTTTTTTTCCCGACTTCCTC‑GGTCTTCTA‑CTTTTTGTGCCCCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACg                                             <  2:715531/234‑1 (MQ=255)
                                                                                                                                                                                                 gctgctTTTATTCCCGGTCTTTTTTCCCGACTTCCGC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACg                                             <  1:522794/239‑1 (MQ=255)
                                                                                                                                                                                                  ctgctTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAACCGATCGCCAGGATTAATG‑TCGcaaag                                             >  1:599014/1‑236 (MQ=255)
                                                                                                                                                                                                   tttttttATTTCCGGGTTTTTTGACCGCCGTCCCC‑GTGCTTTTC‑CTTTGTTGGCCGCCGGTGCAAAA‑GCATCTGTCGGGGA‑AGGTGATGCCGCATTTGAGCTTTTCTCGCATGCCGGGGGGC‑GGTTTTAGCG‑CTGGGACCGGTGGCGGTAAGCCTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGcc                                           <  2:261016/236‑1 (MQ=255)
                                                                                                                                                                                                    tttttttttctccggtttttttACCCCCTTCCCCGGGTTTTTTT‑TTTTTGTGCCCGCCGGTGCAAAA‑GCATCTGCCAGTGA‑ATTTGATGCCGCATTTGCGCTTTTCTCGCAGGCCTGGGGGC‑GTTTTTGGCG‑CCGGGACCGGTGGCGGTAATCCTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGcc                                           <  2:867231/223‑1 (MQ=255)
                                                                                                                                                                                                               cGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGGGGGGTCCAGCGGAAGGAGGATTGGCGCGAACCCCTTGTTCATAAAAAACAATC‑CCAGGATTAATGTTCAAAACGCCggttatagggg                                >  2:304309/1‑230 (MQ=255)
                                                                                                                                                                                                                  gCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATGGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGCCGGATTATGTGGTTa                             <  1:517579/238‑1 (MQ=255)
                                                                                                                                                                                                                      ttttACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGGGGGTACCAGGGAAAGGGTGAGGAGGGCGACCGCCTTGAGCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACACCGGGTTATGTGGGTATGCCa                        >  2:203637/1‑239 (MQ=255)
                                                                                                                                                                                                                           ccGATTTCCTC‑GGTCTTCTT‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGCCGGATTATGTGGTTATGCCATTTTc                   <  1:90846/239‑1 (MQ=255)
                                                                                                                                                                                                                                             ctA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCCCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGCCGGCTTATGTGGTTATGCCATTTTTCCGGGGTTTTCGTTTTg  >  2:842637/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
CTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCG‑GCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTC‑GGTCTTCTA‑CTTTTATTGCCGCCGGTGCAAAA‑GCATCTGACAGTGA‑AGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGC‑GGTTTTAGCG‑CCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCAT‑AAAGACGATCGCCAGGATTAATG‑TCGAAACGCCGGATTATGTGGTTATGCCATTTTCCGGCGTTTTTCGTTTTG  >  NC_000913/4368721‑4369194

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