Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,270,605 G→A 9.7% intergenic (+60/‑121) mepS → / → rtn murein DD‑endopeptidase, space‑maker hydrolase, mutational suppressor of prc thermosensitivity, outer membrane lipoprotein, weak murein LD‑carboxypeptidase/resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,270,6050GA 9.7% 185.3 / 9.3 62intergenic (+60/‑121)mepS/rtnmurein DD‑endopeptidase, space‑maker hydrolase, mutational suppressor of prc thermosensitivity, outer membrane lipoprotein, weak murein LD‑carboxypeptidase/resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base A (2/4);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

GTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTT  >  NC_000913/2270463‑2270749
                                                                                                                                              |                                                                                                                                                
gtggtgTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTg                                                                                                                                            >  3:54736/1‑149 (MQ=255)
  ggtgTTTTTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCt                                                                                                                                          >  2:338996/1‑149 (MQ=255)
  ggtgTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCt                                                                                                                                          <  1:10127/149‑1 (MQ=255)
  ggtgTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCt                                                                                                                                          <  2:349887/149‑1 (MQ=255)
  ggtgTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCt                                                                                                                                          <  4:68697/149‑1 (MQ=255)
   gtgtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtt                                                                                                                                         >  2:566488/1‑149 (MQ=255)
     gtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttt                                                                                                                                       <  1:210932/149‑1 (MQ=255)
         ttattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGt                                                                                                                                   >  4:607946/1‑149 (MQ=255)
          tattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGtt                                                                                                                                  <  1:338996/149‑1 (MQ=255)
            tttCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCt                                                                                                                                          >  3:669810/1‑139 (MQ=255)
            tttCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGATAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCt                                                                                                                                          <  4:669810/139‑1 (MQ=255)
                  gCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTttcattca                                                                                                                          >  3:662707/1‑149 (MQ=255)
                    atgaatgaACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatc                                                                                                                        >  2:519742/1‑149 (MQ=255)
                     tgaatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                       <  1:418346/149‑1 (MQ=255)
                         tgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAgaga                                                                                                                   >  4:184439/1‑149 (MQ=255)
                            aCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGaaaa                                                                                                                <  3:402110/149‑1 (MQ=255)
                             ccGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                                                                               <  1:481035/149‑1 (MQ=255)
                                tACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                                                                            <  4:54736/149‑1 (MQ=255)
                                tACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                                                                            <  1:397143/149‑1 (MQ=255)
                                   tGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtt                                                                                                                                         >  4:141867/1‑117 (MQ=255)
                                   tGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtt                                                                                                                                         <  3:141867/117‑1 (MQ=255)
                                     gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGtttt                                                                                                                                >  4:21298/1‑124 (MQ=255)
                                     gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGtttt                                                                                                                                <  3:21298/124‑1 (MQ=255)
                                     gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                       >  4:271677/1‑133 (MQ=255)
                                     gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                       <  3:271677/133‑1 (MQ=255)
                                         aagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCcgcg                                                                                                   >  4:77871/1‑149 (MQ=255)
                                              ttACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTg                                                                                              >  3:382081/1‑149 (MQ=255)
                                                        gCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                     >  1:666313/1‑100 (MQ=255)
                                                        gCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                     <  2:666313/100‑1 (MQ=255)
                                                        gCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCta                                                                                    >  3:185809/1‑149 (MQ=255)
                                                           cGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                     <  1:373751/97‑1 (MQ=255)
                                                           cGCCGGGGTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                     >  2:373751/1‑97 (MQ=255)
                                                             ccGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTc                                                                               >  4:422577/1‑149 (MQ=255)
                                                                     tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCAtt                                                                       <  3:607946/149‑1 (MQ=255)
                                                                        gCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGtt                                                                    <  4:662707/149‑1 (MQ=255)
                                                                                     aaCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgt                                                       >  3:402247/1‑149 (MQ=255)
                                                                                     aaCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgt                                                       >  2:443523/1‑149 (MQ=255)
                                                                                      aCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                                                                               <  3:359361/92‑1 (MQ=255)
                                                                                      aCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                                                                               >  4:359361/1‑92 (MQ=255)
                                                                                      aCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgtt                                                      <  3:363590/149‑1 (MQ=255)
                                                                                      aCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgtt                                                      <  1:519742/149‑1 (MQ=255)
                                                                                         gTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                   <  1:566488/149‑1 (MQ=255)
                                                                                         gTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                   <  4:442657/149‑1 (MQ=255)
                                                                                             ggATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTaa                                               <  3:184439/149‑1 (MQ=255)
                                                                                                  cAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCa                                          <  3:422577/149‑1 (MQ=255)
                                                                                                   aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                   <  1:52925/139‑1 (MQ=255)
                                                                                                   aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                   >  2:52925/1‑139 (MQ=255)
                                                                                                   aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAg                                         >  4:666091/1‑149 (MQ=255)
                                                                                                    aTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAgg                                        >  2:269961/1‑149 (MQ=255)
                                                                                                       ccTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGaaa                                     <  4:406753/149‑1 (MQ=255)
                                                                                                            gCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAAc                                <  4:382081/149‑1 (MQ=255)
                                                                                                               aTCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGa                             <  3:80936/149‑1 (MQ=255)
                                                                                                                  cTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATc                          <  4:402247/149‑1 (MQ=255)
                                                                                                                    gACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCAt                        >  3:151675/1‑149 (MQ=255)
                                                                                                                              cTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAcgc              >  4:81319/1‑149 (MQ=255)
                                                                                                                                   aGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTcc         >  2:18501/1‑149 (MQ=255)
                                                                                                                                       cTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACg                               <  2:188454/123‑1 (MQ=255)
                                                                                                                                       cTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACg                               >  1:188454/1‑123 (MQ=255)
                                                                                                                                        tGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAtt    <  4:185809/149‑1 (MQ=255)
                                                                                                                                          cTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAtttt  <  3:81319/149‑1 (MQ=255)
                                                                                                                                          cTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAtttt  >  2:248827/1‑149 (MQ=255)
                                                                                                                                          cTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAtttt  >  1:272914/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTT  >  NC_000913/2270463‑2270749

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