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% 226.7 / ‑3.9 66I38I (ATC→ATTrhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/26);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01
*NC_0009133,966,5310CT98.5% 228.6 / ‑3.9 66L39L (CTG→TTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (39/26);  total (40/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAA  >  NC_000913/3966292‑3966760
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                        
cctctCCTGGGACCAGGCAGCGAAAAGACATGAGGCGATGACCTTAAACAGGCATGGATGACCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCAc                                                                                                                                                                                                                              <  2:482811/250‑1 (MQ=255)
aCTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCAc                                                                                                                                                                                                                              >  1:244876/1‑251 (MQ=255)
 ctctCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGCGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCCTTTTGAAGCAGCACg                                                                                                                                                                                                                             >  1:803912/1‑251 (MQ=255)
    cccttggACCAGCCAGCGAAAAGACAGACGGCGATGACCGTAAACAGGCATGTATTATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCaa                                                                                                                                                                                                                          <  2:607717/245‑1 (MQ=255)
          gACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTg                                                                                                                                                                                                                    <  1:220885/251‑1 (MQ=255)
           aCCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAATgg                                                                                                                                                                                                                   >  2:17147/1‑251 (MQ=255)
             cAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTCTGCGTAAGCAGGAAATTATTTTTTCCATTTTGAATCAGCACGCAAAGAGTGGCg                                                                                                                                                                                                                 >  2:727723/1‑251 (MQ=255)
             cAGGCAGCGAAAAGAAATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGACATACAATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATAAGCAGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGAGGCTGGCAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCAAGCAAAGAGTGGCg                                                                                                                                                                                                                 >  1:745993/1‑251 (MQ=255)
               ggCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGaa                                                                                                                                                                                                               <  1:510920/251‑1 (MQ=255)
                  aGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGa                                                                                                                                                                                                                >  1:587365/1‑247 (MQ=255)
                     gAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGatatt                                                                                                                                                                                                        >  2:248616/1‑250 (MQ=255)
                        aaGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTTCCATTTTGAAGCACCACGCAAAGAGGGGCGAAAATATCttt                                                                                                                                                                                                     >  2:430139/1‑251 (MQ=255)
                          gACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTgg                                                                                                                                                                                                   <  1:533515/251‑1 (MQ=255)
                                aGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAGCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATgg                                                                                                                                                                                             <  1:494606/251‑1 (MQ=255)
                                aGTCGATGACCGAAAACAGGCATGGATGATCCTGACATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGCAACTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGAAGGACATTATTTTTGCCATTTTGAAGCAGCACGAAAGGAGTGGCGAAGATATATTTGGTGTTgg                                                                                                                                                                                             >  1:721115/1‑251 (MQ=255)
                                   tCGATGACCGTAAACAGGCATGGATGACCCTGCCATACCATTCACACCATTAAGTTCGAACTTTACCCCAAGTTTAAGAACTCACACCACTATAAATCTTAACCACTTAGAGACTCCGCCGGTTTCTCAACTGATCACTCTCGGCGGAAATATGGGGCTGGAAACCCTGGCTCGTATGCGTCAGCGGGACATTATTTTTGCCCTTTTGAAGCAGCACGCAAAGAGTGGCGCAGaggtcgttaggtatggcg                                                                                                                                                                                            >  1:326960/1‑234 (MQ=255)
                                   gcGATGACCGTAAACAGGTATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCg                                                                                                                                                                                           <  2:492503/250‑1 (MQ=255)
                                    cGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGt                                                                                                                                                                                          >  1:669291/1‑251 (MQ=255)
                                            gTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgaga                                                                                                                                                                                  <  1:402000/251‑1 (MQ=255)
                                             tAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTgg                                                                                                                                                                                     >  1:745037/1‑247 (MQ=255)
                                              aaaCAGGCATGGATGAACCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAAGACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGata                                                                                                                                                                                <  2:62803/251‑1 (MQ=255)
                                                    gaagtgatgctCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTCCTGTAGATATTGCAg                                                                                                                                                                          <  2:169426/246‑1 (MQ=255)
                                                           tgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatggat                                                                                                                                                                   >  1:574178/1‑251 (MQ=255)
                                                             atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAGGATATCTTTGGTGTTGGGGTACTGGGGATATTGAAGGATGGAttt                                                                                                                                                                 >  2:715332/1‑251 (MQ=255)
                                                             atCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGAttt                                                                                                                                                                 <  1:550934/251‑1 (MQ=255)
                                                                  gccattacATTCTCAACATTAAGTTCGAGATTGACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCATGATTGATTTGGtt                                                                                                                                                             <  2:734444/243‑1 (MQ=255)
                                                                     aTACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGAAGGATGGATTTGGTTtcct                                                                                                                                                         >  2:638072/1‑251 (MQ=255)
                                                                          aTTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTg                                                                                                                                                                <  1:428566/239‑1 (MQ=255)
                                                                                aacaTTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGc                                                                                                                                                >  1:114919/1‑249 (MQ=255)
                                                                                        gTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTCCTGGAGATATTGCAGGATGGATTTGGTTTCCTCAGTTCCGCAgcacgatcc                                                                                                                                      >  2:114516/1‑243 (MQ=255)
                                                                                          tCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCtc                                                                                                                                    >  2:196404/1‑250 (MQ=255)
                                                                                           cGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCGTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCGCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTa                                                                                                                                    <  1:722673/250‑1 (MQ=255)
                                                                                               aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCg                                                                                                                               <  1:17147/251‑1 (MQ=255)
                                                                                               aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAGCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCg                                                                                                                               <  1:8018/251‑1 (MQ=255)
                                                                                                  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGcc                                                                                                                             <  1:347074/250‑1 (MQ=255)
                                                                                                    ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCCGGATGGATTTg                                                                                                                                                                >  1:43321/1‑213 (MQ=255)
                                                                                                    ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTg                                                                                                                                                                <  2:43321/213‑1 (MQ=255)
                                                                                                         aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGGTTTGGTTTCCTCCGTTCCGCAGACAGCTCCCACCTCCCCGGTCCtgt                                                                                                                     >  2:671741/1‑250 (MQ=255)
                                                                                                         aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtg                                                                                                                      >  1:791141/1‑250 (MQ=255)
                                                                                                            ttAAGAACTCACACCACTATGACTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGGTGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatgc                                                                                                                  >  2:215532/1‑250 (MQ=255)
                                                                                                            ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCCCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAATTTGAAGCAGCACGCAAAGAGTGGCGAAGATATATTTGGTGACGGCGTAATGGAGATATTGCAGGATGGATTTGGTTTCTCCCGTTCCGCAGTAAGATCCTACCTCGCCGGTCCtgatg                                                                                                                   >  2:542529/1‑250 (MQ=255)
                                                                                                              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                 <  2:194724/250‑1 (MQ=255)
                                                                                                                  aCTCAAACCACTATGAAACTTACCGAATTAACGACTACGCCGGTTTCTGCGCAGATCACTCTCAGCGAAAATATGGAGCTGGAAAACCTGGCTCGAATGCGTACGCAGGAAATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGCTGGCGCACTGGGGATATTGAAGGTAGGCTTTGGTTTCCTCGGTTACGCAGACCGCCCCTACCTCCCCgatcggaatgtcatct                                                                                                             >  2:537519/1‑235 (MQ=255)
                                                                                                                   cTCACACAACTATGAATCTTACCGAATTACAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCCCGCAAAAAGTGGCGAAGCTATCTTTGGATAAGGCGGACTGGAGATATTGACGGCAGGATTTGGTTTCCTCCGTTCCGCAGACCGCTCCTACCCAGCCGGTCCtgattgcatcctc                                                                                                           >  2:88630/1‑242 (MQ=255)
                                                                                                                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACCCTCGGCGAAAATATGGGGCTGGCAAACCTGGCTCGAATGCGTAAGCAGGACATTATTTTTCCACTTTTGAAGCAGCCCGCGAAGAGGGGCGAAGATATCTTTGGGGATGGAGTACTGGAGATATTGGAGGATGGGTTTGGTTTCCTCCGTGCCGCAGACAGCTCCTACCTCGCCGGTCCTTGTGGCATCTACGTTAccc                                                                                                    >  1:62859/1‑251 (MQ=255)
                                                                                                                           acTATGAAGCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAATCCTGGCTCGTATGCGTAAGCAGGACATGATTTGTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATGTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGTTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                                                                                                   <  1:112660/251‑1 (MQ=255)
                                                                                                                                   tCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATgg                                                                                                                                                                                             >  2:196766/1‑153 (MQ=255)
                                                                                                                                   tCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATgg                                                                                                                                                                                             <  1:196766/153‑1 (MQ=255)
                                                                                                                                          gAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGTATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAACCAAATccgccg                                                                                    >  1:474163/1‑251 (MQ=255)
                                                                                                                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgag                                                                                                                                                                                   <  1:179905/155‑1 (MQ=255)
                                                                                                                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGgag                                                                                                                                                                                   >  2:179905/1‑155 (MQ=255)
                                                                                                                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGt                                                                                   >  1:436711/1‑251 (MQ=255)
                                                                                                                                                  gAATACGCCGGTTTCTGAGCTGCTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCCGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGCCAGCTCCTACCTCGCCGGTCCTGATGCCCTCTACGTTTCCCCTAGCCAAATCCGCCGttttcacc                                                                            >  2:327458/1‑246 (MQ=255)
                                                                                                                                                  gAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTACGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                                                            >  2:314137/1‑251 (MQ=255)
                                                                                                                                                           ggTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCATGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTg                                                                   >  1:22386/1‑251 (MQ=255)
                                                                                                                                                           ggTTTCTGAGCGGATCACTCTCGGCGAAAATATGGGGCGGGAAACCCTGGCTCGTAGGCGTAAGCAGGACATTATTTTTGCCATTTGGAAGCAGCACGCAAAGAGTGGCGAAGATAGCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCATACAGCTCTTCCCTCGCCGGTCCTTAGGACAGCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTg                                                                   <  1:715332/251‑1 (MQ=255)
                                                                                                                                                                                    gAAAATATGGGGCTGGAAACCCTGGCTCCTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGAGGCCAAGATATCTTTGGTTATGGCGTACTGGAGATATTGCCGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCCAATCCGCCGTTTCAACCTCCGCCCTGGTGATACCCTCTCTGGCAAGAtttc                                         >  2:416004/1‑249 (MQ=255)
                                                                                                                                                                                                     aacccTGGCTGGTAGGCGTAAGCGGGACATTATTTTTGCTTTTTGAAAGCAGCAGGCAAAGAGGGGCGAAGATTTCTTGGGTGAGGGCGTACGGGAGATATTGAGGGATGTATTTGGTTCCCTCCGTCCCGCAGACAGCTCCTCCCTCGCCGTTCCTTAGGACATTTACGTTCCCCCTAGCAAAATCCGCCGTTTCAACCCCCGAACGGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                        <  1:114516/248‑1 (MQ=255)
                                                                                                                                                                                                         cTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCCGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGCACg                    >  2:13831/1‑251 (MQ=255)
                                                                                                                                                                                                          tGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGc                   >  2:599632/1‑251 (MQ=255)
                                                                                                                                                                                                            gctggttggCGTAAGCAGGATATTAGTTTTGCCTTTTCGAATCAGCACGCAAAGAGTGGCAATATTATCTTGGGTGATGGCTGGCTGAATATATTGCAGGATGGATTGGTTTTCCTCCGTTCCGCAGACAGCTCTTACCCCGCCGGTCCTGATGACACCTACTGTTCCCTTAGCCAAACCCGCCGTTTCACCCCCCGCACTTTTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAACGTGAACGCTa                 <  1:676039/243‑1 (MQ=255)
                                                                                                                                                                                                            gCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGCCAGCTCCTACCTCGCCGGTCCTGCTGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGACGGGGAACGCt                  >  2:470783/1‑250 (MQ=255)
                                                                                                                                                                                                               cttATGCGAAAGCGGGACATTTTTTTTCCCTTTTTAAAGCGGCACGAAAAGAGGGGGGAAGATATTTTTGGTGATGGCGTACGGGAGATATGGCAGGAGGGATTTGGTTTCCCCCTTCCCGCAGACAGCCCCTCCCCCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCCCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt                <  1:248616/247‑1 (MQ=255)
                                                                                                                                                                                                                  aTGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGGCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGGAAGAAGGGTAACGCTAtttttc           >  2:176891/1‑249 (MQ=255)
                                                                                                                                                                                                                         aGCAGGACATTATTTTTGCCATATTGAAGCAGCACGCAAAGAGTGGCGAGGATATCTTTGGGGATGGCGTACTGGAGATAGTGCTGGGTGGATTTGGTTTCCTCCGTTCCGCAGACAACTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCCAATCCGCCGTTTCAACCTCCGCACTGGGGTACCCATCTCTGGAAACATTCGACAGCCGAAAGAAAGTGAACGCTAtttttctcatacg    >  2:591576/1‑242 (MQ=255)
                                                                                                                                                                                                                             ggACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa  >  1:729097/1‑249 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                        
ACTCTCCTGTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAA  >  NC_000913/3966292‑3966760

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