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→ATTrhotranscription 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

TGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACG  >  NC_000913/3966287‑3966768
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                               
tGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagc                                                                                                                                                                                                                                           >  1:18544/1‑251 (MQ=255)
    tACTCTCCTGTGACCAGTCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCATGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCACGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAagcagac                                                                                                                                                                                                                                       >  2:357569/1‑249 (MQ=255)
     aCTCTCCTGTGACCAGGCAGCGAAAAGACATGATTCGATGACCGGAAACAGGAATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGTGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATTGGGCTGGAAAACCTGGCTCGTATTCGTAAGCAGGACATTATTTTTGCCATCTTGAagcagccc                                                                                                                                                                                                                                      <  2:465076/251‑3 (MQ=255)
     aCTCGCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGGAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTGAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACTCCGTTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAATCTGTCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCAc                                                                                                                                                                                                                                      <  1:413955/251‑1 (MQ=255)
            tgtgACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTAACCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCCGGACATTATTTTTGCCATTTTGAAGCAGCACGCaa                                                                                                                                                                                                                                  >  2:247007/1‑248 (MQ=255)
            agtgACCAGGCAGCGAAACGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGACAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCCTTTTGACGCAGCACGCaa                                                                                                                                                                                                                                  >  2:13710/2‑248 (MQ=255)
               gACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGt                                                                                                                                                                                                                             <  1:430462/250‑1 (MQ=255)
                    ggCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGa                                                                                                                                                                                                                        >  1:128460/1‑250 (MQ=255)
                              aGACAGGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTg                                                                                                                                                                                                             <  2:382521/251‑1 (MQ=255)
                                        cGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGt                                                                                                                                                                                                   >  2:98505/1‑251 (MQ=255)
                                        cGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGt                                                                                                                                                                                                   <  2:581171/251‑1 (MQ=255)
                                           tGACCGTAAACAGGCATGGATGATCCTGCCATACCATTAACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGGGATGGCGTACt                                                                                                                                                                                               >  1:522818/1‑251 (MQ=255)
                                                     cAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCACAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTg                                                                                                                                                                                     >  1:267654/1‑251 (MQ=255)
                                                         cATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGCTATTGCAgg                                                                                                                                                                                 >  2:181809/1‑251 (MQ=255)
                                                            ggATGATCCTGCAATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGACGGACTGGAGATATTGCAgggtg                                                                                                                                                                              >  2:406854/1‑248 (MQ=255)
                                                             gatgatCCTGCCATACCATTCACACCATTAAGTTCGAGATTTACCCCACGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGCGTGGCGAAGATATCTTTGGTGCTGGCGTACTGGAGATATTGCAggaggg                                                                                                                                                                             >  2:136827/1‑248 (MQ=255)
                                                               tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGGTATTGCAggctggat                                                                                                                                                                           >  2:352278/1‑251 (MQ=255)
                                                                gatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatggat                                                                                                                                                                           <  1:171348/250‑1 (MQ=255)
                                                                gatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAtt                                                                                                                                                                          <  1:595562/251‑1 (MQ=255)
                                                                 atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGGttt                                                                                                                                                                         >  2:132111/1‑251 (MQ=255)
                                                                  cccTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCGTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGAGATCTTTGGTGATGGCGTACGGGAGATATTGCAGGATGGATTTg                                                                                                                                                                        <  2:373902/250‑1 (MQ=255)
                                                                   ccTGCCATACCATTCACAACATTAAGTTCGAGATTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGGTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCAAGCAAAGAGGGGCGAAGATATTTTTGGTGATGGGGTACTGGGGATAtttgagggatggtttggt                                                                                                                                                                      >  2:162502/1‑235 (MQ=255)
                                                                    cTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGGGGCGAAGATATCTTTGGTGATGGCGTAATGGAGATATTGCAGGATGGATTTGGt                                                                                                                                                                      >  1:499388/1‑251 (MQ=255)
                                                                        cATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcc                                                                                                                                                                  <  1:98505/251‑1 (MQ=255)
                                                                           aCCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTTCTGGAGATATTGCAGGATGGATTTGGTTtcctcc                                                                                                                                                               <  2:486584/251‑1 (MQ=255)
                                                                           aCCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTtcctc                                                                                                                                                                <  1:421252/250‑1 (MQ=255)
                                                                                           aGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCTAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTc                                                                                                                                              <  2:18544/251‑1 (MQ=255)
                                                                                               cGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTa                                                                                                                                           <  1:19776/250‑1 (MQ=255)
                                                                                                  gaTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTc                                                                                                                                       <  1:346556/251‑1 (MQ=255)
                                                                                                   aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCg                                                                                                                                      >  1:368737/1‑251 (MQ=255)
                                                                                                        ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGACTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTCCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGt                                                                                                                                 >  2:461472/1‑251 (MQ=255)
                                                                                                           cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCt                                                                                                                              >  1:229849/1‑251 (MQ=255)
                                                                                                              gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                           <  1:455706/251‑1 (MQ=255)
                                                                                                              gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                                           >  2:503250/1‑251 (MQ=255)
                                                                                                                  aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                        <  2:214023/250‑1 (MQ=255)
                                                                                                                  aaGAACTCACACCACTATGAAGCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCTAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGTAGATATTGCAGGATGGTTTGGGTTTCCTCCGTTCCGGAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                        <  1:38434/250‑1 (MQ=255)
                                                                                                                   aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGAGGGATTTTGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGCTGACAt                                                                                                                      >  1:132792/1‑251 (MQ=255)
                                                                                                                         cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGt                                                                                                                >  1:337453/1‑251 (MQ=255)
                                                                                                                         cacaCCACGATGAATCTGACCGAATTAAAGAATACGCCGGTGTCTGAACTGATCACTCTCGGCTAAAATATGTGGCTGGAAAACCTGTCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTGGAAGCAGCACTCAAAGAGTGGCGAAGATATCTGTGGTGATGGCTTATGGGAGATATGGTAGGATGGAGTTGGATTCCTCCGTACTGCAAACATCTCCTGCCTCTGCGGTCCTGACGACTTCTAc                                                                                                                  <  2:620857/249‑1 (MQ=255)
                                                                                                                           caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                                                                                              <  2:270834/251‑1 (MQ=255)
                                                                                                                                cTATGCATCTTACCGAATTCAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGAAATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGt                                                                                                                                                                                                           >  2:620606/1‑155 (MQ=255)
                                                                                                                                cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAGCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACAGTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGt                                                                                                                                                                                                           <  1:620606/155‑1 (MQ=255)
                                                                                                                                    gAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                                                                     >  1:537786/1‑251 (MQ=255)
                                                                                                                                        cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAAt                                                                                                 >  2:129354/1‑251 (MQ=255)
                                                                                                                                           aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATcc                                                                                               <  1:503250/250‑1 (MQ=255)
                                                                                                                                               aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGt                                                                                          >  1:118584/1‑251 (MQ=255)
                                                                                                                                                  tAAAGAATACGCCGGTTTATGGGCTGATCACTCTCGGCGAAAATATGGGGCCGGAAAACCTGGCTCGTATGCGTAAGCAGGAAATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGAGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCAGTTCCGCAGACAGCTCCTACCTCGCCGGCCCTGATGACATTCACGTTTCACCTAGCCAAACCAGCCGTTTc                                                                                       >  2:56632/1‑251 (MQ=255)
                                                                                                                                                          aCGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGGATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGTTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                                                                               <  2:442706/251‑1 (MQ=255)
                                                                                                                                                            gctggTTTCGGAGCTGATCACTCGCGGTGAAAAAATTGGGTTGGAAAACCTGGCTCGTATGCGTAAGCAGGCATTAATTTTTGCGATTGTGAGGCAGCAGGCAAAGAGTGGGAAAGATATCTTTGGTGTGGGTTGACTGGACATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACCCCTCCTACCCCCCCGTTCCTGATGCCAGCTACAGTTCCCCTATCCCAGTCCGCCGTTTCAACCTCCGCa                                                                             <  2:406212/248‑1 (MQ=255)
                                                                                                                                                            gCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTAATGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGAAAGCTCATACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGc                                                                              >  2:497875/1‑250 (MQ=255)
                                                                                                                                                                      gAGCTGATAACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCCGCACGCACAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGACCTGATGACATCTACGTTTCCCCTAGCCCAATCCGCCGTTTCAACCTCCGCACTGGAGATAc                                                                   >  1:582350/1‑251 (MQ=255)
                                                                                                                                                                       ggctgatacctctCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAcc                                                                  <  2:175873/242‑1 (MQ=255)
                                                                                                                                                                           gATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtct                                                              <  1:129354/251‑1 (MQ=255)
                                                                                                                                                                               cctcccGGCGAAAATATGGGGCTGGAAAACCTGGCGCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                          <  2:128460/246‑1 (MQ=255)
                                                                                                                                                                                       cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                  <  1:450925/251‑1 (MQ=255)
                                                                                                                                                                                       agAAAATATGGGGCTGGAAAACCGGGCCCGTATGCGTAATCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTTCTGGAGATATTGCAGGATGGATTTGGTGTCCTCCGTTCCCCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCGACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAGCTCTGGTAAGATTc                                                  <  2:403617/250‑1 (MQ=255)
                                                                                                                                                                                        gAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                 >  1:323417/1‑251 (MQ=255)
                                                                                                                                                                                         aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCCCCGGTCCTGATGACCTCTCTGTTTCCCCTAGCCACATCCGCCGTTTCAACCTCCGCCCTGGTGCTACCATCTCTgagaaggttcgc                                                >  2:107813/1‑240 (MQ=255)
                                                                                                                                                                                         aaaaTATGGGGATGGAACACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGAAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGGTTCCTCCGTTACGCAGACAGCTCCTACCTCGCCGGTCCAGATGACATATACGTTTCCCCTAACCAAATCCGCCGTTTCAACCTACGCACTGGTGTTACCATCTCTGGTAAGATTcgc                                                >  1:207367/1‑251 (MQ=255)
                                                                                                                                                                                           aaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGCCATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                              <  2:348615/251‑1 (MQ=255)
                                                                                                                                                                                                      ggAAAACCTGGCTCGTATGCTTAAGCAGGAAATTATTTTTTCATTTTTGAAGCAGCACGAAAAGATTGGAGAAGATATTTTTGTTGATGGCGTACTGGAGATATGGCAGGAGGGTTTTGGTTCCCCCCTTTCCGCAGCCAGCTCCTCCCTCGCCGGTCCTGATGCCATCTACTTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAGCTCTGGTAAGATTCGCCCGCCGAaagaa                                    <  1:57054/250‑1 (MQ=255)
                                                                                                                                                                                                          aaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                >  1:519352/1‑250 (MQ=255)
                                                                                                                                                                                                          aaCCGGGCCCGTATGTGTAAGCAGGAAATTTTTTTTGCCTTTTGGAAGCAGCACGAAAAGAGGGGCGAAATTATTTTTGGTGATGGCTTCCGGGAGATATGGCGGGATGGATTTGTTTTCCTCCTTTCCCCAACCAGCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTCCCCCTAGCAAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                               <  2:212891/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGGGGAGAAGATATCTTGGGTGATGGCGTACTGGAGATATTGTAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATGCGCCCGCCGAAAGAAGGTGAACGCTAtttt                     <  1:162502/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGTAAGCAGGCCTTTTTTTTTGCCATTTTGAAGAAGCACGCAAAGAGTGGCAAAGATATCTTTGGTGATGGCTTCCTGGAGATATTGCATGTTGGTTTTGTTTCCCTCCTTCCCGCAGACACCTCCTCCCTCGCCGGTCCGGATGACATCTCCGTTTCCCCTAGCAAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                <  1:352278/251‑1 (MQ=255)
                                                                                                                                                                                                                             aGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg            >  2:564927/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGAAAAGAGTGGCGAAGATATTTTTGGGGAGGGCGTACTGGAGATATTGCAGGATGGTTTTGGTTTCCCCGTTCCCGCAGACAGCTCCTACCCCGCCGGTCCTGATGACATCTACGTTCCCCCTAGCCAAATCCGCCGTTTCACCCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg            <  1:279258/251‑1 (MQ=255)
                                                                                                                                                                                                                                        attTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaccg  >  1:207487/1‑247 (MQ=255)
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                               
TGCGTACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACG  >  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: