Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 745,368 Δ1 bp 6.9% coding (204/735 nt) ybgL → UPF0271 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913745,3680C. 6.9% 304.9 / 16.8 87coding (204/735 nt)ybgLUPF0271 family protein
Reads supporting (aligned to +/- strand):  ref base C (44/37);  new base . (3/3);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAAAGAAGCA  >  NC_000913/745225‑745515
                                                                                                                                               |                                                                                                                                                   
aCGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgc                                                                                                                                                >  1:1071017‑M2/1‑144 (MQ=255)
      gTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgcc                                                                                                                                               <  1:1147412‑M2/144‑7 (MQ=255)
      gTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgcc                                                                                                                                               >  2:1147412‑M2/1‑138 (MQ=255)
           ctctGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTgc                                                                                                                                     >  2:1155151/1‑149 (MQ=255)
                ccAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCa                                                                                                                                >  3:573996/1‑149 (MQ=255)
                ccAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCa                                                                                                                                >  1:579640/1‑149 (MQ=255)
                    tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGaaa                                                                                                                            >  3:277959/1‑149 (MQ=255)
                        gCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGt                                                                                                                        <  1:1155151/149‑1 (MQ=255)
                             tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAg                                                                                                                                         <  3:491980/127‑1 (MQ=255)
                             tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAg                                                                                                                                         >  4:491980/1‑127 (MQ=255)
                             tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGaaa                                                                                                                            <  4:244812/140‑1 (MQ=255)
                             tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGaaa                                                                                                                            >  3:244812/1‑140 (MQ=255)
                             tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACg                                                                                                                   >  3:1422963/1‑149 (MQ=255)
                                 tttCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCa                                                                                                               >  1:1228863/1‑149 (MQ=255)
                                      tGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGc                                                                                                          >  2:927725/1‑149 (MQ=255)
                                       gCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCt                                                                                                         <  2:1228863/149‑1 (MQ=255)
                                         aGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGt                                                                                                       <  4:573996/149‑1 (MQ=255)
                                         aGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGt                                                                                                       <  2:579640/149‑1 (MQ=255)
                                           gCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTAt                                                                                                     <  4:277959/149‑1 (MQ=255)
                                                  gcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGaaa                                                                                                                            <  2:27887/119‑1 (MQ=255)
                                                  gcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGaaa                                                                                                                            >  1:27887/1‑119 (MQ=255)
                                                      aTCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGgcgc                                                                                          <  1:927725/149‑1 (MQ=255)
                                                          tGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGt                                                                                                       <  2:252076/132‑1 (MQ=255)
                                                          tGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTg                                                                                      >  4:1161506/1‑149 (MQ=255)
                                                          tGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATACAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGt                                                                                                       >  1:252076/1‑132 (MQ=255)
                                                            cAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTgccgcc                                                                                                                                 >  2:538496/1‑104 (MQ=255)
                                                            cAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTgccgcc                                                                                                                                 <  1:538496/104‑1 (MQ=255)
                                                               gCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAAc                                                                                 >  3:631352/1‑149 (MQ=255)
                                                                 ttGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTa                                                                                                                     <  1:913272/111‑1 (MQ=255)
                                                                 ttGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTa                                                                                                                     >  2:913272/1‑111 (MQ=255)
                                                                         gtgAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCt                                                                                                                                       <  1:367293/85‑1 (MQ=255)
                                                                         gtgAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCt                                                                                                                                       >  2:367293/1‑85 (MQ=255)
                                                                         gtgAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGt                                                                       <  4:631352/149‑1 (MQ=255)
                                                                            aaGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGt                                                                                                                        <  3:449358/97‑1 (MQ=255)
                                                                            aaGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGt                                                                                                                        >  4:449358/1‑97 (MQ=255)
                                                                            aaGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCt                                                                                                         <  3:897535/112‑1 (MQ=255)
                                                                            aaGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCATTTACGCCCAGACGCt                                                                                                         >  4:897535/1‑112 (MQ=255)
                                                                                  tAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCt                                                                                                         <  3:867751/106‑1 (MQ=255)
                                                                                  tAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCt                                                                                                         >  4:867751/1‑106 (MQ=255)
                                                                                  tAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAggc                                                              >  1:608545/1‑149 (MQ=255)
                                                                                            tgtCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGc                                                    >  4:813651/1‑149 (MQ=255)
                                                                                             gtCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCt                                                                                       >  2:788586/1‑113 (MQ=255)
                                                                                             gtCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCt                                                                                       <  1:788586/113‑1 (MQ=255)
                                                                                                gcgATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCa                                                >  4:1155959/1‑149 (MQ=255)
                                                                                                 cgATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCAt                                               <  4:1422963/149‑1 (MQ=255)
                                                                                                  gATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgc                                                                                    >  4:1072575‑M2/1‑46 (MQ=255)
                                                                                                  gATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgc                                                                                    <  3:1072575‑M2/111‑66 (MQ=255)
                                                                                                  gATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATg                                              >  1:719028/1‑149 (MQ=255)
                                                                                                    ttGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatccacaggcaatattggcagaggaaaccag                                            <  3:1271520‑M2/149‑106 (MQ=255)
                                                                                                       gcgcTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGgcggcg                                                          <  4:204185/132‑1 (MQ=255)
                                                                                                       gcgcTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGgcggcg                                                          >  3:204185/1‑132 (MQ=255)
                                                                                                           tCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTTGCGCGCTGGCAACGATTGCCCGTGCGCAAGgcggcg                                                          >  1:28498/1‑128 (MQ=255)
                                                                                                           tCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGgcggcg                                                          <  2:28498/128‑1 (MQ=255)
                                                                                                             aCCCGAGTTGTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCAc                                   >  1:1256484/1‑149 (MQ=255)
                                                                                                               ccGAGTTTTCCCGACAGGGAAAATTCGGGTGgcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatc                                                                          >  2:626102‑M2/1‑32 (MQ=255)
                                                                                                               ccGAGTTTTCCCGACAGGGAAAATTCGGGTGgcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatc                                                                          <  1:909464‑M2/108‑77 (MQ=255)
                                                                                                               ccGAGTTTTCCCGACAGGGAAAATTCGGGTGgcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatc                                                                          <  1:626102‑M2/108‑77 (MQ=255)
                                                                                                               ccGAGTTTTCCCGACAGGGAAAATTCGGGTGgcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatc                                                                          >  2:909464‑M2/1‑32 (MQ=255)
                                                                                                                cGAGTTTTCCCGACAGGGAAAATTTTGGTCGCagcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatccacaggcaatattggcagaggaaaccagcgttaatagctc                                >  4:1271520‑M2/1‑32 (MQ=255)
                                                                                                                 gAGTTTTCCCGACAGGGAAAATTCGGGTGgcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatccacaggcaatattggcagaggaaaccagcgttaatag                                     <  1:995043‑M2/143‑114 (MQ=255)
                                                                                                                 gAGTTTTCCCGACAGGGAAAATTCGGGTGgcgccaatcgcgacaccattttttattgcttcacgcacgcaagcctgcatgatttgcgcatcgcctgcatgaaatccacaggcaatattggcagaggaaaccagcgttaatag                                     >  2:995043‑M2/1‑30 (MQ=255)
                                                                                                                    ttttCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCAc                                   >  4:563290/1‑142 (MQ=255)
                                                                                                                    ttttCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCAc                                   <  3:563290/142‑1 (MQ=255)
                                                                                                                    ttttCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCAtgt                            >  2:1283703/1‑149 (MQ=255)
                                                                                                                        cccGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCAtgtt                           >  1:415819/1‑146 (MQ=255)
                                                                                                                        cccGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCAtgtt                           <  2:415819/146‑1 (MQ=255)
                                                                                                                        cccGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTa                        <  3:1155959/149‑1 (MQ=255)
                                                                                                                         ccGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTAc                       >  1:610811/1‑149 (MQ=255)
                                                                                                                              aGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCa                  >  3:664912/1‑149 (MQ=255)
                                                                                                                              aGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCa                  >  1:872146/1‑149 (MQ=255)
                                                                                                                               gggAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAg                 <  2:719028/149‑1 (MQ=255)
                                                                                                                                ggAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGgc                                                                                            <  2:267755/73‑1 (MQ=255)
                                                                                                                                ggAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGgc                                                                                            >  1:267755/1‑73 (MQ=255)
                                                                                                                                ggAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAAc                    <  1:1123512/145‑1 (MQ=255)
                                                                                                                                ggAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAAc                    >  2:1123512/1‑145 (MQ=255)
                                                                                                                                 gAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATg                                                     >  1:901987/1‑111 (MQ=255)
                                                                                                                                 gAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATg                                                     <  2:901987/111‑1 (MQ=255)
                                                                                                                                 gAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAggc               >  3:847687/1‑149 (MQ=255)
                                                                                                                                   aaaTTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAggcgg             <  3:813651/149‑1 (MQ=255)
                                                                                                                                   aaaTTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAggcgg             <  4:664912/149‑1 (MQ=255)
                                                                                                                                        ttGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGaaa        <  4:847687/149‑1 (MQ=255)
                                                                                                                                          ggTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAaaga      <  2:1256484/149‑1 (MQ=255)
                                                                                                                                            tCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAaagaag    <  2:610811/149‑1 (MQ=255)
                                                                                                                                              gCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAAAGAAGca  >  2:1082738/1‑149 (MQ=255)
                                                                                                                                              gCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAAAGAAGca  >  4:793923/1‑149 (MQ=255)
                                                                                                                                              gCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAAAGAAGca  <  2:872146/149‑1 (MQ=255)
                                                                                                                                              gCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAAAGAAGca  >  2:867924/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
ACGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGACAGGGAAAATTTTGGTCGCAGCGCCATGCAGCTGCCGCCAGAAACCGTTTACGCCCAGACGCTGTATCAAATTGGCGCGCTGGCAACGATTGCCCGTGCGCAAGGCGGCGTAATGCGTCATGTCAAACCGCACGGCATGTTGTACAACCAGGCGGCGAAAGAAGCA  >  NC_000913/745225‑745515

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