Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 874,830 T→A 5.7% F618I (TTC→ATC)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913874,8300TA 5.7% 225.2 / 3.6 70F618I (TTC→ATC) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base A (2/2);  total (33/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 = 8.41e-01

TCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGA  >  NC_000913/874683‑874975
                                                                                                                                                   |                                                                                                                                                 
tCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGAttt                                                                                                                                                  <  3:307766/149‑1 (MQ=255)
 cTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTc                                                                                                                                                 <  7:17257/149‑1 (MQ=255)
    cGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTa                                                                                                                                              <  6:51152/149‑1 (MQ=255)
    cGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTa                                                                                                                                              <  2:294564/149‑1 (MQ=255)
     ggAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTAc                                                                                                                                             <  8:199501/149‑1 (MQ=255)
      gAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACt                                                                                                                                            <  3:158923/149‑1 (MQ=255)
           ggggATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCata                                                                                                                                       <  6:264893/149‑1 (MQ=255)
              gATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTa                                                                                                                                              >  2:338631/1‑139 (MQ=255)
              gATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTa                                                                                                                                              <  1:338631/139‑1 (MQ=255)
                tCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAAt                                                                                                                                  <  4:235131/149‑1 (MQ=255)
                  aTTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATg                                                                                                                                 <  7:83760/148‑1 (MQ=255)
                         cGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGGTATGGGTTAATTTAGACATTTTAGTCAGCTTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGccc                                                                                                                         <  4:247415/149‑1 (MQ=255)
                         cGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGccc                                                                                                                         >  3:5584/1‑149 (MQ=255)
                          gCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCt                                                                                                                        >  6:245336/1‑149 (MQ=255)
                          gCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATATCTTACTGCATATTAATGTTTCGCCCt                                                                                                                        >  4:149710/1‑149 (MQ=255)
                           cTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTa                                                                                                                                              <  1:13675/126‑1 (MQ=255)
                           cTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTa                                                                                                                                              >  2:13675/1‑126 (MQ=255)
                           cTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTg                                                                                                                       <  7:164118/149‑1 (MQ=255)
                              aaGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTgg                                                                                                                      <  5:185360/147‑1 (MQ=255)
                                 ttCATTACTATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGt                                                                                                                  >  8:8054/1‑148 (MQ=255)
                                    aTTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAg                                                                                                                   <  5:171966/144‑1 (MQ=255)
                                    aTTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAg                                                                                                                   >  6:171966/1‑144 (MQ=255)
                                    aTTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTaaa                                                                                                              >  4:196803/1‑149 (MQ=255)
                                        ccATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAAc                                                                                                             >  8:329403/1‑146 (MQ=255)
                                        ccATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAAc                                                                                                             <  7:329403/146‑1 (MQ=255)
                                         cATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAAc                                                                                                         <  1:331216/149‑1 (MQ=255)
                                          aTTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGaa                                                                                                          >  8:3416/1‑147 (MQ=255)
                                          aTTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGaa                                                                                                          <  7:3416/147‑1 (MQ=255)
                                               tGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACt                                                                                                   >  8:101174/1‑149 (MQ=255)
                                                 aaGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACt                                                                                                   <  4:26723/147‑1 (MQ=255)
                                                 aaGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGAGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACtt                                                                                                  >  3:83838/1‑148 (MQ=255)
                                                  aGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTc                                                                                                <  7:181069/149‑1 (MQ=255)
                                                           gTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGtt                                                                                       >  6:6695/1‑149 (MQ=255)
                                                            tttATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGttt                                                                                      <  1:54719/149‑1 (MQ=255)
                                                                       gCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAATGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTAATGAGCGTTTTACCACcatca                                                                           >  2:308644/1‑149 (MQ=255)
                                                                        cTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACcatcat                                                                          <  7:185518/149‑1 (MQ=255)
                                                                        cTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATATCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACcatcat                                                                          >  2:38546/1‑149 (MQ=255)
                                                                           ggTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATgaa                                                                       <  5:245336/149‑1 (MQ=255)
                                                                                 cAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGaaaa                                                                 >  4:109648/1‑149 (MQ=255)
                                                                                  aGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAAt                                                                <  1:95439/149‑1 (MQ=255)
                                                                                  aGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATATCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGaaaa                                                                 <  1:38546/148‑1 (MQ=255)
                                                                                   gATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATg                                                               >  6:275644/1‑149 (MQ=255)
                                                                                       attaAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTc                                                                                                <  6:303281/112‑1 (MQ=255)
                                                                                       attaAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTc                                                                                                >  5:303281/1‑112 (MQ=255)
                                                                                       attaAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCt                                                           >  3:309609/1‑149 (MQ=255)
                                                                                           aaaCGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAg                                                       >  3:186000/1‑149 (MQ=255)
                                                                                           aaaCGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATATCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGaa                                                        <  3:149710/148‑1 (MQ=255)
                                                                                                gCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGaa                                                  <  4:150771/149‑1 (MQ=255)
                                                                                                    tGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGc                                              >  8:92186/1‑149 (MQ=255)
                                                                                                    tGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGc                                              >  3:235361/1‑149 (MQ=255)
                                                                                                      ggTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCct                                            >  2:49190/1‑149 (MQ=255)
                                                                                                      ggTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCct                                            <  5:14040/149‑1 (MQ=255)
                                                                                                           aTTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCtctg                                       >  4:132835/1‑149 (MQ=255)
                                                                                                             tCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGtt                                     >  3:36456/1‑149 (MQ=255)
                                                                                                             tCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGtt                                     >  5:28762/1‑149 (MQ=255)
                                                                                                               aGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGtt                                     >  4:25394/1‑146 (MQ=255)
                                                                                                                 aCATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTgaga                                 <  7:116431/149‑1 (MQ=255)
                                                                                                                   aTTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATga                                                                        <  4:70174/108‑1 (MQ=255)
                                                                                                                   aTTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATga                                                                        >  3:70174/1‑108 (MQ=255)
                                                                                                                   aTTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGAt                                >  5:153936/1‑148 (MQ=255)
                                                                                                                    ttttAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAAc                                                 <  3:75351/130‑1 (MQ=255)
                                                                                                                    ttttAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAAc                                                 >  4:75351/1‑130 (MQ=255)
                                                                                                                          tCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGaaa                        >  3:224034/1‑149 (MQ=255)
                                                                                                                             gcgtgcgtCGCTGAAGGATGAATTCTTACTGCGTATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCg                     >  8:154843/1‑149 (MQ=255)
                                                                                                                                      gCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGc            <  4:309609/149‑1 (MQ=255)
                                                                                                                                      gCTGAAGGATGATTTCTTACTGCATATTAAAGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAATGAAAATGGCCTGAAGGCGAACAGCCGCTGAGATGAGATCACTGAAACCCTGagccacc            >  7:115889/1‑143 (MQ=255)
                                                                                                                                          aaGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGaa         <  7:8054/148‑1 (MQ=255)
                                                                                                                                          aaGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGaa         <  7:101174/148‑1 (MQ=255)
                                                                                                                                             gatgatTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTaa     <  3:109648/149‑1 (MQ=255)
                                                                                                                                              atgatTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAAt    <  7:92186/149‑1 (MQ=255)
                                                                                                                                                 atTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGa  <  4:5584/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAACGGCAATGGGTGAATTCAGACATTTTAGTCAGCGTGCGTCGCTGAAGGATGATTTCTTACTGCATATTAATGTTTCGCCCTGGCAGTTAAACGAACCACACTTTCATGAGCGTTTTACCACCATCATGAAAGAAAATGGCCTGAAGGCGAACAGCCTCTGTGTTGAGATCACTGAAACCGTGATCGAGCGAATTAATGA  >  NC_000913/874683‑874975

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: