Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 347,833 G→A A204V (GCC→GTC)  prpR ← DNA‑binding transcriptional dual regulator PrpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913347,8330GA97.6% 139.7 / ‑3.9 41A204V (GCC→GTC) prpRDNA‑binding transcriptional dual regulator PrpR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/26);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGG‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGTCGAGGCGCAAATTAAAGGTTTTTTGAA  >  NC_000913/347589‑348079
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
cGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGAAAAATATTCCCGATGAATCGCCGGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGt                                                                                                                                                                                                                                                  <  1:164382/251‑1 (MQ=255)
 gCCCCGCAGTTGACGGCAAAAAACGGATGCGGCTGTTGGCCCGGTCGCGCAGCGTGGCGGGCAAAATATTCCCGATAAATCGCCTGGGCCGCAAGTGTTTTCCCCGTCCCCTTTCCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGTCGTTGCGTGTGA‑CATCGTg                                                                                                                                                                                                                                                 <  2:250004/251‑1 (MQ=255)
   cccGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCAGGGCCCCCACACCTTTGCCCGACCCCGTTTCCCCCTCACTCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGCAAAGAGTACGCCGGGCGTTGCGGGTGA‑CATCGTggg                                                                                                                                                                                                                                               >  1:385993/1‑249 (MQ=255)
                             gCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTTCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACAACGCCGCACTGGAGCGGGCATACAGCCAAATAGTCTGCCGTACTTGTTCCATCTGTGGGGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTc                                                                                                                                                                                                                     >  1:251852/1‑251 (MQ=255)
                                acttttcacccTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAACCGCCTGTGCCGCCAGTTTTTGGACCGTCCCGGTTTCCCCCTAAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTTCGCAGGGCGTTGCGG‑TGGACATCGTGAGTGTTAGGGCGTAACGACATGCGCGTCata                                                                                                                                                                                                                  <  2:19463/243‑1 (MQ=255)
                                      tGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAg                                                                                                                                                                                                             <  1:310873/250‑1 (MQ=255)
                                            gTCGCGCATCGTGGCGGGAAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCAt                                                                                                                                                                                                        <  1:42814/249‑1 (MQ=255)
                                                cgcATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCCTCTGTGGTGATTGACCGAGCATATCGCCCCGCACGTAACGAGTACGCAGGGCGGTGCGGGTGA‑CATCGTGAGGGTTATGGCGTAACGCCCTGCGCGGCCTATCCAGCGCCTCGCTg                                                                                                                                                                                                   >  2:153837/1‑250 (MQ=255)
                                                      gTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCt                                                                                                                                                                                            <  1:240759/251‑1 (MQ=255)
                                                          cGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGgcg                                                                                                                                                                                        >  1:165496/1‑251 (MQ=255)
                                                                      tcccgagtaagcGCCTGGGCCCCCAGCTCTTTGCCCGTCCCCGTTCCCCCCTCAACCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTCCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTggcgg                                                                                                                                                                             <  2:225948/239‑1 (MQ=255)
                                                                      ccccGGTGAGTCGCCTGGGCCCCCAGCTCTTTGCCCGTCCCCGTTGCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAGCGCCTGGCGCACTGTggcggc                                                                                                                                                                            <  2:138377/250‑1 (MQ=255)
                                                                           aTGAAGCGCCTGGGCAGCAAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCACCACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAAt                                                                                                                                                                       <  2:9341/251‑1 (MQ=255)
                                                                               atcgcatgggcagcctggtatttGCCCGTCCACGTTTACACCTCAATCCACACCGACGCACTGGAGCGGCCATACCGCAAAATAGTCTGCCGTCCTTGTTCCATCTGTGGTGATTGACCACGCATATCGCCCAGCCCGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTACCGAACTGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAAt                                                                                                                                                                       >  2:201485/21‑247 (MQ=255)
                                                                               aTCGCCTGGGCCGCCAGCTTTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGAt                                                                                                                                                                   <  2:223200/251‑1 (MQ=255)
                                                                                  gCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATaa                                                                                                                                                                 >  1:140467/1‑250 (MQ=255)
                                                                                  gCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTCCCCCCTCAATCAACACCGCCGCACTGTAGCGGGCATACAGCAGAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATaaa                                                                                                                                                                <  2:98495/251‑1 (MQ=255)
                                                                                     gggcccgccAGCTCTTGGCCCGTCCCCTTTTCCCCCCCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAAt                                                                                                                                                             <  2:385993/247‑1 (MQ=255)
                                                                                                          gTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATc                                                                                                                                                                                                                                                    <  1:287959/143‑1 (MQ=255)
                                                                                                          gTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATc                                                                                                                                                                                                                                                    >  2:287959/1‑143 (MQ=255)
                                                                                                                 tgtCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGTAATAGATAAAAATTCCGGTCAGTCCGGCTTCTTCTGCCaa                                                                                                                                  <  2:360515/248‑1 (MQ=255)
                                                                                                                     ccccTCAACCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATAGCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCgg                                                                                                                             <  2:330182/251‑1 (MQ=255)
                                                                                                                           aatcaaCACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATc                                                                                                                        >  1:198015/1‑250 (MQ=255)
                                                                                                                           aatcaaCACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCa                                                                                                                       <  1:153837/251‑1 (MQ=255)
                                                                                                                               aaCACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGccc                                                                                                                   >  1:334530/1‑251 (MQ=255)
                                                                                                                                 cacCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGGTATGGCGTAACGACATGCGCGTCATATCCATCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTg                                                                                                                 <  1:196881/251‑1 (MQ=255)
                                                                                                                                  cccgccgCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATAGCCAGCGCATCGCTGAACGCCCGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGGCATTCCGGCTTCTTCTGCCAAAGCGGTAATCAGCCCTgc                                                                                                                <  2:199835/250‑1 (MQ=255)
                                                                                                                                                             aaaaTAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGc                                                                                     >  1:118995/1‑251 (MQ=255)
                                                                                                                                                                      tGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCTGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCTGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTc                                                                            >  1:332961/1‑251 (MQ=255)
                                                                                                                                                                         cGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCGGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGGTAGCGTTTAGCTCGtt                                                                         <  2:141781/251‑1 (MQ=255)
                                                                                                                                                                                gTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCAGATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGGCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCGGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTg                                                                   <  2:209052/250‑1 (MQ=255)
                                                                                                                                                                                         gtggtgATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCgcgtgcg                                                         <  1:343737/251‑1 (MQ=255)
                                                                                                                                                                                           ggtgATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCgcgtgcgt                                                        <  1:233815/250‑1 (MQ=255)
                                                                                                                                                                                                  gACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCg                                                <  1:57226/251‑1 (MQ=255)
                                                                                                                                                                                                               gCCCAGCGCGTAACGAGTACGCAGGGCGGTGCGGGTGA‑CATCGTGAGTGGTATGGCGTAACGACATGCGGGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGTCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCATGCCTTCGCCGACCACCGCTTCTGTGCCTTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGGAGCTAc                                   <  2:65850/251‑1 (MQ=255)
                                                                                                                                                                                                                   aGCACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTg                               <  1:314580/251‑1 (MQ=255)
                                                                                                                                                                                                                     cACGTAACGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGt                             >  1:137883/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaCGAGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGTCGAgg                        <  1:145618/251‑1 (MQ=255)
                                                                                                                                                                                                                              aGTACGCAGGGCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGTCGAGgcgc                     >  1:102676/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggCGTTGCGGGTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGTCGAGGCGCAAATTAAAgg           >  2:320912/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ggTGA‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGTCGAGGCGCAAATTAAAGGTTTTTTGaa  <  2:140467/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
CGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGG‑CATCGTGAGTGTTATGGCGTAACGACATGCGCGTCATATCCAGCGCATCGCTGAACGCCTGGCGCACGGTGGCGGCGGAATAGATAAAAATTCCGGTCATTCCGGCTTCTTCTGCCAAATCGGTAATCAGCCCTGCGCCGACCACCGCTTCGGTGCCGTTAGCTTTTAGCTCGTTAATCTGCCCGCGTGCGTCTTCCTCGGTAATGTAGCTACGTTGGTCGAGGCGCAAATTAAAGGTTTTTTGAA  >  NC_000913/347589‑348079

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