Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 656,862 G→A 5.2% S102S (TCG→TCA)  pagP → phospholipid:lipid A palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913656,8620GA 5.2% 256.1 / 2.5 79S102S (TCG→TCA) pagPphospholipid:lipid A palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (2/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

CTGCCATCACCTGGCATGCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTA  >  NC_000913/656714‑657006
                                                                                                                                                    |                                                                                                                                                
cTGCCATCACCTGGCATGCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCg                                                                                                                                                  <  3:223122/149‑1 (MQ=255)
   ccATCACCTGGCATGCACGTTTCTCTTACGACAAAGAAAAAACCGATCGCTATACCGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTg                                                                                                                                                <  8:194755/148‑1 (MQ=255)
   ccATCACCTGGCATGCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTg                                                                                                                                                <  2:176910/148‑1 (MQ=255)
        aCCTGGCATGCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACaa                                                                                                                                          <  6:92211/149‑1 (MQ=255)
          cTGTCATGCACGTTTCGCTTACGACAAACAAAAAACGGATCGATATAACGCGCGACCGTGGGGTTGCGGTTTTGGCCTGTCGGGTTGTGATGAAAAAGGAAACTGGCATGTCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAAt                                                                                                                                        <  2:3360/149‑1 (MQ=255)
             gCATGCACGTTTCGATTACGACAAAGAAATAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTTGCCTGTGTCGTTTGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATgg                                                                                                                                      <  1:242039/148‑1 (MQ=255)
                 gCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAAc                                                                                                                                  <  5:37808/148‑1 (MQ=255)
                 gCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAAc                                                                                                                                  >  6:37808/1‑148 (MQ=255)
                  cACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAAcc                                                                                                                                 >  7:311528/1‑148 (MQ=255)
                     gTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATgg                                                                                                                                      >  7:163001/1‑140 (MQ=255)
                     gTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATgg                                                                                                                                      <  8:163001/140‑1 (MQ=255)
                     gTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGaa                                                                                                                                   <  4:294436/143‑1 (MQ=255)
                     gTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGaa                                                                                                                                   >  3:294436/1‑143 (MQ=255)
                     gTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAAc                                                                                                                                  >  5:338069/1‑144 (MQ=255)
                     gTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAAc                                                                                                                                  <  6:338069/144‑1 (MQ=255)
                      tttCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTg                                                                                                                            >  7:37392/1‑149 (MQ=255)
                         cGCTTACGAGAAAGAAAAAACCGATCGCTATAACGAGCGACGGTGGGGGGGCGGTTTTGGCCTGGCGCGTTGGGTTGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGcc                                                                                                                          <  6:232219/148‑1 (MQ=255)
                           cTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGGTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGa                                                                                                                       >  2:330462/1‑149 (MQ=255)
                               cGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATAc                                                                                                                    >  3:132610/1‑148 (MQ=255)
                               cGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATAc                                                                                                                    <  4:132610/148‑1 (MQ=255)
                                 aCAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGa                                                                                                                 <  8:271091/149‑1 (MQ=255)
                                   aaaGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGcc                                                                                                                          >  2:241886/1‑138 (MQ=255)
                                   aaaGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGcc                                                                                                                          <  1:241886/138‑1 (MQ=255)
                                       aaaaaaCCGATCGCTATAACGAGCGACAGTGGGGTGGCGGATTTGGCCTGTCGAGATGGGATCAAAAAGGAAACTGGCATGGCCTGTAAGCCATGGCAATTAAGGACTCGTGGAACAAATGGGAACCGATTGCAGGATACGGATGGGa                                                                                                            >  1:274028/1‑148 (MQ=255)
                                              cGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATACCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCt                                                                                                    >  3:224739/1‑149 (MQ=255)
                                                       taACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCg                                                                                            <  2:186530/148‑1 (MQ=255)
                                                            aGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCg                                                                                      <  8:311528/149‑1 (MQ=255)
                                                             gCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCATGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCg                                                                                      >  7:158152/1‑148 (MQ=255)
                                                             gCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCATGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCg                                                                                      >  7:152487/1‑148 (MQ=255)
                                                                   gTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGaaa                                                                                                          >  1:368002/1‑122 (MQ=255)
                                                                   gTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGaaa                                                                                                          <  2:368002/122‑1 (MQ=255)
                                                                   gTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGaaa                                                                               <  3:18927/149‑1 (MQ=255)
                                                                    tGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGaaaa                                                                              <  4:224002/149‑1 (MQ=255)
                                                                       ggTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAAttt                                                                           <  8:37392/149‑1 (MQ=255)
                                                                        gTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACAGATTGCCGGATACGGATGGGAAAGTACCTGTCGACCGATGGCTGATGAAAAtttt                                                                          >  6:303558/1‑149 (MQ=255)
                                                                          ggCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCa                                                                        >  1:360833/1‑149 (MQ=255)
                                                                           gCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGAt                                                                                                                      <  5:114583/102‑1 (MQ=255)
                                                                           gCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGAt                                                                                                                      >  6:114583/1‑102 (MQ=255)
                                                                           gCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCAt                                                                       >  7:185558/1‑149 (MQ=255)
                                                                            cGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCAtt                                                                      <  5:140247/149‑1 (MQ=255)
                                                                               ttttGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGTGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAg                                                                   <  5:303558/149‑1 (MQ=255)
                                                                               ttttGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAg                                                                   <  2:136322/149‑1 (MQ=255)
                                                                               ttttGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATACCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAg                                                                   <  4:224739/149‑1 (MQ=255)
                                                                                    gCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCCTGGAACAAATGGTAACCGATTGGGGGATACGTCTGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCt                                                               <  2:274028/148‑1 (MQ=255)
                                                                                       tGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTg                                                              >  3:221256/1‑146 (MQ=255)
                                                                                       tGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTg                                                              <  4:221256/146‑1 (MQ=255)
                                                                                         tCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGa                                                           <  2:360833/147‑1 (MQ=255)
                                                                                         tCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGAtt                                                         >  4:141954/1‑149 (MQ=255)
                                                                                           gcgTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTgg                                                                                        <  6:113284/116‑1 (MQ=255)
                                                                                           gcgTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTgg                                                                                        >  5:113284/1‑116 (MQ=255)
                                                                                                   aTGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGc                                               <  7:96392/149‑1 (MQ=255)
                                                                                                   aTGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCATGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGc                                               <  8:158152/149‑1 (MQ=255)
                                                                                                   aTGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCATGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGc                                               <  8:152487/149‑1 (MQ=255)
                                                                                                    tGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAAt                                                                             <  5:82443/118‑1 (MQ=255)
                                                                                                    tGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAAt                                                                             >  6:82443/1‑118 (MQ=255)
                                                                                                    tGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCg                                              >  5:381319/1‑149 (MQ=255)
                                                                                                     gAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCg                                              >  3:330578/1‑148 (MQ=255)
                                                                                                       aaaaGGAAACTGGCATGGCCTGTATGCCATGGCATCTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTaa                                           >  4:370467/1‑149 (MQ=255)
                                                                                                         aaGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACg                                         >  1:339864/1‑149 (MQ=255)
                                                                                                         aaGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACg                                         >  7:272753/1‑149 (MQ=255)
                                                                                                          aGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGATGGCGTACCgc                                        >  4:285145/1‑148 (MQ=255)
                                                                                                           ggAAACAGCCATCGACTCTATCCCATGCCATAAAAGCACACGACGAACAAAAAGGAACCGATTACACCAAACGGAATGCAAAGAACCTCGCGACCGCTGGCCGATGAAAATTTTCATTTAGGTCTGGGATTCACAGCTGGCCTAAccg                                        >  8:300576/1‑146 (MQ=255)
                                                                                                             aaaCTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCAc                                     >  7:394779/1‑149 (MQ=255)
                                                                                                                cTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCAcgc                                   >  4:75099/1‑148 (MQ=255)
                                                                                                                   gcctGGCCTGTAAGCCATGGCATTTAAGGCCTCGGGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGGCTGGGATTCACCGCTGGCGTAACGGCACGCGa                                 <  4:331517/144‑1 (MQ=255)
                                                                                                                        gCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTg                           <  1:330462/148‑1 (MQ=255)
                                                                                                                         ccTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGc                                               <  7:279022/127‑1 (MQ=255)
                                                                                                                         ccTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGc                                               >  8:279022/1‑127 (MQ=255)
                                                                                                                         ccTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGa                         <  7:300576/149‑1 (MQ=255)
                                                                                                                             tATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTa                     >  8:295427/1‑149 (MQ=255)
                                                                                                                                  cATGGCATTTAAGGACTCGTGGAACAAATGGGACCCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATcc                <  3:75099/149‑1 (MQ=255)
                                                                                                                                  cATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATcc                >  8:218767/1‑149 (MQ=255)
                                                                                                                                  cATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATcc                >  4:357138/1‑149 (MQ=255)
                                                                                                                                    tGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACATGGCGACCGCTGGCGGATGAAAATTTTGATTTAGGTCTGGGATTCACCGGTGGCGTAACGGCCCGCGATAACTGGAATTACATACct              >  6:298130/1‑149 (MQ=255)
                                                                                                                                    tGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGACGGATACGGATGGGAAAGTACCTGGAGAACTCTGGCGGATGAAAATTTTAATTTAGGTCTGGGATTAACCGCAGGCGTAACGGCAAGCCATAACTGGAATTAAATccc               >  5:254143/1‑148 (MQ=255)
                                                                                                                                             aGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGtt     >  5:175889/1‑149 (MQ=255)
                                                                                                                                             aGGACTCGTGGAACAAAAGGGAACCGTTTTCCGGATAAGGATGGGTAATTTCCTTGCGACCGCGTGCGGATGAAAATTTTCATTTAGGTGTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGtt     <  3:285145/149‑1 (MQ=255)
                                                                                                                                              ggACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGCCCGCTGGGGGATGAAAATTTTGATTTAGGTGTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGTAATTACATGCCTCTCCCGGtac    >  4:2433/1‑147 (MQ=255)
                                                                                                                                               gACTCGTGGAACAAATTGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCCCTGGCGGATGAAAATTTTCATTTAGGTCTTGTATTAACAGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCGCCCGGTTc    >  4:25844/1‑148 (MQ=255)
                                                                                                                                                aCTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTa  <  8:185558/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTGCCATCACCTGGCATGCACGTTTCGCTTACGACAAAGAAAAAACCGATCGCTATAACGAGCGACCGTGGGGTGGCGGTTTTGGCCTGTCGCGTTGGGATGAAAAAGGAAACTGGCATGGCCTGTATGCCATGGCATTTAAGGACTCGTGGAACAAATGGGAACCGATTGCCGGATACGGATGGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTA  >  NC_000913/656714‑657006

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