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,5300CT81.6% 86.7 / 13.1 38I38I (ATC→ATTrhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base T (14/17);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 9.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
*NC_0009133,966,5310CT81.6% 85.8 / 13.3 38L39L (CTG→TTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base T (14/17);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 9.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACG  >  NC_000913/3966287‑3966768
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                             
tGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagc                                                                                                                                                                                                                                         >  1:353171/1‑251 (MQ=255)
ggCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagc                                                                                                                                                                                                                                         <  2:191431/250‑1 (MQ=255)
          ccTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGGAAACAGGCATGGATTATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCa                                                                                                                                                                                                                                 <  1:438857/249‑1 (MQ=255)
               gACCAGGCAGCGAAAAGAAATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAAGCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGt                                                                                                                                                                                                                           <  1:405199/250‑1 (MQ=255)
                     gCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTCCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGaa                                                                                                                                                                                                                     >  2:457710/1‑250 (MQ=255)
                            aaaGACATGAGTCGATGACCGTAAACAGGCATGGAGGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCTCGCAAAGAGTGGCGAAGATATCtt                                                                                                                                                                                                             <  2:52566/251‑1 (MQ=255)
                             aaGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCAGTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACAGTAGTTTTGCCATTTTGAAGCTGCACGCAAAGAGTGGCGAATATATCttt                                                                                                                                                                                                            <  1:79095/251‑1 (MQ=255)
                                                  aaaCAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGc                                                                                                                                                                                                                                  >  2:291250/1‑208 (MQ=255)
                                                  aaaCAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGc                                                                                                                                                                                                                                  <  1:291250/208‑1 (MQ=255)
                                                            ggATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggat                                                                                                                                                                              <  1:240790/250‑1 (MQ=255)
                                                              atgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTGTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAgtatgga                                                                                                                                                                           <  1:95847/251‑1 (MQ=255)
                                                                    cTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGt                                                                                                                                                                     >  1:354510/1‑251 (MQ=255)
                                                                       ccATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtc                                                                                                                                                                  >  1:104544/1‑251 (MQ=255)
                                                                       ccATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCGTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtc                                                                                                                                                                  <  1:325862/251‑1 (MQ=255)
                                                                            ccATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCg                                                                                                                                                             <  2:227193/251‑1 (MQ=255)
                                                                                   caacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGc                                                                                                                                                       >  1:484130/1‑250 (MQ=255)
                                                                                               cGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTTCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTa                                                                                                                                           <  1:227166/250‑1 (MQ=255)
                                                                                                   aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAAATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGGTCCGCAGACAGCTCCTACCTcc                                                                                                                                      >  2:351139/1‑250 (MQ=255)
                                                                                                          ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc                                                                                                                               <  1:9196/251‑1 (MQ=255)
                                                                                                                 tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatga                                                                                                                         <  2:354510/250‑1 (MQ=255)
                                                                                                                     aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCt                                                                                                                    >  1:183779/1‑251 (MQ=255)
                                                                                                                               acTNTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                                                                                                          <  2:447323/251‑1 (MQ=255)
                                                                                                                                  aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAg                                                                                                       <  2:82138/251‑1 (MQ=255)
                                                                                                                                          tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                                                                                <  1:60485/250‑1 (MQ=255)
                                                                                                                                                  tAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGCTGGATTTGGTTTCCTCCGTTCCGCAGACAGCTACTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc                                                                                       >  1:427148/1‑251 (MQ=255)
                                                                                                                                                     aGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAc                                                                                    >  1:398130/1‑251 (MQ=255)
                                                                                                                                                      gcaTACGCCGGGTTCTGAGCTGCTCACTCTCGGCGAAGATATGGGGCTGGAAAACCTGGCTCGTACGCGTAAGCAGGACAATATTTTGGCCATCCTGAAGCAGCACGCAAAGAGAGGCGAAGATATCTTTGGTGATGGCGTACTGGAGCTATTGCAGGCTGGATTTGGTTTCCTCCGTTCCGCAGACAGCGCATACCTGGCCGGTACTGATGACATATAAGCTTCCCCTAACCAAATCCGCCGGTTCAAcc                                                                                   >  2:141688/3‑251 (MQ=255)
                                                                                                                                                            gCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGGGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc                                                                              >  1:395596/1‑250 (MQ=255)
                                                                                                                                                             ccGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCAc                                                                            >  1:306872/1‑251 (MQ=255)
                                                                                                                                                                         atGATCCCTCTCGGCGAAAATCTGGGGCTGGAAAACCTGGCTCATATGCGTAAGCAGCACATTATTTTTGCCATTTTGAAGCAGCACGCCAAGAGTGGCGAAGCTTTCTTTGGTGATGCCGTACCGGAGATATTGCAGGATTGATTTGGTTTCCTCCGTTCCGCAGCCAGCTCCTACCTCGCCGGTCCTGATGACATCCACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCCCTGGTGATACCAt                                                                >  2:318298/2‑251 (MQ=255)
                                                                                                                                                                                       cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                  >  2:291724/1‑251 (MQ=255)
                                                                                                                                                                                             tatGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                                            >  2:408822/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt                                 >  1:61437/1‑251 (MQ=255)
                                                                                                                                                                                                                        gcttAAGAAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGTTGGTTTTGTTTCCCTCCGTTCCGCAGACAGCTCTTACCTCCCCGGTCCTGATGACATCTACTTTTCCCCTAGCCAATCCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGTTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                  <  1:263138/247‑1 (MQ=255)
                                                                                                                                                                                                                                aGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa          >  1:167822/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ttattTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa    >  1:235454/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ttattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTACTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc   >  2:284751/1‑251 (MQ=255)
                                                                                                                                                                                                                                        attTTTGCCATTTTGAAGCAGCACGAAAAGAGTGGCGAAGATATCTTTGGTGAGGGGGTACGGGAGATATTGCAGGAGGGATTTGGTTCCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTCCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg  <  1:408822/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                             
TGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACG  >  NC_000913/3966287‑3966768

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