Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,877 C→A 9.1% T154N (ACT→AAT)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,8770CA 9.1% 341.5 / 17.9 125T154N (ACT→AAT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (51/62);  new base A (4/7);  total (56/69)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 127 total aligned reads displayed.

AAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGC  >  NC_000913/3966729‑3967026
                                                                                                                                                    |                                                                                                                                                     
aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                                      <  1:2006312/150‑1 (MQ=255)
aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                                      >  2:462756/1‑150 (MQ=255)
 aagaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTg                                                                                                                                                     >  1:1250738/1‑150 (MQ=255)
 aagaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATg                                                                                                                                                     <  2:429146/150‑1 (MQ=255)
      ggTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGAt                                                                                                                                                <  1:792078/150‑1 (MQ=255)
      ggTGAACGCTAATTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAACGCTCGTCTGCGTATGTAACGTGGTAACGGTTCTACTGAAGAt                                                                                                                                                <  1:1210404/150‑1 (MQ=255)
           aCGCTATTTTGCGCTGCTGAAAGTTAAGGAAGTTAACTTCGACAAACCTGAAAACGCCCTCAACAAAATCCTCTTTGAGAACTTAACCCCGCGGCACGCAAACTCTCGTCTGCGTATGGAACGTGGAAACGGTTCTACTGAAGATTTAAc                                                                                                                                           <  1:1743746/150‑1 (MQ=255)
           aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                                           >  1:472245/1‑150 (MQ=255)
           aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                                           <  1:1907112/150‑1 (MQ=255)
           aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                                           >  1:802676/1‑150 (MQ=255)
           aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                                           <  1:1338587/150‑1 (MQ=255)
             gCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg                                                                                                                                         >  2:496463/1‑150 (MQ=255)
                    tGCGCTGCTGAAAGTTAACGAAGTTAAATTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCAGCACGCAAACTCTCGTCAGCGTATGGAACATGGTAACGCTTCTACTGAAGAATTAACTGCTCGCGt                                                                                                                                  >  1:1318363/1‑150 (MQ=255)
                     gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTTACGGTTCTACTGAAGATTTAACTGCTCGCGTa                                                                                                                                 >  2:223883/1‑150 (MQ=255)
                     gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTa                                                                                                                                 <  2:1079092/150‑1 (MQ=255)
                         tgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTgg                                                                                                                             <  1:1446346/150‑1 (MQ=255)
                          gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAATGCTCGCGTACTGGa                                                                                                                            >  2:617731/1‑150 (MQ=255)
                          gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAAGGTTCTAATGAAGATTTAACTGCTCGCGTAATGGa                                                                                                                            >  1:1177094/1‑150 (MQ=255)
                            tgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATc                                                                                                                          <  2:199542/150‑1 (MQ=255)
                            tgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGGGGTACAGGTTCTATTGAAGATTTATCTGCTCGCGTACTGGATc                                                                                                                          >  2:1335875/1‑150 (MQ=255)
                            tcgaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAACCTCTTTGAGAACTTAACCCCGCTGCACGCACACTCGCGTCTGCGAATGGAACATGGATACGGTTCTACTGAAGATTTAACTGCTCGCGTCATGGATc                                                                                                                          >  2:1831960/4‑150 (MQ=255)
                             gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                         <  2:674647/150‑1 (MQ=255)
                             gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                         >  2:1922757/1‑150 (MQ=255)
                             gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                         >  2:1161611/1‑150 (MQ=255)
                             cgaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                                         >  1:479532/3‑150 (MQ=255)
                               aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                                                                       >  1:1853772/1‑150 (MQ=255)
                                 gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                                                                                     <  1:835007/150‑1 (MQ=255)
                                 gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                                                                                     <  2:251382/150‑1 (MQ=255)
                                   tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATc                                                                                                                   <  1:174161/150‑1 (MQ=255)
                                     aCGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAc                                                                                                                 >  2:271767/1‑150 (MQ=255)
                                      cGAAGTTAACTTAGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCAcc                                                                                                                <  1:1814238/150‑1 (MQ=255)
                                        aaGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGAATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATAACCTa                                                                                                              >  2:942062/1‑150 (MQ=255)
                                          gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                                                                            <  1:225283/150‑1 (MQ=255)
                                            tAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCgg                                                                                                          >  1:1059702/1‑150 (MQ=255)
                                            tAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCgg                                                                                                          >  2:1236826/1‑150 (MQ=255)
                                              aCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                                        >  1:445359/1‑150 (MQ=255)
                                              aCTTCGACAAACCTGAAAACGCCCGCAACAAAATCATCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTATGCGTATGGAACGTGGTAACGGTTTTAATGAAGATTTAAGTGATCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                                        <  2:399416/150‑1 (MQ=255)
                                               cTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                                                                                       <  1:658901/150‑1 (MQ=255)
                                                ttCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACGCTCGTCTGCGTATGGACCATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGt                                                                                                      <  1:1156455/150‑1 (MQ=255)
                                                ttCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACACAAACTCTCGTCTGCGTATGGAACGGGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCgg                                                                                                      >  2:1040400/1‑149 (MQ=255)
                                                  cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                                                                                    <  2:1345527/150‑1 (MQ=255)
                                                  cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGTTTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                                                                                    <  1:31073/150‑1 (MQ=255)
                                                  cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCACGTCTGCGTATGGAACATGGTAACGGTTCTACTGAATATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                                                                                    <  1:601100/150‑1 (MQ=255)
                                                    aCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTc                                                                                                  <  2:526567/150‑1 (MQ=255)
                                                          cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTg                                                                                            >  1:1941202/1‑150 (MQ=255)
                                                          cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCGATCGGTCGTGGGCAGCGTg                                                                                            >  1:677344/1‑150 (MQ=255)
                                                             aaaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTc                                                                                         <  1:1416344/150‑1 (MQ=255)
                                                                   cccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtg                                                                                   >  2:1082699/1‑150 (MQ=255)
                                                                    ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                                  >  1:1648435/1‑150 (MQ=255)
                                                                    ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCCTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                                  <  2:1653366/150‑1 (MQ=255)
                                                                     cGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgtg                                                                                 <  1:1988793/150‑1 (MQ=255)
                                                                         acaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAc                                                                             >  2:952215/1‑150 (MQ=255)
                                                                         acaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAc                                                                             >  1:1407335/1‑150 (MQ=255)
                                                                          caaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcc                                                                            >  2:1454675/1‑150 (MQ=255)
                                                                            aaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCAGGCAACACCTATCGGTCGTGGTCAGCGTGGACAGATTGTGGCAacgc                                                                          >  1:1044036/1‑150 (MQ=255)
                                                                             aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                                         >  1:190117/1‑150 (MQ=255)
                                                                             aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                                         <  1:1200636/150‑1 (MQ=255)
                                                                              aTCCTCTTTGAGAACTTAACCCCGCCGCTCGCAATCTCTGGTCTCCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGAGCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcctccg                                                                        <  2:292883/150‑1 (MQ=255)
                                                                                ccTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaa                                                                      <  2:1140623/150‑1 (MQ=255)
                                                                                 ctctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaaa                                                                     >  1:1599780/1‑150 (MQ=255)
                                                                                 ctctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaaa                                                                     <  2:1059702/150‑1 (MQ=255)
                                                                                   ctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGc                                                                   >  2:740119/1‑150 (MQ=255)
                                                                                     ttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg                                                                 <  1:1486048/150‑1 (MQ=255)
                                                                                      tGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                                <  1:770080/150‑1 (MQ=255)
                                                                                       gagaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAATATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGTTCAGCGTTGTCTGATTGTGGCACCGCCGAAAGCCGGt                                                               <  1:1475652/150‑1 (MQ=255)
                                                                                        agaACTTCACCCCGCTGCACGCACACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa                                                              <  1:92041/150‑1 (MQ=255)
                                                                                            cTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAc                                                          >  1:1222463/1‑150 (MQ=255)
                                                                                               aaCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGTCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAt                                                       <  1:1597675/150‑1 (MQ=255)
                                                                                                 ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTTTCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgg                                                     >  1:440046/1‑149 (MQ=255)
                                                                                                     gCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgct                                                 <  2:649354/150‑1 (MQ=255)
                                                                                                     gCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCGCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgct                                                 <  1:1114502/150‑1 (MQ=255)
                                                                                                       tGCACGCAATCTCTCGTCTGCGTATGGAACGGGGTAACGGTTCTTCTGAAGATTTAAGTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTTGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctgc                                               <  2:770295/150‑1 (MQ=255)
                                                                                                        gCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCa                                              <  2:816647/150‑1 (MQ=255)
                                                                                                          aCGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGa                                            >  1:1567305/1‑150 (MQ=255)
                                                                                                             cAAACTCTCGTCTGCGAATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTAGCGTACTGGATCTGGCATCACCTATCGGTATTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCATAACa                                         >  1:264852/1‑150 (MQ=255)
                                                                                                                aCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCATCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTTCAGAACAttt                                      >  1:60154/1‑149 (MQ=255)
                                                                                                                 ctctCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGc                                     >  2:1969532/1‑150 (MQ=255)
                                                                                                                 ctctCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGc                                     <  1:692680/150‑1 (MQ=255)
                                                                                                                     cGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCag                                 >  2:106636/1‑150 (MQ=255)
                                                                                                                       tCTGCGTATGGAACGTGGTTACTGTTCTACTTAAGATTTATGTGCTCTCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCagag                               <  1:191352/150‑1 (MQ=255)
                                                                                                                         tGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa                             >  2:57221/1‑150 (MQ=255)
                                                                                                                             tATGGAACGTGGTAACGGTTCTAATGAAGATTAAACTGCCCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGc                         <  2:949592/150‑1 (MQ=255)
                                                                                                                              aTGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt                        >  1:1644013/1‑150 (MQ=255)
                                                                                                                              aTGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt                        <  2:1762395/150‑1 (MQ=255)
                                                                                                                              aTGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGGCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt                        <  1:87562/150‑1 (MQ=255)
                                                                                                                               tGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt                       <  2:447719/150‑1 (MQ=255)
                                                                                                                                ggaaaatGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGACGAAAGCCGGTAATACCATGCTTCTGCAGAACATTGCTCAGAGCATTGCTTa                      <  2:1250115/144‑1 (MQ=255)
                                                                                                                                    cattgTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGTCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAac                  <  1:158980/146‑1 (MQ=255)
                                                                                                                                      tGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                                <  1:1092348/102‑1 (MQ=255)
                                                                                                                                      tGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                                >  2:1092348/1‑102 (MQ=255)
                                                                                                                                       ggTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGGTCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAACGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAaccac               >  2:1254135/1‑150 (MQ=255)
                                                                                                                                            cGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACACCCACCCggt          >  2:64208/1‑149 (MQ=255)
                                                                                                                                              gTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAtt        >  1:1243304/1‑150 (MQ=255)
                                                                                                                                               ttCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtg       <  1:1672899/150‑1 (MQ=255)
                                                                                                                                               ttCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtg       <  1:1502026/150‑1 (MQ=255)
                                                                                                                                                 cTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtg     <  2:810646/150‑1 (MQ=255)
                                                                                                                                                 cTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGAGTGGTCTGATTGTGGCACAGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTAAAACCACCCGGATtgtg     >  1:184023/1‑150 (MQ=255)
                                                                                                                                                  tACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATtgtgt    >  1:2062662/1‑150 (MQ=255)
                                                                                                                                                    cTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGc  <  1:588934/150‑1 (MQ=255)
                                                                                                                                                    cTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGc  >  2:1342612/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGC  >  NC_000913/3966729‑3967026

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