Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,238,156 A→G 5.5% S175G (AGT→GGT)  yjbH → DUF940 family extracellular polysaccharide protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,238,1560AG 5.5% 227.3 / 2.3 73S175G (AGT→GGT) yjbHDUF940 family extracellular polysaccharide protein
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAG  >  NC_000913/4238010‑4238296
                                                                                                                                                  |                                                                                                                                            
gTGCGCGGGCTATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAg                                                                                                                                             <  8:392503/148‑1 (MQ=255)
   cgcgGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGAt                                                                                                                                         <  8:412654/149‑1 (MQ=255)
       ggATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATaaa                                                                                                                                      >  8:362329/1‑148 (MQ=255)
            tCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGATAAATCCGCTCTGTTCAGCCAGTGATAAATATTGc                                                                                                                                >  8:355100/1‑149 (MQ=255)
            tCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGc                                                                                                                                >  2:221159/1‑149 (MQ=255)
                 ggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGGTCGATTTTACGCTCGGCCTGGGCGGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATcg                                                                                                                           <  6:106615/149‑1 (MQ=255)
                   tACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATcgc                                                                                                                          >  8:351745/1‑148 (MQ=255)
                     cGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTACCGCGAt                                                                                                                       <  5:303710/149‑1 (MQ=255)
                           tGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGc                                                                                                                 >  3:364136/1‑149 (MQ=255)
                            gTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGc                                                                                                                 >  3:392473/1‑148 (MQ=255)
                             tttGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATAtt                                                                                                                                  <  8:30802/130‑1 (MQ=255)
                             tttGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATAtt                                                                                                                                  >  7:30802/1‑130 (MQ=255)
                             tttGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGc                                                                                                                 <  8:82962/147‑1 (MQ=255)
                             tttGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGc                                                                                                                 <  5:394134/147‑1 (MQ=255)
                               tGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTAc                                                                                                              <  4:392808/148‑1 (MQ=255)
                                 aTGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaaa                                                                                                           <  5:252068/149‑1 (MQ=255)
                                  tGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaaa                                                                                                           >  1:50608/1‑148 (MQ=255)
                                   gctgAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAAc                                                                                                          <  7:351745/145‑1 (MQ=255)
                                     ggAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAgg                                                                                                       >  8:16025/1‑149 (MQ=255)
                                           tCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGa                                                                                                  <  1:221159/148‑1 (MQ=255)
                                               gttgCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGAt                                                                                                 <  6:22104/145‑1 (MQ=255)
                                                ttgCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAt                                                                                             <  2:363123/148‑1 (MQ=255)
                                                ttgCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAt                                                                                             <  2:78110/148‑1 (MQ=255)
                                                ttgCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATc                                                                                            >  3:358454/1‑149 (MQ=255)
                                                     aGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACgg                                                                                       >  5:221474/1‑149 (MQ=255)
                                                         aaGCCTGGGGGCCGGTCGATTTTACGCTCGGCCGGGGCTGGGGGGATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGc                                                                                   <  2:18536/149‑1 (MQ=255)
                                                            gcTGGGGGCCGTTCGATTTTCCGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATaaa                                                                                                                                      <  6:104340/94‑1 (MQ=255)
                                                             cTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCCGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGc                                                                                   >  8:96496/1‑145 (MQ=255)
                                                              tGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAAc                                                                                                                    <  8:242472/111‑1 (MQ=255)
                                                              tGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAAc                                                                                                                    >  7:242472/1‑111 (MQ=255)
                                                                ggggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGt                                                                            <  4:30370/149‑1 (MQ=255)
                                                                    ccGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCa                                                                        <  4:392473/149‑1 (MQ=255)
                                                                     cGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAt                                                                                             <  8:228709/127‑1 (MQ=255)
                                                                     cGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAt                                                                                             >  7:228709/1‑127 (MQ=255)
                                                                     cGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCAc                                                                       <  4:358454/149‑1 (MQ=255)
                                                                      gTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACg                                                                      <  4:364136/149‑1 (MQ=255)
                                                                      gTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACg                                                                      >  3:245734/1‑149 (MQ=255)
                                                                       ttCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACgg                                                                     >  1:336688/1‑149 (MQ=255)
                                                                       ttCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACgg                                                                     <  7:16025/149‑1 (MQ=255)
                                                                           aTTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCt                                                                 >  8:171632/1‑149 (MQ=255)
                                                                                  gCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGTGATCTATCTACGGTAGCCAGATGTTCCa                                                                        <  4:194168/135‑1 (MQ=255)
                                                                                  gCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCa                                                                        >  3:194168/1‑135 (MQ=255)
                                                                                  gCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCAc                                                          >  8:367279/1‑149 (MQ=255)
                                                                                  gCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTACTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCa                                                                        >  6:132509/1‑135 (MQ=255)
                                                                                  gCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTACTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCa                                                                        <  5:132509/135‑1 (MQ=255)
                                                                                       gCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGttt                                                     <  7:171632/149‑1 (MQ=255)
                                                                                       gCCTGGGAGGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCGCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGttt                                                     <  2:282738/149‑1 (MQ=255)
                                                                                                    gTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGaa                                         >  7:21065/1‑148 (MQ=255)
                                                                                                    gTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAAt                                        <  7:362329/149‑1 (MQ=255)
                                                                                                       tttGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTCTTTGGCGGCGTGGAATAc                                      >  8:116499/1‑148 (MQ=255)
                                                                                                          gggCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCAACAAACAGGCGGGATCTATCGACGGTAGCCAGATTTTCCACGGTctt                                                                 >  8:169178/1‑116 (MQ=255)
                                                                                                               ccAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGCCGccc                             >  7:347798/1‑149 (MQ=255)
                                                                                                               ccAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGccc                             <  6:221474/149‑1 (MQ=255)
                                                                                                                cAGCGGTAATGTGAAAAATCCGCTCCGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCt                            <  7:96496/149‑1 (MQ=255)
                                                                                                                      tAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAAc                      >  5:16653/1‑149 (MQ=255)
                                                                                                                           tgAAAAATCCGCTCTGTTCGGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTg                 >  7:238454/1‑149 (MQ=255)
                                                                                                                            gAAAAATCCGCTCTGTTCAGCCGGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGccc                             <  7:402931/136‑1 (MQ=255)
                                                                                                                            gAAAAATCCGCTCTGTTCAGCCGGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGccc                             >  8:402931/1‑136 (MQ=255)
                                                                                                                            gAAAAATCCGCTCTGTTCAGCCGGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGccc                             >  6:122975/1‑136 (MQ=255)
                                                                                                                            gAAAAATCCGCTCTGTTCAGCCGGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGccc                             <  5:122975/136‑1 (MQ=255)
                                                                                                                             aaaaaTCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTgcg               >  8:38965/1‑149 (MQ=255)
                                                                                                                               aaaTCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTgcgc              >  4:35450/1‑148 (MQ=255)
                                                                                                                               aaaTCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTgcgc              <  3:35450/148‑1 (MQ=255)
                                                                                                                                  tCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTg           >  4:309619/1‑148 (MQ=255)
                                                                                                                                     gCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATAc                                      >  5:317565/1‑118 (MQ=255)
                                                                                                                                     gCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATAc                                      <  6:317565/118‑1 (MQ=255)
                                                                                                                                          gTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCATGGCAACCGCTGCGCCTGAAACTGGAg  >  4:293302/1‑149 (MQ=255)
                                                                                                                                           ttCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAg  >  2:323746/1‑148 (MQ=255)
                                                                                                                                           ttCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAg  >  2:169689/1‑148 (MQ=255)
                                                                                                                                           ttCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAg  <  1:169689/148‑1 (MQ=255)
                                                                                                                                           ttCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTACCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGa   <  2:50608/147‑1 (MQ=255)
                                                                                                                                                ccAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAAc                      >  8:264617/1‑123 (MQ=255)
                                                                                                                                                ccAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAAc                      <  7:264617/123‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAG  >  NC_000913/4238010‑4238296

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