Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,530 2 bp→TT coding (114‑115/1260 nt) rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,5300CT98.5% 215.5 / ‑6.0 68I38I (ATC→ATT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/33);  total (34/34)
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.89e-01
*NC_0009133,966,5310CT98.5% 211.0 / NA 68L39L (CTG→TTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/33);  minor base G (0/1);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACG  >  NC_000913/3966287‑3966768
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                               
tGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAagc                                                                                                                                                                                                                                           >  1:18544/1‑251 (MQ=255)
    tACTCTCCTGTGACCAGTCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCATGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCACGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAagcagac                                                                                                                                                                                                                                       >  2:357569/1‑249 (MQ=255)
     aCTCTCCTGTGACCAGGCAGCGAAAAGACATGATTCGATGACCGGAAACAGGAATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGTGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATTGGGCTGGAAAACCTGGCTCGTATTCGTAAGCAGGACATTA‑TTTTTGCCATCTTGAagcagccc                                                                                                                                                                                                                                      <  2:465076/251‑3 (MQ=255)
     aCTCGCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGGAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTGAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACTCCGTTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAATCTGTCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCAc                                                                                                                                                                                                                                      <  1:413955/251‑1 (MQ=255)
            tgtgACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTAACCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCCGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCaa                                                                                                                                                                                                                                  >  2:247007/1‑248 (MQ=255)
            agtgACCAGGCAGCGAAACGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGACAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCCTTTTGACGCAGCACGCaa                                                                                                                                                                                                                                  >  2:13710/2‑248 (MQ=255)
               gACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGt                                                                                                                                                                                                                             <  1:430462/250‑1 (MQ=255)
                    ggCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGa                                                                                                                                                                                                                        >  1:128460/1‑250 (MQ=255)
                              aGACAGGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTg                                                                                                                                                                                                             <  2:382521/251‑1 (MQ=255)
                                        cGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGt                                                                                                                                                                                                   >  2:98505/1‑251 (MQ=255)
                                        cGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGt                                                                                                                                                                                                   <  2:581171/251‑1 (MQ=255)
                                           tGACCGTAAACAGGCATGGATGATCCTGCCATACCATTAACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGGGATGGCGT‑ACt                                                                                                                                                                                               >  1:522818/1‑251 (MQ=255)
                                                     cAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCACAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTg                                                                                                                                                                                     >  1:267654/1‑251 (MQ=255)
                                                         cATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGCTATTGCAgg                                                                                                                                                                                 >  2:181809/1‑251 (MQ=255)
                                                            ggATGATCCTGCAATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGACGG‑ACTGGAGATATTGCAgggtg                                                                                                                                                                              >  2:406854/1‑248 (MQ=255)
                                                             gatgatCCTGCCATACCATTCACACCATTAAGTTCGAGATTTACCCCACGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGCGTGGCGAAGATATCTTTGGTGCTGGCGT‑ACTGGAGATATTGCAggaggg                                                                                                                                                                             >  2:136827/1‑248 (MQ=255)
                                                               tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGGTATTGCAggctggat                                                                                                                                                                           >  2:352278/1‑251 (MQ=255)
                                                                gatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAggatggat                                                                                                                                                                           <  1:171348/250‑1 (MQ=255)
                                                                gatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGAtt                                                                                                                                                                          <  1:595562/251‑1 (MQ=255)
                                                                 atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGGttt                                                                                                                                                                         >  2:132111/1‑251 (MQ=255)
                                                                  cccTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCGTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGAGATCTTTGGTGATGGCGT‑ACGGGAGATATTGCAGGATGGATTTg                                                                                                                                                                        <  2:373902/250‑1 (MQ=255)
                                                                   ccTGCCATACCATTCACAACATTAAGTTCGAGATT‑ACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGGTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCAAGCAAAGAGGGGCGAAGATATTTTTGGTGATGGGGT‑ACTGGGGATAtttgagggatggtttggt                                                                                                                                                                      >  2:162502/1‑235 (MQ=255)
                                                                    cTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGGGGCGAAGATATCTTTGGTGATGGCGT‑AATGGAGATATTGCAGGATGGATTTGGt                                                                                                                                                                      >  1:499388/1‑251 (MQ=255)
                                                                        cATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTtcc                                                                                                                                                                  <  1:98505/251‑1 (MQ=255)
                                                                           aCCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑TCTGGAGATATTGCAGGATGGATTTGGTTtcctcc                                                                                                                                                               <  2:486584/251‑1 (MQ=255)
                                                                           aCCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTtcctc                                                                                                                                                                <  1:421252/250‑1 (MQ=255)
                                                                                           aGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCTAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTc                                                                                                                                              <  2:18544/251‑1 (MQ=255)
                                                                                               cGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTa                                                                                                                                           <  1:19776/250‑1 (MQ=255)
                                                                                                  gaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTc                                                                                                                                       <  1:346556/251‑1 (MQ=255)
                                                                                                   aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCg                                                                                                                                      >  1:368737/1‑251 (MQ=255)
                                                                                                        ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGACTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑CCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGt                                                                                                                                 >  2:461472/1‑251 (MQ=255)
                                                                                                           cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCt                                                                                                                              >  1:229849/1‑251 (MQ=255)
                                                                                                              gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                           <  1:455706/251‑1 (MQ=255)
                                                                                                              gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                           >  2:503250/1‑251 (MQ=255)
                                                                                                                  aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                        <  2:214023/250‑1 (MQ=255)
                                                                                                                  aaGAACTCACACCACTATGAAGCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCTAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGTAGATATTGCAGGATGGTTTGGGTTTCCTCCGTTCC‑GGAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                        <  1:38434/250‑1 (MQ=255)
                                                                                                                   aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGAGGGATTTTGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGCTGACAt                                                                                                                      >  1:132792/1‑251 (MQ=255)
                                                                                                                         cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGt                                                                                                                >  1:337453/1‑251 (MQ=255)
                                                                                                                         cacaCCACGATGAATCTGACCGAATTAAAGAATACGCCGGTGTCTGAACTGATCACTCTCGGCTAAAATATGTGGCTGGAAAACCTGTCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTGGAAGCAGCACTCAAAGAGTGGCGAAGATATCTGTGGTGATGGCTT‑ATGGGAGATATGGTAGGATGGAGTTGGATTCCTCCGTACT‑GCAAACATCTCCTGCCTCTGCGGTCCTGACGACTTCTAc                                                                                                                  <  2:620857/249‑1 (MQ=255)
                                                                                                                           caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                                                                                              <  2:270834/251‑1 (MQ=255)
                                                                                                                                cTATGCATCTTACCGAATTCAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGAAATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGt                                                                                                                                                                                                           >  2:620606/1‑155 (MQ=255)
                                                                                                                                cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAGCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACAGTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGt                                                                                                                                                                                                           <  1:620606/155‑1 (MQ=255)
                                                                                                                                    gAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                                                                     >  1:537786/1‑251 (MQ=255)
                                                                                                                                        cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                                                                                                 >  2:129354/1‑251 (MQ=255)
                                                                                                                                           aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATcc                                                                                               <  1:503250/250‑1 (MQ=255)
                                                                                                                                               aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGt                                                                                          >  1:118584/1‑251 (MQ=255)
                                                                                                                                                  tAAAGAATACGCCGGTTTATGGGCTGATCACTCTCGGCGAAAATATGGGGCCGGAAAACCTGGCTCGTATGCGTAAGCAGGAAATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGAGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCAGTTCC‑GCAGACAGCTCCTACCTCGCCGGCCCTGATGACATTCACGTTTCACCTAGCCAAACCAGCCGTTTc                                                                                       >  2:56632/1‑251 (MQ=255)
                                                                                                                                                          aCGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGGATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGTTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                                                                               <  2:442706/251‑1 (MQ=255)
                                                                                                                                                            gctggTTTCGGAGCTGATCACTCGCGGTGAAAAAATTGGGTTGGAAAACCTGGCTCGTATGCGTAAGCAGG‑CATTAATTTTTGCGATTGTGAGGCAGCAGGCAAAGAGTGGGAAAGATATCTTTGGTG‑TGGGTTGACTGGACATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACCCCTCCTACCCCCCCGTTCCTGATGCCAGCTACAGTTCCCCTATCCCAGTCCGCCGTTTCAACCTCCGCa                                                                             <  2:406212/248‑1 (MQ=255)
                                                                                                                                                            gCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑AATGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGAAAGCTCATACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc                                                                              >  2:497875/1‑250 (MQ=255)
                                                                                                                                                                      gAGCTGATAACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCCGCACGCACAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGACCTGATGACATCTACGTTTCCCCTAGCCCAATCCGCCGTTTCAACCTCCGCACTGGAGATAc                                                                   >  1:582350/1‑251 (MQ=255)
                                                                                                                                                                       ggctgatacctctCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAcc                                                                  <  2:175873/242‑1 (MQ=255)
                                                                                                                                                                           gATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtct                                                              <  1:129354/251‑1 (MQ=255)
                                                                                                                                                                               cctcccGGCGAAAATATGGGGCTGGAAAACCTGGCGCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                          <  2:128460/246‑1 (MQ=255)
                                                                                                                                                                                       cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                  <  1:450925/251‑1 (MQ=255)
                                                                                                                                                                                       agAAAATATGGGGCTGGAAAACCGGGCCCGTATGCGTAATCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑TCTGGAGATATTGCAGGATGGATTTGGTGTCCTCCGTTCC‑CCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCGACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAGCTCTGGTAAGATTc                                                  <  2:403617/250‑1 (MQ=255)
                                                                                                                                                                                        gAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                 >  1:323417/1‑251 (MQ=255)
                                                                                                                                                                                         aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCCCCGGTCCTGATGACCTCTCTGTTTCCCCTAGCCACATCCGCCGTTTCAACCTCCGCCCTGGTGCTACCATCTCTgagaaggttcgc                                                >  2:107813/1‑240 (MQ=255)
                                                                                                                                                                                         aaaaTATGGGGATGGAACACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGAAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGGTTCCTCCGTTAC‑GCAGACAGCTCCTACCTCGCCGGTCCAGATGACATATACGTTTCCCCTAACCAAATCCGCCGTTTCAACCTACGCACTGGTGTTACCATCTCTGGTAAGATTcgc                                                >  1:207367/1‑251 (MQ=255)
                                                                                                                                                                                           aaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGCCATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                              <  2:348615/251‑1 (MQ=255)
                                                                                                                                                                                                      ggAAAACCTGGCTCGTATGCTTAAGCAGGAAATTA‑TTTTTTCATTTTTGAAGCAGCACGAAAAGATTGGAGAAGATATTTTTGTTGATGGCGT‑ACTGGAGATATGGCAGGAGGGTTTTGGTTCCCCCCTTTCC‑GCAGCCAGCTCCTCCCTCGCCGGTCCTGATGCCATCTACTTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAGCTCTGGTAAGATTCGCCCGCCGAaagaa                                    <  1:57054/250‑1 (MQ=255)
                                                                                                                                                                                                          aaCCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                >  1:519352/1‑250 (MQ=255)
                                                                                                                                                                                                          aaCCGGGCCCGTATGTGTAAGCAGGAAATTT‑TTTTTGCCTTTTGGAAGCAGCACGAAAAGAGGGGCGAAATTATTTTTGGTGATGGCTT‑CCGGGAGATATGGCGGGATGGATTTGTTTTCCTCCTTTCC‑CCAACCAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTCCCCCTAGCAAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                               <  2:212891/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTATGCGTAAGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGGGGAGAAGATATCTTGGGTGATGGCGT‑ACTGGAGATATTGTAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATGCGCCCGCCGAAAGAAGGTGAACGCTAtttt                     <  1:162502/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGTAAGCAGGCCTTTT‑TTTTTGCCATTTTGAAGAAGCACGCAAAGAGTGGCAAAGATATCTTTGGTGATGGCTT‑CCTGGAGATATTGCATGTTGGTTTTGTTTCCCTCCTTCCC‑GCAGACACCTCCTCCCTCGCCGGTCCGGATGACATCTCCGTTTCCCCTAGCAAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                <  1:352278/251‑1 (MQ=255)
                                                                                                                                                                                                                             aGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg            >  2:564927/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCAGGACATTA‑TTTTTGCCATTTTGAAGCAGCACGAAAAGAGTGGCGAAGATATTTTTGGGGAGGGCGT‑ACTGGAGATATTGCAGGATGGTTTTGGTTTCC‑CCGTTCCCGCAGACAGCTCCTACCCCGCCGGTCCTGATGACATCTACGTTCCCCCTAGCCAAATCCGCCGTTTCACCCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg            <  1:279258/251‑1 (MQ=255)
                                                                                                                                                                                                                                        a‑ttTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaccg  >  1:207487/1‑247 (MQ=255)
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                               
TGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTA‑TTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGT‑ACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCC‑GCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACG  >  NC_000913/3966287‑3966768

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