Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,168,420 A→T 6.9% intergenic (‑114/‑292) yegR ← / → yegS uncharacterized protein/phosphatidylglycerol kinase, metal‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,168,4200AT 6.9% 299.3 / 11.5 87intergenic (‑114/‑292)yegR/yegSuncharacterized protein/phosphatidylglycerol kinase, metal‑dependent
Reads supporting (aligned to +/- strand):  ref base A (37/44);  new base T (3/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

ATTAATGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTC  >  NC_000913/2168278‑2168567
                                                                                                                                              |                                                                                                                                                   
attaatGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTa                                                                                                                                               <  3:184368/149‑1 (MQ=255)
attaatGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTa                                                                                                                                               <  3:201288/149‑1 (MQ=255)
 ttaatGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAg                                                                                                                                              <  2:630481/149‑1 (MQ=255)
    atGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGaa                                                                                                                                           <  1:758768/149‑1 (MQ=255)
     tGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAAt                                                                                                                                          <  2:435181/149‑1 (MQ=255)
      gATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACTTATGTAGGaaat                                                                                                                                         <  1:109633/149‑3 (MQ=255)
       atatAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACTTATGTAGGaaatg                                                                                                                                        <  2:272920/149‑4 (MQ=255)
        tataGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTAc                                                                                                                                       <  4:484495/149‑1 (MQ=255)
           agcagcTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGa                                                                                                                                    <  4:47197/149‑1 (MQ=255)
              agcTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACAt                                                                                                                                 >  3:568648/1‑149 (MQ=255)
              agcTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAAGAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACag                                                                                                                                 <  2:690502/149‑2 (MQ=255)
                 tATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAg                                                                                                                                              <  3:476869/133‑1 (MQ=255)
                 tATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAg                                                                                                                                              >  4:476869/1‑133 (MQ=255)
                      ttttCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGaa                                                                                                                                           >  3:5951/1‑131 (MQ=255)
                      ttttCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGaa                                                                                                                                           <  4:5951/131‑1 (MQ=255)
                            tctcACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAgg                                                                                                                   >  3:763745/1‑149 (MQ=255)
                             ctcACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTAAGGACATTGAGGCCCTTCAggg                                                                                                                  >  2:652038/1‑149 (MQ=255)
                               cACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGACATGATTATTATAATTGAATAGGATTATAATATATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTa                                                                                                                <  4:763745/149‑1 (MQ=255)
                                aCCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGa                                                                                                                                    >  3:97479/1‑128 (MQ=255)
                                aCCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGa                                                                                                                                    <  4:97479/128‑1 (MQ=255)
                                aCCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTaa                                                                                                               >  4:202322/1‑149 (MQ=255)
                                       cccATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTaa                                                                                                               >  4:420717/1‑142 (MQ=255)
                                       cccATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTaa                                                                                                               <  3:420717/142‑1 (MQ=255)
                                       cccATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATg                                                                                                        >  2:498351/1‑149 (MQ=255)
                                        ccATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATgg                                                                                                       >  1:442075/1‑149 (MQ=255)
                                          aTTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATgggg                                                                                                     >  2:81353/1‑149 (MQ=255)
                                           ttAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCAttgtgg                                                                                                    >  3:767174/1‑144 (MQ=255)
                                            tAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGc                                                                                                   <  2:224560/149‑1 (MQ=255)
                                                 aCCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACTTATGTAGGaa                                                                                                                                           <  4:663429/104‑1 (MQ=255)
                                                 aCCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACTTATGTAGGaa                                                                                                                                           >  3:663429/1‑104 (MQ=255)
                                                 aCCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTaa                                                                                              <  1:652038/149‑1 (MQ=255)
                                                      cGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAatatat                                                                                         <  1:600566/149‑1 (MQ=255)
                                                            tctctAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATatattatatt                                                                                   <  3:202322/149‑1 (MQ=255)
                                                                tAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTAc                                                                                                                                       <  4:277826/93‑1 (MQ=255)
                                                                tAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTAc                                                                                                                                       >  3:277826/1‑93 (MQ=255)
                                                                    gtTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGc                                                                           <  1:498351/149‑1 (MQ=255)
                                                                     tttAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAAGAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATg                                                                                                        >  4:653573/1‑119 (MQ=255)
                                                                     cttAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATg                                                                                                        <  3:653573/118‑1 (MQ=255)
                                                                      ttAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCa                                                                         >  4:80037/1‑149 (MQ=255)
                                                                        aaaTATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACt                                                                       <  4:217272/149‑1 (MQ=255)
                                                                           tataATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGtt                                                                    <  4:568648/149‑1 (MQ=255)
                                                                             taataaTGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAgg                                                                                                                   <  3:306962/100‑1 (MQ=255)
                                                                             taataaTGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAgg                                                                                                                   >  4:306962/1‑100 (MQ=255)
                                                                               ataatGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTAGGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGt                                                                >  3:543609/1‑149 (MQ=255)
                                                                                taatGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGtt                                                               <  1:403899/149‑1 (MQ=255)
                                                                                taatGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGtt                                                               <  1:81353/149‑1 (MQ=255)
                                                                                 aatGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTg                                                              >  1:293481/1‑149 (MQ=255)
                                                                                  atGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGAGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGa                                                             <  3:380003/149‑1 (MQ=255)
                                                                                        tGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTa                                                       <  2:442075/149‑1 (MQ=255)
                                                                                         gATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAggg                                                                                                                  <  3:757821/89‑1 (MQ=255)
                                                                                         gATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAggg                                                                                                                  >  4:757821/1‑89 (MQ=255)
                                                                                            tattatAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGta                                                   <  4:543609/149‑1 (MQ=255)
                                                                                             attatAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATg                                                                                                        <  4:57619/95‑1 (MQ=255)
                                                                                             attatAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGCAACTCCATg                                                                                                        >  3:57619/1‑95 (MQ=255)
                                                                                               tataATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTatat                                                >  2:631825/1‑149 (MQ=255)
                                                                                                    ttGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAatatat                                                                                         <  3:17440/103‑1 (MQ=255)
                                                                                                    ttGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAatatat                                                                                         >  4:17440/1‑103 (MQ=255)
                                                                                                        aTAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGtt                                                               >  2:82573/1‑125 (MQ=255)
                                                                                                        aTAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGtt                                                               <  1:82573/125‑1 (MQ=255)
                                                                                                         tAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCGATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCtttttt                                      <  4:767174/149‑1 (MQ=255)
                                                                                                          aGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTg                                     <  1:631825/149‑1 (MQ=255)
                                                                                                            gATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAgg                                                                                                                   >  4:9512/1‑69 (MQ=255)
                                                                                                            gATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAgg                                                                                                                   <  3:9512/69‑1 (MQ=255)
                                                                                                                 taataaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGaaa                                  <  4:407822/145‑1 (MQ=255)
                                                                                                                 taataaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGaaa                                  >  3:407822/1‑145 (MQ=255)
                                                                                                                  aataaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTa                                                       >  1:561040/1‑123 (MQ=255)
                                                                                                                  aataaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTa                                                       <  2:561040/123‑1 (MQ=255)
                                                                                                                  aataaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTAt                             >  4:90029/1‑149 (MQ=255)
                                                                                                                   ataaATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATc                            >  2:222747/1‑149 (MQ=255)
                                                                                                                      aaTGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAAt                                 >  1:22382/1‑141 (MQ=255)
                                                                                                                      aaTGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAAt                                 <  2:22382/141‑1 (MQ=255)
                                                                                                                          ttCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCa                     <  3:90029/149‑1 (MQ=255)
                                                                                                                               tACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGt                >  1:141251/1‑149 (MQ=255)
                                                                                                                                 caacaTTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGtata                                                 >  2:371426/1‑114 (MQ=255)
                                                                                                                                 caacaTTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGtata                                                 <  1:371426/114‑1 (MQ=255)
                                                                                                                                  aacaTTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGc                                                                           >  2:713175/1‑87 (MQ=255)
                                                                                                                                  aacaTTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGc                                                                           <  1:713175/87‑1 (MQ=255)
                                                                                                                                   gtaattCCTACTTATGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGt            >  2:589170/5‑149 (MQ=255)
                                                                                                                                      attCCTACTTATGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTAtt         >  4:629570/2‑149 (MQ=255)
                                                                                                                                          cTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAat                                                           >  3:305876/1‑95 (MQ=255)
                                                                                                                                          cTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAat                                                           <  4:305876/95‑1 (MQ=255)
                                                                                                                                          cTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGACTATTAg                                                      >  3:525493/1‑100 (MQ=255)
                                                                                                                                          cTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAg                                                      <  4:525493/100‑1 (MQ=255)
                                                                                                                                          cTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCt     >  3:768478/1‑149 (MQ=255)
                                                                                                                                            aCATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCt     >  3:135499/1‑147 (MQ=255)
                                                                                                                                            aCATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCt     <  4:135499/147‑1 (MQ=255)
                                                                                                                                             cATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTc  >  2:11007/1‑149 (MQ=255)
                                                                                                                                              aTAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGaaa                                  >  3:140913/1‑116 (MQ=255)
                                                                                                                                              aTAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGaaa                                  <  4:140913/116‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ATTAATGATATAGCAGCTATCTTTTTCATCTCACCTTCCCCCATTAAATACCAACGACACTCTCTAGTGTTTAAATATAATAATGGCATGATTATTATAATTGAATAGGATTATAATAAATGTTCTGTACAACATTTCCTACATAAGTAGGAATTACGGACATTGAGGCCCTTCAGGGTAACTCCATGGGGGCTTTAATATATTATATTGAAGATGCCACTGTTTAGTTGAATATTAGGTATATGCTCTTTTTTGAAATTTATCGGTGGCAGCCGTTAGTATTCGCTGTC  >  NC_000913/2168278‑2168567

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