Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,722,301 A→C 5.0% intergenic (+196/+27) sokA → / ← glyS pseudogene, sRNA antisense regulator blocking mokA/hokA translation/glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,722,3010AC 5.0% 693.7 / 23.4 196intergenic (+196/+27)sokA/glySpseudogene, sRNA antisense regulator blocking mokA/hokA translation/glycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (110/76);  new base C (6/4);  total (116/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 197 total aligned reads displayed.

CCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAA  >  NC_000913/3722153‑3722448
                                                                                                                                                    |                                                                                                                                                   
ccGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGa                                                                                                                                                     <  1:604152/149‑1 (MQ=255)
  gTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATg                                                                                                                                                   <  4:519430/149‑1 (MQ=255)
              gCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCtttt                                                                                                                                         >  3:633751/1‑144 (MQ=255)
               cTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                                                                      <  4:506251/149‑1 (MQ=255)
                    taaGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACg                                                                                                                                 >  3:582363/1‑149 (MQ=255)
                     aaGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACgg                                                                                                                                <  2:828385/149‑1 (MQ=255)
                        tttgtttCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAAc                                                                                                                                  >  3:622957/1‑144 (MQ=255)
                              tcttctCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTtattat                                                                                                                       >  4:763655/1‑149 (MQ=255)
                                   tcCTTGATATAATCCTAAAAAAATCTAATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCATGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCaa                                                                                                                  >  4:974293/1‑149 (MQ=255)
                                     cTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGc                                                                                                                                                 <  2:772358/116‑1 (MQ=255)
                                     cTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGc                                                                                                                                                 >  1:772358/1‑116 (MQ=255)
                                     cTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCa                                                                                                                                                <  3:213384/117‑1 (MQ=255)
                                     cTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatt                                                                                                                         <  3:16493/140‑1 (MQ=255)
                                     cTTGATATAATCCTAAAAAAATCTCATAAAATTAACATATGAGATAATCTTTATTCCGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                >  4:603722/1‑149 (MQ=255)
                                               tCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                                                                         >  3:828129/1‑114 (MQ=255)
                                                ccTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATcc                                                                                                     >  1:286061/1‑149 (MQ=255)
                                                 cTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCg                                                                                                    >  1:1003291/1‑149 (MQ=255)
                                                     aaaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAc                                                                                                >  4:459246/1‑149 (MQ=255)
                                                      aaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                >  2:548148/1‑132 (MQ=255)
                                                      aaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACa                                                                                                                <  1:548148/132‑1 (MQ=255)
                                                          tctcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCtt                                                                                                                                           >  4:625214/1‑100 (MQ=255)
                                                          tctcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTtattat                                                                                                                       <  2:263853/121‑1 (MQ=255)
                                                          tctcATAAAATTAATATATGAGATAATCTTTATTCAACAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCa                                                                                           <  3:763655/149‑1 (MQ=255)
                                                            tcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGc                                                                                                              >  3:799256/1‑128 (MQ=255)
                                                            tcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAgg                                                                                         >  3:505514/1‑149 (MQ=255)
                                                              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                                                                         <  2:1127566/99‑1 (MQ=255)
                                                                 aaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGAt                                                                                                                                                    >  4:398574/1‑85 (MQ=255)
                                                                  aaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGt                                                                                   <  3:1061630/149‑1 (MQ=255)
                                                                    ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatt                                                                                                                         >  1:839638/1‑109 (MQ=255)
                                                                    ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAAc                                                                                      >  1:213704/1‑144 (MQ=255)
                                                                    ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTc                                                                                 >  3:547884/1‑149 (MQ=255)
                                                                       atatatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATc                                                                                                      >  1:751395/1‑125 (MQ=255)
                                                                               gaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtct                                                                      >  1:397099/1‑149 (MQ=255)
                                                                                aTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAata                                                                                                        >  2:752853/1‑114 (MQ=255)
                                                                                aTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtctc                                                                     <  3:596829/149‑1 (MQ=255)
                                                                                  aaTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCa                                                                   <  3:974293/149‑1 (MQ=255)
                                                                                  aaTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCa                                                                   <  4:332208/149‑1 (MQ=255)
                                                                                   aTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAg                                                                  >  4:1144996/1‑149 (MQ=255)
                                                                                    tCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGc                                                                 >  3:327130/1‑149 (MQ=255)
                                                                                    tCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGc                                                                 >  4:886069/1‑149 (MQ=255)
                                                                                      tttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAcgcg                                                                                             <  4:579962/119‑1 (MQ=255)
                                                                                      tttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTcg                                                                                >  4:649827/1‑132 (MQ=255)
                                                                                      tttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCAt                                                               >  3:1136021/1‑149 (MQ=255)
                                                                                         aTTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTcgc                                                                               <  2:16898/130‑1 (MQ=255)
                                                                                         aTTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTcgc                                                                               >  1:16898/1‑130 (MQ=255)
                                                                                             atcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAACATGGTCAGa                                                        >  4:220581/3‑149 (MQ=255)
                                                                                             agcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAg                                                                                                                                               >  2:1121564/1‑62 (MQ=255)
                                                                                             agcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAg                                                                                                                                               <  1:1121564/62‑1 (MQ=255)
                                                                                             agcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATg                                                              >  4:792272/1‑143 (MQ=255)
                                                                                             agcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc                                                       >  2:149302/1‑149 (MQ=255)
                                                                                                agaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACgg                                                     >  1:180065/1‑149 (MQ=255)
                                                                                                agaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACgg                                                     >  1:15105/1‑149 (MQ=255)
                                                                                                 gaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtc                                                                       <  4:744352/130‑1 (MQ=255)
                                                                                                 gaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGt                                                    <  2:1003291/149‑1 (MQ=255)
                                                                                                  aagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTattatt                                                                                                                      >  4:1078542/1‑82 (MQ=255)
                                                                                                    gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGa                                                 >  4:716835/1‑149 (MQ=255)
                                                                                                    gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGa                                                 >  2:457550/1‑149 (MQ=255)
                                                                                                         ttaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACgcgc                                                                                            >  2:612078/1‑101 (MQ=255)
                                                                                                         ttaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGtt                                                                                  >  4:212671/1‑111 (MQ=255)
                                                                                                         ttaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGc                                            >  4:232908/1‑149 (MQ=255)
                                                                                                          taAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCa                                           >  1:933693/1‑149 (MQ=255)
                                                                                                           aaaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCaa                                          >  2:308075/1‑149 (MQ=255)
                                                                                                            aaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATgg                                                                                                                                                  >  1:487060/1‑44 (MQ=255)
                                                                                                            aaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATc                                                                                                      <  4:92966/88‑1 (MQ=255)
                                                                                                            aaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTcc                                                                    >  1:1037920/1‑122 (MQ=255)
                                                                                                              ggTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTc                                       >  1:721475/1‑149 (MQ=255)
                                                                                                               gTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCt                                      <  4:774553/149‑1 (MQ=255)
                                                                                                                 tGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    <  3:716835/149‑1 (MQ=255)
                                                                                                                   cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGt                                  <  2:729072/149‑1 (MQ=255)
                                                                                                                    tGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtc                                 >  2:571390/1‑149 (MQ=255)
                                                                                                                       attattTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcatc                              >  2:471175/1‑149 (MQ=255)
                                                                                                                        ttattTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    >  3:3948/1‑142 (MQ=255)
                                                                                                                        ttattTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtcatca                             >  1:1150029/1‑149 (MQ=255)
                                                                                                                          attTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGt                                                    >  1:519537/1‑124 (MQ=255)
                                                                                                                          attTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGt                                                    <  2:519537/124‑1 (MQ=255)
                                                                                                                            ttAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCa                         >  2:952239/1‑149 (MQ=255)
                                                                                                                                 gATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGtc                                 >  3:231139/1‑136 (MQ=255)
                                                                                                                                  aTAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGa                      <  2:136532/146‑1 (MQ=255)
                                                                                                                                   attaaataGCCTGCCATCTGGCAGGCTTTTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCAt                  >  4:815970/8‑149 (MQ=255)
                                                                                                                                   attaaataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCAt                  >  3:789285/4‑149 (MQ=255)
                                                                                                                                    aaaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtctc                                                                     >  3:939739/1‑97 (MQ=255)
                                                                                                                                    aaaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATc                 >  3:900014/1‑149 (MQ=255)
                                                                                                                                      aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGa                                                        <  2:583994/108‑1 (MQ=255)
                                                                                                                                      aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGa                                                        <  4:195826/108‑1 (MQ=255)
                                                                                                                                      aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGa                                                        >  1:583994/1‑108 (MQ=255)
                                                                                                                                      aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCAc               >  3:996086/1‑149 (MQ=255)
                                                                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTc                                       <  4:592678/122‑1 (MQ=255)
                                                                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGt                                  <  3:104326/127‑1 (MQ=255)
                                                                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGt                                  >  4:104326/1‑127 (MQ=255)
                                                                                                                                          aGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCttt                                    >  1:34632/1‑124 (MQ=255)
                                                                                                                                             cTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCg        >  2:381525/1‑149 (MQ=255)
                                                                                                                                              tGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGa                                                        <  4:482911/100‑1 (MQ=255)
                                                                                                                                              tGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGa                                                        >  3:482911/1‑100 (MQ=255)
                                                                                                                                              tGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCtt                                     >  1:626347/1‑119 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      <  3:459246/149‑1 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      <  3:202505/149‑1 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      >  4:845475/1‑149 (MQ=255)
                                                                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaa      >  1:145975/1‑149 (MQ=255)
                                                                                                                                                  aGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCgaagaa   >  1:963619/1‑149 (MQ=255)
                                                                                                                                                   gATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGaaa  <  2:998786/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTTGATACGCAATTCTTTGTCATCAACCATGACCATCACTTTATCGAAGAAA  >  NC_000913/3722153‑3722448

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