Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,122,658 G→C 5.0% P228A (CCA→GCA)  ftsN ← essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,122,6580GC 5.0% 270.6 / 6.3 82P228A (CCA→GCA) ftsNessential cell division protein
Reads supporting (aligned to +/- strand):  ref base G (31/45);  new base C (2/2);  total (34/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-02

ATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAG  >  NC_000913/4122510‑4122806
                                                                                                                                                    |                                                                                                                                                    
aTTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTgg                                                                                                                                                      <  1:563209/149‑1 (MQ=255)
    tCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCa                                                                                                                                                  <  2:407627/149‑1 (MQ=255)
    tCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCa                                                                                                                                                  <  1:98663/149‑1 (MQ=255)
      gAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGCCGCAGc                                                                                                                                                >  2:708599/1‑149 (MQ=255)
           aaaGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTtgct                                                                                                                                           >  2:819877/1‑149 (MQ=255)
           aaaGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTtgct                                                                                                                                           >  4:10117/1‑149 (MQ=255)
             aGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTtgctgc                                                                                                                                         >  4:759811/1‑149 (MQ=255)
              gCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCg                                                                                                                                        >  2:446342/1‑149 (MQ=255)
                cTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCCACTGGCGCAGCTTGCTGCGGc                                                                                                                                      <  4:861060/149‑1 (MQ=255)
                cTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGc                                                                                                                                      <  2:628805/149‑1 (MQ=255)
                  tCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCtt                                                                                                                                    <  3:683835/149‑1 (MQ=255)
                   cGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTc                                                                                                                                   <  2:84447/149‑1 (MQ=255)
                   cGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGCCGCAGCAAGCTGCgccagt                                                                                                                                   <  1:708599/149‑6 (MQ=255)
                        gCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTg                                                                                                                              >  3:845612/1‑149 (MQ=255)
                         ccAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTgc                                                                                                                             <  1:381072/149‑1 (MQ=255)
                         ccAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTgc                                                                                                                             <  3:583668/149‑1 (MQ=255)
                         ccAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTgc                                                                                                                             >  1:795377/1‑149 (MQ=255)
                                 gCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCgt                                                                                                                     <  3:404876/149‑1 (MQ=255)
                                   aggTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCAGCATCCAGCGGCGTTCGTCATTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGCCTGCGCAGTCgtgt                                                                                                                   <  2:613113/147‑1 (MQ=255)
                                     gTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTgc                                                                                                                 >  3:450982/1‑149 (MQ=255)
                                       aCCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTgcgc                                                                                                               >  3:255902/1‑149 (MQ=255)
                                           tctcTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGa                                                                                                           <  3:860647/149‑1 (MQ=255)
                                           tctcTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGa                                                                                                           >  2:790741/1‑149 (MQ=255)
                                           tctcTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGCCGCAGCAAGCTGCgccag                                                                                                                                    <  1:354694/124‑5 (MQ=255)
                                           tctcTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGCCGCAGCAAGCTGCgccag                                                                                                                                    >  2:354694/1‑120 (MQ=255)
                                             tctGCCTTTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGt                                                                                                         >  2:471522/1‑149 (MQ=255)
                                              ctGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGtt                                                                                                        <  1:446342/149‑1 (MQ=255)
                                              ctGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCACACGACTGCGCAgtcgaa                                                                                                        <  3:702829/149‑6 (MQ=255)
                                               tGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGttt                                                                                                       <  2:453077/149‑1 (MQ=255)
                                               tGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGttt                                                                                                       >  4:326559/1‑149 (MQ=255)
                                                gCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGACAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTg                                                                                                      <  3:846955/149‑1 (MQ=255)
                                                 ccTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTgc                                                                                                     <  4:36417/149‑1 (MQ=255)
                                                  cTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTgca                                                                                                    <  1:259493/149‑1 (MQ=255)
                                                       ccGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGACACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGa                                                                                               <  3:254960/149‑1 (MQ=255)
                                                        cgcgCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGAt                                                                                              <  1:165550/149‑1 (MQ=255)
                                                          cgcCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATcc                                                                                            <  4:463431/149‑1 (MQ=255)
                                                            ccTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTg                                                                                          <  1:114160/149‑1 (MQ=255)
                                                            ccTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTg                                                                                          <  1:222621/149‑1 (MQ=255)
                                                             ctctGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTgg                                                                                         >  4:621451/1‑149 (MQ=255)
                                                              tctGAACGAACCGCGCTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGt                                                                                        <  1:42561/149‑1 (MQ=255)
                                                              tctGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTc                                                                                                                                   >  4:186961/1‑106 (MQ=255)
                                                              tctGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTc                                                                                                                                   <  3:186961/106‑1 (MQ=255)
                                                              tctGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGt                                                                                        >  3:619888/1‑149 (MQ=255)
                                                                tGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTAc                                                                                      <  3:350103/149‑1 (MQ=255)
                                                                  aacgaacCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACgg                                                                                    <  4:619888/149‑1 (MQ=255)
                                                                  aacgaacCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGACGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTg                                                                                          <  1:115306/143‑1 (MQ=255)
                                                                  aacgaacCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGACGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTg                                                                                          >  2:115306/1‑143 (MQ=255)
                                                                   acgaacCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGgc                                                                                   >  4:563715/1‑149 (MQ=255)
                                                                       acCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGgctgct                                                                               <  2:640198/149‑1 (MQ=255)
                                                                       acCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGGCTGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGCTCCTGGTACGggtgct                                                                               <  2:822243/149‑1 (MQ=255)
                                                                        ccGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGCGTTTGCAGCAGATCCTGGTACGGctgctg                                                                              >  2:835236/1‑149 (MQ=255)
                                                                         cGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTga                                                                             >  1:124607/1‑149 (MQ=255)
                                                                         cGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTga                                                                             <  1:395257/149‑1 (MQ=255)
                                                                         cGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTga                                                                             <  2:584133/149‑1 (MQ=255)
                                                                            aCTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGagaa                                                                          >  3:352860/1‑149 (MQ=255)
                                                                               gCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAgaagaa                                                                       >  2:270621/1‑149 (MQ=255)
                                                                                  ccATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGcc                                                                    >  2:52087/1‑149 (MQ=255)
                                                                                  ccATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTAGCGCAGCCTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGcc                                                                    <  4:566120/149‑1 (MQ=255)
                                                                                   cATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCg                                                                   <  4:618122/149‑1 (MQ=255)
                                                                                       cAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGttt                                                               <  1:554030/149‑1 (MQ=255)
                                                                                        aGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTa                                                              <  3:759811/149‑1 (MQ=255)
                                                                                          cggcgTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGa                                                            <  1:724023/149‑1 (MQ=255)
                                                                                           ggcgTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGa                                                            >  1:45659/1‑148 (MQ=255)
                                                                                           ggcgTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGa                                                            <  2:45659/148‑1 (MQ=255)
                                                                                           ggcgTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGAc                                                           <  1:471522/149‑1 (MQ=255)
                                                                                               ttCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGgc                                                       >  4:794213/1‑149 (MQ=255)
                                                                                               ttCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGgc                                                       >  4:74572/1‑149 (MQ=255)
                                                                                               ttCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGgc                                                       <  3:182415/149‑1 (MQ=255)
                                                                                                 cGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAAACTGgcgc                                                     <  4:75366/149‑1 (MQ=255)
                                                                                                  gTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGcgcg                                                    >  2:510901/1‑149 (MQ=255)
                                                                                                  gTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGcgcg                                                    <  1:74840/149‑1 (MQ=255)
                                                                                                   tCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCgg                                                   >  1:838737/1‑149 (MQ=255)
                                                                                                    cTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGt                                                  <  2:823887/149‑1 (MQ=255)
                                                                                                       ttttCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTgc                                               <  3:69480/149‑1 (MQ=255)
                                                                                                            tcCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGTCGCGGTTGCGCCTg                                          >  4:386989/1‑149 (MQ=255)
                                                                                                                       ttttGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCgg                               >  3:788254/1‑149 (MQ=255)
                                                                                                                         ttGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCt                             <  4:112305/149‑1 (MQ=255)
                                                                                                                              gcgTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGa                        >  3:76734/1‑149 (MQ=255)
                                                                                                                                 tCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGt                                                                 >  1:381251/1‑105 (MQ=255)
                                                                                                                                 tCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCATGGTACGGCTGCTGAGAAGAAGCCGGt                                                                 <  2:381251/105‑1 (MQ=255)
                                                                                                                                  cAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTa                    <  1:103563/149‑1 (MQ=255)
                                                                                                                                    gCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTAcg                  >  3:110819/1‑149 (MQ=255)
                                                                                                                                               aaCTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTACGCGTctgctgct       >  3:404878/1‑149 (MQ=255)
                                                                                                                                                    gcgcAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAg  >  3:480987/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTTTCGAGTCAAAGCCTTCGAACGCCAGCTGAGCACGTACCGTCTCTGCCTGTTCCGCGCCTCTGAACGAACCGCACTGCACCATCCAGCGGCGTTCGTCTTTTTTCTCCGCCGTCGGTTTTGGCGCGTCAGCGGCACGAGCAACTGGCGCAGCTTGCTGCGGCTTCGACTGCGCAGTCGTGTGCGCAGGAGTTTGCAGCAGATCCTGGTACGGCTGCTGAGAAGAAGCCGGTTTAGACTGGCGCGGTTGCGCCTGTACCGGTGCGGCTTGCGACGTACGCGTCTGCTGCTGCCAG  >  NC_000913/4122510‑4122806

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