Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,676 G→T R87L (CGC→CTC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,6760GT98.3% 419.6 / ‑4.3 120R87L (CGC→CTC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (70/48);  total (71/49)
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

Only 100 of 120 total aligned reads displayed.

CTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGA  >  NC_000913/3966531‑3966819
                                                                                                                                                 |                                                                                                                                               
cTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGtt                                                                                                                                            >  1:534087/1‑151 (MQ=255)
 tGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGttt                                                                                                                                           >  2:1302874/1‑151 (MQ=255)
  gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGtt                                                                                                                                            <  1:1198025/149‑1 (MQ=255)
  gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTc                                                                                                                                          >  2:143951/1‑151 (MQ=255)
        gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTc                                                                                                                                    <  2:343740/151‑1 (MQ=255)
        gcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTc                                                                                                                                    <  2:741902/151‑1 (MQ=255)
          aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCg                                                                                                                                  >  1:886354/1‑151 (MQ=255)
          aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCg                                                                                                                                  <  1:962163/151‑1 (MQ=255)
          aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCg                                                                                                                                  <  1:1314125/151‑1 (MQ=255)
          aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCg                                                                                                                                  >  1:145394/1‑151 (MQ=255)
           cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGc                                                                                                                                 >  1:684410/1‑151 (MQ=255)
           cGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGc                                                                                                                                 >  1:323617/1‑151 (MQ=255)
              aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACt                                                                                                                              >  2:1576541/1‑151 (MQ=255)
              aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACt                                                                                                                              <  1:1393387/151‑1 (MQ=255)
               aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTg                                                                                                                             >  2:314914/1‑151 (MQ=255)
                  agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTg                                                                                                                          >  2:168036/1‑151 (MQ=255)
                   gTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTg                                                                                                                          >  2:508187/1‑150 (MQ=255)
                        gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATAcc                                                                                                                     >  2:170779/1‑150 (MQ=255)
                        gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATAcc                                                                                                                     >  2:1046581/1‑150 (MQ=255)
                        gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCa                                                                                                                    <  2:207075/151‑1 (MQ=255)
                         aaGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAt                                                                                                                   >  2:1009705/1‑151 (MQ=255)
                          aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                  >  1:1510368/1‑151 (MQ=255)
                          aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                  >  2:1057751/1‑151 (MQ=255)
                          aGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCAtc                                                                                                                  >  1:784640/1‑151 (MQ=255)
                               tCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                             >  1:272295/1‑151 (MQ=255)
                                 tttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGt                                                                                                            >  2:535990/1‑150 (MQ=255)
                                  ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                                                                                             <  2:281035/148‑1 (MQ=255)
                                  ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGt                                                                                                            <  2:645594/149‑1 (MQ=255)
                                   tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa                                                                                                          <  1:602193/150‑1 (MQ=255)
                                   tGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTaa                                                                                                          <  2:1069540/150‑1 (MQ=255)
                                       gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAtt                                                                                                      <  2:822036/150‑1 (MQ=255)
                                       gATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                                                                                                     <  2:79233/151‑1 (MQ=255)
                                          ggCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                  <  1:430737/151‑1 (MQ=255)
                                           gCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                                                                                                 <  2:1443030/151‑1 (MQ=255)
                                            cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccg                                                                                                >  1:959789/1‑151 (MQ=255)
                                               aCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                                                                                              <  2:1149745/150‑1 (MQ=255)
                                               aCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                             >  2:57283/1‑151 (MQ=255)
                                               aCTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                             >  2:251010/1‑151 (MQ=255)
                                                cTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                                                                                             >  1:357069/1‑150 (MQ=255)
                                                 tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                                                                                            >  1:1377332/1‑150 (MQ=255)
                                                 tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                                                                                            >  1:789986/1‑150 (MQ=255)
                                                 tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                                                                                            >  1:808033/1‑150 (MQ=255)
                                                   gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                           >  1:208135/1‑149 (MQ=255)
                                                   gagaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                         <  2:620966/151‑1 (MQ=255)
                                                    agaTATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                           >  1:1232504/1‑148 (MQ=255)
                                                      atatTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaaga                                                                                        >  1:1411499/1‑149 (MQ=255)
                                                         ttGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                   >  1:236515/1‑151 (MQ=255)
                                                          tGCAGGATGGTTTTGGTTTCCTCCGTCCCGCAGACACCTCCTCCCTCGCCGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCTCCGTTTCAACCCCCGCACGGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg                                                                                   <  1:1173212/150‑1 (MQ=255)
                                                            cAGCATGGATTGGTTTTCCTCCTTCCCGCAGCCAGCCCCCACCACGCGGCCCCTAATGAAATATACGTTCCCCATAGCCAATTCCTCCGTCTCAACCTCCGCACTGGTGACCCCATCTCTGGAAAGATCCGCCCGCCGAAAGAAGGTGAAc                                                                                <  1:99317/151‑1 (MQ=255)
                                                             aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg                                                                               >  2:722572/1‑151 (MQ=255)
                                                               gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCt                                                                             >  2:228603/1‑151 (MQ=255)
                                                                 tggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt                                                                           >  1:336367/1‑151 (MQ=255)
                                                                   gatTTGTTTTCCTCCGTTCCGCAGACACCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCTCCGTTTCACCCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                         <  1:361417/151‑1 (MQ=255)
                                                                      ttGTTTCCCCCCGTCCCGCAAACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                       <  1:43818/150‑1 (MQ=255)
                                                                       tGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                     >  1:1525994/1‑151 (MQ=255)
                                                                       tGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                     >  2:286259/1‑151 (MQ=255)
                                                                        ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                    <  2:893837/151‑1 (MQ=255)
                                                                        ggTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACGGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                    <  2:919891/151‑1 (MQ=255)
                                                                         gTTTCCTCCGTTCCGCAGACACCTCCTCCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                   <  1:1461665/151‑1 (MQ=255)
                                                                          tttCCTCCTTCCCGCAGACAGCTCCTCCCTCGCCGTTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                   <  1:1296323/150‑1 (MQ=255)
                                                                             cctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                               >  2:1442000/1‑151 (MQ=255)
                                                                                ccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                              <  2:501016/149‑1 (MQ=255)
                                                                                 cGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                           <  2:1216967/151‑1 (MQ=255)
                                                                                      cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                      <  2:1040574/151‑1 (MQ=255)
                                                                                      cGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                      >  2:1166717/1‑151 (MQ=255)
                                                                                       gcagctaGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTGCCCTAGCCAAGTCCTCCGTTTAAACCTCCGCACTGGTGATACCATCTCTGGAAAGATTCGCCCGCCGAAAGAAGGCGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                     >  1:873163/7‑151 (MQ=255)
                                                                                         aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                    >  1:597090/1‑150 (MQ=255)
                                                                                          gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                    >  1:57900/1‑149 (MQ=255)
                                                                                            cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                <  1:1513126/151‑1 (MQ=255)
                                                                                             aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                >  2:964375/1‑150 (MQ=255)
                                                                                             aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                                >  1:646049/1‑150 (MQ=255)
                                                                                              ccTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                              <  1:1047718/150‑1 (MQ=255)
                                                                                                   tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                         >  1:44879/1‑151 (MQ=255)
                                                                                                    cccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCTTTTCACCCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                        <  1:1614306/150‑1 (MQ=255)
                                                                                                        cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                     >  2:392478/1‑150 (MQ=255)
                                                                                                         gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                   >  1:1063703/1‑151 (MQ=255)
                                                                                                         gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCTCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                   <  2:287690/151‑1 (MQ=255)
                                                                                                         gCCGGTCCTGATGACATCTACGTTCCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                 <  1:798812/121‑1 (MQ=255)
                                                                                                         cccGTTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                   <  1:865720/150‑1 (MQ=255)
                                                                                                          ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                  >  2:1043620/1‑151 (MQ=255)
                                                                                                           cGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                 <  1:1616563/151‑1 (MQ=255)
                                                                                                           cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                 >  1:743407/1‑151 (MQ=255)
                                                                                                           cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                 >  2:176911/1‑151 (MQ=255)
                                                                                                            ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                <  2:1095826/151‑1 (MQ=255)
                                                                                                                 tgatgaCATCTACGTTTCCCCTAGCCAAATCCTCTGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                           >  1:427223/1‑151 (MQ=255)
                                                                                                                     gaCATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                       <  2:1307217/151‑1 (MQ=255)
                                                                                                                      aCATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                      <  2:600586/151‑1 (MQ=255)
                                                                                                                       cATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                     >  2:910463/1‑151 (MQ=255)
                                                                                                                       cATCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                     >  2:27299/1‑151 (MQ=255)
                                                                                                                        aTCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                     >  1:929497/1‑150 (MQ=255)
                                                                                                                        aTCTACGTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                     >  2:637536/1‑150 (MQ=255)
                                                                                                                              gTTTCCCCTAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa               >  2:268608/1‑150 (MQ=255)
                                                                                                                                      tAGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt      >  2:465068/1‑151 (MQ=255)
                                                                                                                                       aGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt       <  2:720068/149‑1 (MQ=255)
                                                                                                                                       aGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg     >  2:1046209/1‑151 (MQ=255)
                                                                                                                                       aGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg     >  2:1590849/1‑151 (MQ=255)
                                                                                                                                       aGCCAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg     >  2:470650/1‑151 (MQ=255)
                                                                                                                                       aGCAAAATCCGCCGTTCCACCCCCCCCACTGGTGATATGTTCCCTGGTAAGATCCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt      <  1:319634/150‑1 (MQ=255)
                                                                                                                                         ccAAATCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg     >  1:858356/1‑149 (MQ=255)
                                                                                                                                           aaaTCCTCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga  >  2:1017988/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGA  >  NC_000913/3966531‑3966819

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