Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,796 A→G intergenic (+1239/‑405) ilvX → / → ilvM uncharacterized protein IlvX/acetolactate synthase II subunit IlvM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,7960AG97.0% 173.5 / ‑3.3 66intergenic (+1239/‑405)ilvX/ilvMuncharacterized protein IlvX/acetolactate synthase II subunit IlvM
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (24/40);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

ACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGG  >  NC_000913/3951553‑3952041
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
acacTGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTTCGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCTCATCGCCCACACTCGCCCGGAAAATTTCATCTCCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTc                                                                                                                                                                                                                                                   <  1:373605/250‑1 (MQ=255)
acacTGCGCGCAGCTGCGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGGGCGCAAGTcg                                                                                                                                                                                                                                                  >  1:62104/1‑251 (MQ=255)
               gCGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGAt                                                                                                                                                                                                                                    >  1:191923/1‑250 (MQ=255)
                cGTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATAc                                                                                                                                                                                                                                  <  1:213212/251‑1 (MQ=255)
                 gTGATGAACATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCCGTTGGCGCGCAAGTCGCGCGACCAAACGATAcc                                                                                                                                                                                                                                 >  1:315695/1‑251 (MQ=255)
                   gatgaACATGCCTGTCGGTACGACCAGCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCTGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGt                                                                                                                                                                                                                               <  2:14752/251‑1 (MQ=255)
                        aCATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCt                                                                                                                                                                                                                          <  1:296083/251‑1 (MQ=255)
                        aCATGCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATAACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGGTGGCGCGCAAGTCGCGCGACCCAAccgcaccgtcgtct                                                                                                                                                                                                                          >  2:20118/1‑238 (MQ=255)
                         cctgccTGGCGTTACGACCAGCCCGGTGACGCTATCTACGCGCCGTAGTTGTTAAAACAACTGTCAGTACGTAAACCTTCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATTTGGTCTGCGCAACACATCGCCCCCACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGTGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTg                                                                                                                                                                                                                        <  2:177726/247‑1 (MQ=255)
                                ggCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGGTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTcc                                                                                                                                                                                                                 <  2:288172/251‑1 (MQ=255)
                                        gACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACgg                                                                                                                                                                                                         <  1:28529/251‑1 (MQ=255)
                                         cccATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGc                                                                                                                                                                                                        <  2:149572/250‑1 (MQ=255)
                                         aCCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCGGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGc                                                                                                                                                                                                        <  2:69403/251‑1 (MQ=255)
                                                  gTGACGCTATCTACGCTCCGTTGTTGTTAAAACCACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCCCATCTCCCCCACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTGGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCatg                                                                                                                                                                                               <  1:470692/251‑1 (MQ=255)
                                                        cTATCTACGCGCCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGTCTGTATCTCCCGTGACGGCTCTTTCctgatgca                                                                                                                                                                                          >  2:63307/1‑248 (MQ=255)
                                                                  gcCGTTGTTGTTAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGGTGGCGCGCAAGTCGCGCGACCAAACGATAACGGCGGCTGTATCTACGGGGACGGCTCTTTCATGATGAATGAGGAagagg                                                                                                                                                                               >  2:52726/1‑250 (MQ=255)
                                                                          tgttAAAACAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCAc                                                                                                                                                                        <  1:99127/250‑1 (MQ=255)
                                                                           ggtAAAACTACTGTCAGATTGTAACCCTGCGGATTGCGTCGTGACGACATATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACCCTGGCCCGGAGAATTTCATCCCCGCCAGCGGCTTAGTCGCCATGGGGGTTGGGTGACCGGCGTTGGGTGGCGCGCAAGGCGCGCGACCAAACTATATCGGTGTCGGTATCTCCGGTGACGGCTCTTTCGTGATGAATGTGCAAGAGCTGGGCAcc                                                                                                                                                                       <  1:188275/248‑1 (MQ=255)
                                                                               aaaCAACTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGGGCACAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGGACCGTaa                                                                                                                                                                   >  1:202389/1‑250 (MQ=255)
                                                                                aacaacTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCCCACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTcc                                                                                                                                                                                                                 <  2:451648/203‑1 (MQ=255)
                                                                                aacaacTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTCGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGGAAGTCGCGCGACCAAACGACACCGGCGTCCGTATCTcc                                                                                                                                                                                                                 >  1:451648/1‑203 (MQ=255)
                                                                                    acTGTCAGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGAGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCTGCGGTTGGCGCGCAAGGAGCGCGACCAAACCATACCGTCGTCTGTCTCTCCGGTCACGGCCCTTTCCTGATGAATGGTCTAGAGCTggggcccggataacgcc                                                                                                                                                             >  2:86215/1‑237 (MQ=255)
                                                                                       gtcagATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTTGCGCGCCAGACGCGCGACCCAACCACACCCACGTCTTGTTATCCGGCGCCGGCTCTTTCATGATGACTGTGCAAGTGCTGGGccacgcacaaagtccgc                                                                                                                                                          >  2:94338/5‑235 (MQ=255)
                                                                                       gtcagATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGCGCACCGTAAAACGCAAGc                                                                                                                                                         <  2:409858/247‑1 (MQ=255)
                                                                                              cGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCCAGTCGCGCGACCCAACGATACCGGCGGCTGTATCTCCGGGGACGGCTCTTTCAATATGAATGTGCAAGAGCTGGGCACCGTAAAACGccaacaggtacc                                                                                                                                                  >  2:467456/1‑240 (MQ=255)
                                                                                                    ccTGGGGATTGCTGCGTGACCCCAGAGGTGGGGCAGCGCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCTGCTTAGGCACCATTGGTTGTGGTTTACCGGCGGCGGTTGGCGCGCAATTCGCGCGACCAAACGATACCGTCGTCTGTAGCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCGGGGCACCGTAAAACGCAAGCAGTTACCGTTGat                                                                                                                                            <  2:78858/251‑2 (MQ=255)
                                                                                                         ggATTGCGTCGTGACCACAGATGCGGGGCAGCACCAGATGCGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAAAGATACCGTCGTCTGTTTCCCCGGTGACAGCCCTTTCATGATGAATGTGCAACAGCTGGGCACCGTAAAACGCAAGCAGTTACCCTTGTAAAt                                                                                                                                         >  1:175486/1‑249 (MQ=255)
                                                                                                             tGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCt                                                                                                                                    <  1:114142/250‑1 (MQ=255)
                                                                                                              gCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTAACCGTTGAAAATCGTCtt                                                                                                                                   <  2:401999/250‑1 (MQ=255)
                                                                                                                gtcgtGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACt                                                                                                                                >  1:195108/1‑251 (MQ=255)
                                                                                                                gtcgtGAACACCGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCCCAGCGCCCACTCTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTGTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCGTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAAGCGTCTTACt                                                                                                                                <  2:266724/251‑1 (MQ=255)
                                                                                                                       ccACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAAc                                                                                                                         <  1:252089/251‑1 (MQ=255)
                                                                                                                           aGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAAc                                                                                                                     <  1:381738/251‑1 (MQ=255)
                                                                                                                             aggtgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCGCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACgg                                                                                                                   <  2:202389/249‑1 (MQ=255)
                                                                                                                              tgtgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCACGCCAATCGCGCCACCAAAAGATACCGTCGTCTGTATCTCCGGTCACCGCTCTTTCATGATGAATGGGCAAGAGCTGGGCAACGTAAAACGACAACAGTTATCATTGAAAATCGACTTACTCCATACCCAACggg                                                                                                                  >  2:346111/1‑250 (MQ=255)
                                                                                                                                 ggggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAg                                                                                                               <  1:63307/251‑1 (MQ=255)
                                                                                                                                  gggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAgg                                                                                                              <  1:148946/251‑1 (MQ=255)
                                                                                                                                   ggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACCGTTatgg                                                                                                             >  2:149849/1‑248 (MQ=255)
                                                                                                                                     cagcaCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGAt                                                                                                           <  2:335708/251‑1 (MQ=255)
                                                                                                                                         aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGtt                                                                                                       <  1:238240/251‑1 (MQ=255)
                                                                                                                                         aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAAAGATAACGTCGTCTGTATCTCCGGGGACGGCTCTTTCATGATTAATGTGCAAGAGCTGGGCCCCGTAAAACGCAAGCAGTTACCGGTGAAAATCGTCCTACTCGATAAACAACGGGTAGGGAttggt                                                                                                       >  2:34521/1‑247 (MQ=255)
                                                                                                                                             gATGTGGGCTGCGCAGCACATCGCCCCCACTNGCCCGGAAAATTTCATCACCTCCAGCGGATTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTCCTCGATAGCCAACGGTTAGGGATGGTTCGAc                                                                                                   <  2:227896/251‑1 (MQ=255)
                                                                                                                                              gtgtgGGCTGCGCAGCACAGCGCCCACACGCNCCCGTAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGATTTTGGTGTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGGCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCTAGAGCTGGGCACCGTAAAACGCAAGCAGTTAGCGTTGAAAATCGTCTTACTCCATAACCAACGGTTATGGATGGTTCGACa                                                                                                  <  2:234096/250‑1 (MQ=255)
                                                                                                                                                           agctcaTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTg                                                                                     <  1:241764/247‑1 (MQ=255)
                                                                                                                                                               caTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAATCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCCGTTAACGTTGAAAATCGTCTTACTCGAAAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttt                                                                                 >  1:96026/1‑251 (MQ=255)
                                                                                                                                                                aTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGAACAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTGCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGttttt                                                                                <  2:459278/251‑1 (MQ=255)
                                                                                                                                                                  cGCCCACACTCGCCCGGAAAATTTCATCACNCCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttttc                                                                               <  2:224874/250‑2 (MQ=255)
                                                                                                                                                                            cGCCCGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGa                                                                      <  2:195108/249‑1 (MQ=255)
                                                                                                                                                                               ccGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCACGTCGCGCGAGCCAACGCTAACGTCGTCTGTATCTCCGGGGAAGGCTCCTTAATGATGAATGTGCAAGAGCCGGGCACCGTAAAACGCAAGCAGCTACCGTTGAAAATCGTCTTACTCGATAACCAACGGGTAGGGGTGGGTCGACAANGGCAGAAACTGGTTTTTCAGGaaaggtct                                                                   >  2:234853/1‑243 (MQ=255)
                                                                                                                                                                                cGGAACATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGCTACCGTCGTCTGTATCTCCGGTGACGGCTCTTTTATGTTGAATGTGCAAGAGCTTGGCACCGTAACACGCCAGCCGTTACCGTTGAAAATCGTATTAATCGATCAACAACGGTTCGGGATGGTTTGACCAANGCCGCAACTGTTTTTTCCAGAACGATaata                                                                 >  2:222677/1‑247 (MQ=255)
                                                                                                                                                                                cGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGc                                                                >  1:54549/1‑251 (MQ=255)
                                                                                                                                                                                  gAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGGCGCGCGACCAAAAGATACCGTCGTCTGTATTTCCCGGGACGGCTCTTTCATGTTGAATGTGCAAGAGCTGGGCACCGTAAAAAGCAAGCAGTTACCGTTGAAAAACGTCTTACTCGATAAACAACGGTTTGGGATGGGTTGTCAATGGCCGCAACTTTTTTTCCAGGAACGATATAGcc                                                               >  2:126472/1‑249 (MQ=255)
                                                                                                                                                                                     aaTTTCAGCACCTCCAGCGGCTTAGGCACCATGGGGTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCTACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAac                                                           <  1:234853/251‑1 (MQ=255)
                                                                                                                                                                                             caCCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTa                                                   <  1:248868/251‑1 (MQ=255)
                                                                                                                                                                                               ccTCCAGCGGCTTAGGCACCATGGGTTTTGGTTGAACGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCATCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTAcc                                                 <  1:349136/251‑2 (MQ=255)
                                                                                                                                                                                                     gCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTGGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGAGGGTTCGCCAATGGCAGCAACTGTTTTGCCAGGAACGATATAGCGAAACCACCCTTACCGATAAc                                           <  2:172151/251‑1 (MQ=255)
                                                                                                                                                                                                        gcttAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAAcccc                                        <  2:176540/249‑1 (MQ=255)
                                                                                                                                                                                                                 ccATGGTTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGAGGAATGTGCAAGAGCTGGGCACCGTAAACCGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTc                               <  2:277833/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTGGGTTTTGGTTTCCCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCAt                             <  2:79566/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTGGGGTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATCCCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTGCCTCAt                             <  2:189470/251‑1 (MQ=255)
                                                                                                                                                                                                                     gggTTTTGGTTTACCGGCGGCGGTTTGCGCGCAAGTCCCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACCCAAGCAGTTACCGGTGAAAATCGTCTTAATCGATAACCAACGGGTAGGGCTGGTTCGACAATGGCAGCAAATGTTTTTCCAGGAACGAAATAACGGAACCACCCTTACCGCTAAACCCGATTTCCTCaagt                           >  2:100892/1‑248 (MQ=255)
                                                                                                                                                                                                                       gTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGtt                          >  1:376604/1‑250 (MQ=255)
                                                                                                                                                                                                                              tttACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCCAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGCTATAGCGAAAACAACCTTACNCATAACCCCGATTTCCTCCTGtttaccaccc                  >  2:226047/1‑243 (MQ=255)
                                                                                                                                                                                                                                          ggTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATcc      <  1:40625/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTgg  <  1:31058/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTgg  >  1:150998/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
ACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGG  >  NC_000913/3951553‑3952041

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: